CF961G-Partitions【斯特林数】

正题

题目链接:https://www.luogu.com.cn/problem/CF961G


题目大意

nnn个物品分成kkk个组,每个物品权值为wiw_iwi。一个子集SSS的权值为∣S∣∑x∈Swx|S|\sum_{x\in S}w_xSxSwx

求所有划分方法的权值和。


解题思路

考虑对于每个数wiw_iwi的贡献,可以看为同一集合内每一个数都对其进行过贡献。wiw_iwi对自身的贡献为{nk}\begin{Bmatrix}n\\k\end{Bmatrix}{nk}(也就是每一种划分方案都是)。对于j≠ij\neq ij=i的贡献,我们将wiw_iwiwjw_jwj绑定在一个组中,方案数就是{n−1k}\begin{Bmatrix}n-1\\k\end{Bmatrix}{n1k}。也就是对于每一个会被计算{nk}+(n−1){n−1k}\begin{Bmatrix}n\\k\end{Bmatrix}+(n-1)\begin{Bmatrix}n-1\\k\end{Bmatrix}{nk}+(n1){n1k}

考虑用斯特林数的通项计算(先是通项的推导)。
定义f(x)=xn,g(x)={nx}x!f(x)=x^n,g(x)=\begin{Bmatrix}n\\x\end{Bmatrix}x!f(x)=xn,g(x)={nx}x!
那么有f(x)=xn=∑k=0x(xk){nk}k!=∑k=0x(xk)g(k)f(x)=x^n=\sum_{k=0}^x\binom{x}{k}\begin{Bmatrix}n\\k\\\end{Bmatrix}k!=\sum_{k=0}^x\binom{x}{k}g(k)f(x)=xn=k=0x(kx){nk}k!=k=0x(kx)g(k)
根据二项式反演有g(k)=∑i=1k(−1)k−i(ni)f(i)g(k)=\sum_{i=1}^k(-1)^{k-i}\binom{n}{i}f(i)g(k)=i=1k(1)ki(in)f(i)
也就是{nk}k!=∑i=1k(−1)k−i(ni)ik\begin{Bmatrix}n\\k\end{Bmatrix}k!=\sum_{i=1}^k(-1)^{k-i}\binom{n}{i}i^k{nk}k!=i=1k(1)ki(in)ik
化简得{nm}=1m!∑k=0m(−1)k(mk)(m−k)n\begin{Bmatrix}n\\m\end{Bmatrix}=\frac{1}{m!}\sum_{k=0}^m(-1)^{k}\binom{m}{k}(m-k)^n{nm}=m!1k=0m(1)k(km)(mk)n

然后计算即可。


codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll XJQ=1e9+7,N=2e5+10;
ll n,m,ans,inv[N],fac[N];
ll power(ll x,ll b){ll ans=1;while(b){if(b&1)ans=ans*x%XJQ;x=x*x%XJQ;b>>=1;}return ans;
}
void init(){fac[0]=inv[0]=1;for(ll i=1;i<=m;i++){inv[i]=inv[i-1]*power(i,XJQ-2)%XJQ;fac[i]=fac[i-1]*i%XJQ;}
}
ll C(ll n,ll m)
{return fac[n]*inv[n-m]%XJQ*inv[m]%XJQ;}
ll S(ll n,ll m){ll ans=0,z=-1;for(ll i=0;i<=m;i++)z*=-1,(ans+=z*C(m,i)*power(m-i,n)%XJQ)%=XJQ;ans=ans*inv[m]%XJQ;return (ans+XJQ)%XJQ;
}
int main()
{scanf("%lld%lld",&n,&m);init();ll k1=S(n,m),k2=S(n-1,m)*(n-1)%XJQ;for(ll i=1;i<=n;i++){ll x;scanf("%lld",&x);(ans+=x*(k1+k2)%XJQ)%=XJQ;}printf("%lld",ans);
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/320802.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【高精】【快速幂】穿越丛林(ssl 2314)

穿越丛林 ssl 2314 题目大意&#xff1a; 求2n2^n2n 原题&#xff1a; 题目描述&#xff1a; ljj 是一位富有冒险心又很喜欢研究数学的孩纸&#xff0c;有一天&#xff0c;他到一个丛林冒险&#xff0c;这里的树长有像0、4、6、8、9这样形状的洞&#xff0c;他要想穿过丛…

谈谈surging引擎的tcp、http、ws协议和如何容器化部署

1、前言分布式已经成为了当前最热门的话题&#xff0c;分布式框架也百花齐放&#xff0c;群雄逐鹿。从中心化服务治理框架&#xff0c;到去中心化分布式服务框架&#xff0c;再到分布式微服务引擎&#xff0c;这都是通过技术不断积累改进而形成的结果。esb,网关&#xff0c;ngi…

2018陕西省赛K题[watermelon_planting]

题意:有一个序列a[]&#xff0c;描述的是另一个序列ans[]每个位置单位时间的增量。每个单位时间每个位置都会增加一个单位对应增量。时间总长m&#xff0c;每个单位时间包含有两种操作中的一个&#xff1a;1.询问ans[]在[l,r]区间的和&#xff1b;2.修改&#xff1a;a[]在[l,r]…

【并查集】银河英雄传说 (luogu 1196/ssl 1225)

银河英雄传说 luogu 1196 ssl 1225 题目大意&#xff1a; 有n列船&#xff0c;每列一开始有一艘船&#xff0c;可以将某一艘船所在的列所有船接到另外一列&#xff0c;然后会问某两艘船是否在一列&#xff0c;如果在那中间有多少艘船 原题&#xff1a; 题目描述 公元五八…

P4827-[国家集训队]Crash 的文明世界【树形dp,换根法,斯特林数】

正题 题目链接:https://www.luogu.com.cn/problem/P4827 题目大意 一颗nnn个点的树&#xff0c;定义dis(i,j)dis(i,j)dis(i,j)表示树上i,ji,ji,j两点的距离&#xff0c;对于每个iii求∑j1ndis(i,j)m\sum_{j1}^ndis(i,j)^mj1∑n​dis(i,j)m 解题思路 根据斯特林数的性质我们有…

Helm - Kubernetes服务编排的利器

Helm介绍在Kubernetes中部署容器云应用&#xff08;容器或微服务编排&#xff09;是一项有挑战性的工作&#xff0c;Helm就是为了简化在Kubernetes中安装部署容器云应用的一个客户端工具。通过Helm能够帮助开发者定义、安装和升级Kubernetes中的容器云应用。同时&#xff0c;也…

【DP】优美三角剖分

优美三角剖分 题目大意&#xff1a; 有一个图形&#xff0c;通过组合可以变成新的图形&#xff0c;现在要求第n个图形 原题&#xff1a; 题目描述 小X同学为了搞好和小C同学的关系&#xff0c;特意寻找了一些优美的图像作为礼物。 这是一些由无穷无尽三角形组成的极为优美…

Loj#114-k大异或和【线性基】

正题 题目链接:https://loj.ac/problem/114 题目大意 nnn个数&#xff0c;求所有非空子集的kkk小异或和。 解题思路 我们考虑处理一下线性基&#xff0c;我们知道对于did_idi​是一定在第iii位为111。 对于一个i>ji>ji>j且did_idi​的第jjj位为111时&#xff0c;我…

.NET Core微服务之基于MassTransit实现数据最终一致性(Part 1)

一、预备知识&#xff1a;数据一致性关于数据一致性的文章&#xff0c;园子里已经有很多了&#xff0c;如果你还不了解&#xff0c;那么可以通过以下的几篇文章去快速地了解了解&#xff0c;有个感性认识即可。&#xff08;1&#xff09;左正&#xff0c;《保证分布式系统数据一…

2017-2018 ACM-ICPC Nordic Collegiate Programming Contest (NCPC 2017)

D.Distinctive Character 看到样例&#xff0c;第一个反应贪心。先写了个按这一位1和0的数目多少&#xff0c;确定0还是1的东西。感觉不够真&#xff0c;又写了个尽量加到相似的比较小的串上的贪心。在和前边的那个组合一下&#xff0c;换了换顺序。。。好吧就过了13组样例。。…

【二分】买礼物的艰辛

买礼物的艰辛 题目大意&#xff1a; 有n个物品的价格&#xff0c;还有m个人&#xff0c;m个人中的每个人都要购买一定的连续物品&#xff0c;问买下n个物品&#xff0c;花钱最大的人最少花多少钱 原题&#xff1a; 题目描述 小X同学给小C同学选了N件礼物&#xff0c;决定顺…

P4449-于神之怒加强版【莫比乌斯反演】

正题 题目链接:https://www.luogu.com.cn/problem/P4449 题目大意 TTT组询问给出n,mn,mn,m求∑i1n∑j1mgcd(i,j)k\sum_{i1}^n\sum_{j1}^mgcd(i,j)^ki1∑n​j1∑m​gcd(i,j)k 解题思路 ∑i1n∑j1mgcd(i,j)k\sum_{i1}^n\sum_{j1}^mgcd(i,j)^ki1∑n​j1∑m​gcd(i,j)k ∑d1ndk∑i…

Asp.NetCoreWebApi图片上传接口(二)集成IdentityServer4授权访问(附源码)

写在前面本文地址&#xff1a;http://www.cnblogs.com/yilezhu/p/9315644.html作者&#xff1a;yilezhu上一篇关于Asp.Net Core Web Api图片上传的文章使用的是mongoDB进行图片的存储&#xff0c;文章发布后&#xff0c;张队就来了一句&#xff0c;说没有使用GridFS。的确博主只…

Codeforces40E[Number Table]

Codeforces40E[Number Table] 解法一&#xff0c;如果有一维很小&#xff0c;可以考虑状压dp之类的&#xff0c;显然非常不优秀。 解法二&#xff0c;如果n*m较小&#xff0c;可以考虑&#xff0c;设出每个位置是否为-1&#xff0c;解xor线性方程组。再bitset优化一下之类的。然…

【DP】天平问题

天平问题 解题思路&#xff1a; 有n个砝码&#xff0c;问可以称出多少种重量&#xff0c;可以在左边或者右边放&#xff0c;也可以不放 原题&#xff1a; 题目描述 小C为了试验小X&#xff0c;便为物竞的小X出了一道物理相关的题&#xff1a;现在给出n个质量的砝码&#x…

.NET Core开发日志——从ASP.NET Core Module到KestrelServer

ASP.NET Core程序现在变得如同控制台(Console)程序一般&#xff0c;同样通过Main方法启动整个应用。而Main方法要做的事情很简单&#xff0c;创建一个WebHostBuilder类&#xff0c;调用其Build方法生成一个WebHost类&#xff0c;最后启动之。实现代码一目了然&#xff1a;要想探…

BZOJ5358: [Lydsy1805月赛]口算训练

题解&#xff1a;判断d是否整除&#xff0c;可以转化为求这段区间内d的因子的指数是否均大于d中的指数。容易想到把每个数字都分解为素因子形式&#xff0c;对每个素数出现的次数求个前缀和即可。然而&#xff0c;这样时间空间都不行。注意到对于一个数x&#xff0c;小于sqrt(x…

jzoj6800-NOIP2020.9.19模拟spongebob【枚举】

正题 题目链接:https://gmoj.net/senior/#contest/show/3222/0 题目大意 nnn个ai,bia_i,b_iai​,bi​&#xff0c;求一个xxx使得最小化∑i1n∣aixbi∣\sum_{i1}^n|a_ixb_i|i1∑n​∣ai​xbi​∣ 解题思路 每个∣aixbi∣|a_ixb_i|∣ai​xbi​∣可以视为一个分两段的函数&#…

【bfs】极其简单的最短路问题

极其简单的最短路问题 题目大意&#xff1a; 求最短路&#xff0c;权值只有1或2 原题&#xff1a; 题目描述 小C终于被小X感动了&#xff0c;于是决定与他看电影&#xff0c;然而小X距离电影院非常远&#xff0c;现在假设每条道路需要花费小X的时间为1&#xff0c;由于有数…

GraphQL 的前世今生

GraphQL是什么GraphQL是一种新的API标准&#xff0c;它提供了一种更高效、强大和灵活的数据提供方式。它是由Facebook开发和开源&#xff0c;目前由来自世界各地的大公司和个人维护。GraphQL本质上是一种基于api的查询语言&#xff0c;现在大多数应用程序都需要从服务器中获取数…