huantwofat

#include#include#include#includeusing namespace std;int gcd(int a,int b){ while(b^=a^=b^=a%=b) ; return a;}void polya(int n,int m){ int s=0,i,j,p; for(i=1;i<=n;i++)

很简单的组合问题。

#include#include#include#include#includeusing namespace std;__int64 p;__int64 jieji[200005];__int64 ppow(__int64 n,__int64 m){ __int64 s=1,k=n%p; while(m>0) { if(m&1)

void init(){ s[0]=0; s[1]=0; s[2]=1; int i; for (i=3;i<=100;i++) s[i]=(i-1)*(s[i-1]+s[i-2])%mod;}

下面是BC官方题解：这题的本质是求 x1+x2+x3+……+xm=n 且 x1>x2,x3,x4,……,xm 这个方程有多少整数解。可以枚举 x1.设当前 x1=u,那么问题转化为求 x1+x2+x3+……+xm=n−u 且 u>x2,x3,x4,……,xm，可以用生成函数来解决。每一个项的生成函数是f(x)=1+x+x2+……+xu−1=1−xu1−x那么m−1个相加的生成函数是G