poj 2409 Let it Bead Polya计数

旋转可以分为n种置换,对应的循环个数分别是gcd(n,i),个i=0时不动,有n个

翻转分为奇偶讨论,奇数时有n种置换,每种有n/2+1个

偶数时有n种置换,一半是n/2+1个,一半是n/2个

啃论文,PPT,各种书好久才看懂Polya定理,最近做数学题做的严重怀疑自己的智商。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include<algorithm>
#include<map>
#include<cstring>
using namespace std;
typedef long long ll;
ll gcd(ll a,ll b) {return a%b==0?b:gcd(b,a%b);}
ll quickpow(ll m,ll n)
{
    ll ans=1;
    while(n)
    {
        if(n&1) ans=ans*m;
        n=(n>>1);
        m=m*m;
    }
    return ans;
}
int main()
{
    ll c,n;
    while(~scanf("%lld%lld",&c,&n))
    {
        if(c+n==0) break;
        ll ans=quickpow(c,n);
        for(int i=1;i<n;i++) ans+=quickpow(c,gcd(n,i));
        if(n&1) ans+=n*quickpow(c,n/2+1);
        else ans+=(n/2*quickpow(c,n/2+1)+n/2*(quickpow(c,n/2)));
        printf("%lld\n",ans/(2*n));
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

TommyTT

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值