HDU6363

先要知道几个结论gcd(2^a-1,2^b-1)=2^(gcd(a,b))-1;还得知道gcd(F[cnt1],F[cnt2])=F[gcd(cnt1,cnt2)]，这两个结论都可用数学归纳法来推导，推导第二个还要知道F数组相邻两个元素互质，那么此题对于一个gcd来说，他的贡献是2^(F(gcd))-1这么多。而对于把n分成k份使得他的gcd为gcd的倍数，n本书相同，k个书架不同，那么方案数为c(n/gcd+k-1,k-1)，使用隔板法的思想，而要求得gcd 就为当前枚举的gcd，则使用莫比乌斯函数进行容斥，是gcd的倍数-是gcd两倍的方案数，-是gcd3倍的方案数+是gcd4倍的方案数。ans为枚举所有的gcd然后莫比乌斯反演得到的每个gcd的方案数和，最后总方案书是c(n+k-1,k-1)，乘以他的逆元即可。总复杂度为T*n*lnn，但实际上只跑了143ms，使用int运算还是十分快的。

#include
#include
#define maxl 2000010
#define mod 1000000007

int n,k;
int mu[maxl],p[maxl],F[maxl];
int jc[maxl],jcinv[maxl];
bool no[maxl];
int ans;

inline void shai()
{
	int t,j;
	mu[1]=1;no[1]=true;
	for(int i=2;i>=1;
	}
	return ans;
}

inline int c(int n,int r)
{
	if(n<0 || r<0 || r>n) return 0;
	if(r==0 || n==r) return 1;
	return (1ll*((1ll*jc[n]*jcinv[n-r])%mod)*jcinv[r]%mod)%mod;
}

inline void init()
{
	shai();
	jc[0]=1;
	for(int i=1;i=1;i--)
		jcinv[i]=(1ll*jcinv[i+1]*(i+1))%mod;
	F[0]=0;F[1]=1;
	for(int i=2;i