链接:https://www.nowcoder.net/acm/contest/75/C
来源:牛客网
最近古yu迷上了《恋与制作人》,天天嚷嚷着白起我老公(死gay死gay的),白起又向古yu提出了一个问题:
给你n种卡牌(数量无限),将其摆成长度为m的圆环的方案数。
(如果这个环可以通过若干次旋转或者反方向读((abc)(cba)是一样的),则认为他们是一样的)
给你n种卡牌(数量无限),将其摆成长度为m的圆环的方案数。
这时古yu一脸茫然,大哭特苦不能守护白起了。古yu希望你能帮助他解决这个问题去守护白起。
这里的n就相当于polya定理的m种颜色。。。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<iomanip>
#include<vector>
#include<map>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const ll mod=1000000007;
ll quickmod(ll a,ll b,ll c)
{
a=a%c;
ll res=1;
while(b)
{
if(b&1)
res=(res*a)%c;
a=a*a%c;
b>>=1;
}
return res;
}
ll gcd(ll a,ll b)
{
if(b==0)
return a;
return gcd(b,a%b);
}
int main()
{
ll t,n,m,i;
ll res;
scanf("%lld",&t);
while(t--)
{
scanf("%lld%lld",&n,&m);
res=0;
for(i=1;i<=m;i++)
{
res=(res+quickmod(n,gcd(m,i),mod))%mod;
}
if(m%2==0)
res=(res+((quickmod(n,m/2,mod)+quickmod(n,m/2+1,mod))*m/2)%mod)%mod;
else
res=(res+(quickmod(n,m/2+1,mod)*m)%mod)%mod;
//res=res/(2*m);
printf("%lld\n",res*quickmod(2*m,mod-2,mod)%mod);//除法求余转化为求逆元
}
}
代码中注释掉的res=res/(2*m);是这样想的,resm每一步res肯定比mod小了,在除以2*m肯定比mod小,那就不用再求余了,
这是错的!!!本来3/5和1/5不一样