POJ 1286 2409 Polya计数

NeighThorn

于 2016-09-11 09:51:44 发布

阅读量346

点赞数 1

分类专栏： poj题解

本文链接：https://blog.csdn.net/neighthorn/article/details/52502454

版权

poj题解专栏收录该内容

94 篇文章 0 订阅

订阅专栏

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#define LL long long
using namespace std;
int c,s;
inline LL gcd(LL a,LL b){
    return b==0?a:gcd(b,a%b);
}
inline LL power(int a,int b){
    LL ans=1;
    while(b){
        if(b&1)
            ans=ans*a;
        a=a*a,b>>=1;
    }
    return ans;
}
inline LL polya(void){
    LL sum=0;
    for(int i=1;i<=s;i++)
        sum+=power(c,gcd(s,i));
    if(s&1)
        sum+=s*power(c,(s>>1)+1);
    else
        for(int i=1;i<=s;i++){
            if(i&1)
                sum+=power(c,(s>>1)+1);
            else
                sum+=power(c,s>>1);
        }
    return sum/2/s;
}
signed main(void){
    c=3;
    while(scanf("%d",&s)&&s!=-1)
        printf("%lld\n",s==0?0:polya());
    return 0;   
}