来源:http://acm.hdu.edu.cn/showproblem.php?pid=3923
题意:有m种元素,要围成长度为n的圈,问有多少种方法。
思路:可以转化成有m种颜色,要围成长度为n的项链,有多少种方法。裸的Ploya定理,但是最后除的地方用到了乘法逆元。这是今天新的收获,学习了。
因为我们得到的sum值是一直不断的模Mod的结果,所以最后我们不能直接用sum/2*n % Mod来计算结果,因为此时的sum是模后的值,并不是真正的sum。这时我们需要算出2*n在Mod中的乘法逆元为x,则sum * x % Mod即为答案。
代码:
#include <iostream>
#include <cstdio>
#include <string.h>
using namespace std;
#define Mod 1000000007
__int64 gcd(__int64 a,__int64 b){
if(b == 0)
return a;
return gcd(b,a%b);
}
__int64 mi(int x,int n){
__int64 s = 1;
for(__int64 i = 1;i <= n;++i)
s = ((s%Mod) * (x%Mod))%Mod;
return s;
}
void extend_Eulid(__int64 a,__int64 b,__int64 &x,__int64 &y){
if(b == 0){
x = 1;
y = 0;
return;
}
extend_Eulid(b,a%b,x,y);
__int64 temp = x;
x = y;
y = temp - a/b * y;
}
int main(){
//freopen("1.txt","r",stdin);
//freopen("2.txt","w",stdout);
__int64 numcase;
scanf("%I64d",&numcase);
for(__int64 k = 1;k <= numcase;++k){
__int64 color,len;
scanf("%I64d%I64d",&color,&len);
__int64 sum = 0;
for(__int64 i = 1;i <= len;++i){
__int64 x = gcd(len,i);
sum += mi(color,x);
}
if(len % 2 == 0)
sum += ((len/2 * mi(color,len/2))%Mod + (len/2 * mi(color,len/2+1)) % Mod) % Mod;
else
sum += (len * mi (color,(len-1)/2+1)) % Mod;
__int64 xx,yy;
extend_Eulid(2*len,Mod,xx,yy);
xx = (xx % Mod + Mod) % Mod;
printf("Case #%I64d: %I64d\n",k,(sum%Mod * (xx%Mod)) % Mod);
}
return 0;
}
分享到:
相关推荐
HDU的1250,主要是利用高精度加法,但是代码有点繁琐,效率不是很高
HDU1059的代码
杭电ACMhdu1163
hdu1001解题报告
HDU的一题........HDU DP动态规
hdu 1574 passed sorce
hdu acm 教案 搜索入门 hdu acm 教案 搜索入门
hdu2101AC代码
搜索 dfs 解题代码 hdu1241
ACM HDU题目分类,我自己总结的大概只有十来个吧
hdu acm 教案 动态规划(1) hdu acm 教案 动态规划(1)
hdu 5007 Post Robot 字符串枚举。 暴力一下就可以了。
HDU最全ac代码
hdu-acm源代码(上百题)hdu-acm源代码、hdu-acm源代码hdu-acm源代码
hdu动态规划算法集锦
hdu 1166线段树代码
hdu题目分类
HDU图论题目分类
自己做的HDU ACM已经AC的题目