数学:polya
文章平均质量分 79
skajre
这个作者很懒,什么都没留下…
展开
-
POJ 1286 Necklace of beads
注意一点,n有可能为0,所以要特殊处理,否则穿进polya函数中会导致除数为0,造成Runtime error。 // // main.cpp // Richard // // Created by 邵金杰 on 16/8/3. // Copyright © 2016年 邵金杰. All rights reserved. // #include typedef long long L原创 2016-08-03 13:55:17 · 333 阅读 · 0 评论 -
HDU 3923 Invoker
普通的ploy,套模版就行,只是有一个注意点,就是因为取模了,所以最后不能直接除2*n。 // // main.cpp // Richard // // Created by 邵金杰 on 16/8/3. // Copyright © 2016年 邵金杰. All rights reserved. // #include using namespace std; typedef lo原创 2016-08-03 16:01:50 · 430 阅读 · 0 评论 -
暴力polya模版
LL polya(LL m,LL n) { LL result=0; for(int i=1;i<=n;i++) result+=Pow(m,gcd(n,i)); if(n%2) result+=Pow(m,(n+1)/2)*n; else result+=(Pow(m,n/2)+Pow(m,n/2+1))*n/2;原创 2016-08-03 16:09:41 · 204 阅读 · 0 评论 -
POJ 2409 Let it Bead
// // main.cpp // Richard // // Created by 邵金杰 on 16/8/3. // Copyright © 2016年 邵金杰. All rights reserved. // #include #include using namespace std; typedef long long LL; LL gcd(LL a,LL b) { r原创 2016-08-03 14:04:08 · 374 阅读 · 0 评论 -
POJ 2154 Color
#include using namespace std; int n,p; int euler(int x) { int ret=1; for(int i=2;i*i<=x;i++) { if(x%i==0) { x/=i;ret*=i-1; while(x%i==0) {x/=i;ret*=原创 2016-08-03 19:26:17 · 200 阅读 · 0 评论 -
POJ 2888 Magic Bracelet
// // main.cpp // Richard // // Created by 邵金杰 on 16/8/3. // Copyright © 2016年 邵金杰. All rights reserved. // #include #include using namespace std; const int maxn=50000; int primes[maxn+1],is_pr原创 2016-08-04 17:25:10 · 186 阅读 · 0 评论