polya模板

最新推荐文章于 2022-02-20 10:56:43 发布

Dimension_

最新推荐文章于 2022-02-20 10:56:43 发布

阅读量153

点赞数 1

分类专栏： acm polya

本文链接：https://blog.csdn.net/weixin_43847469/article/details/86011484

版权

acm 同时被 2 个专栏收录

8 篇文章 0 订阅

订阅专栏

polya

1 篇文章 0 订阅

订阅专栏

超好用的polya模板233
百用百ac。感谢ss大佬
n个珠子 m个颜色 n很大的时候才用的到phi，一般用前面的就够了。
polya还可以用来算循环节
如
http://acmoj.shu.edu.cn/problem/448/
求循环的个数
在这里插入图片描述

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
ll fast_pow(ll a, ll b){
    ll res = 1;
    while(b>0){
        if(b&1) res = res*a;
        a = a * a;
        b >>= 1;
    }
    return res;
}
int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}
int main(){
    int n,m;
    while(~scanf("%d%d",&n,&m)){
        if(n+m==0) break;
        if(n==0||m==0){
            printf("0\n");
            continue;
        }
        ll ans = 0;
        for(int i=1;i<=n;++i){
            ans += fast_pow(m,gcd(n,i));
        }
        ans /= n;
        ans += ((fast_pow(m,(n+1)/2)+fast_pow(m,n/2+1))/2LL);
        ans /= 2;
        printf("%lld\n",ans);
    }
    return 0;
}

#include<bits/stdc++.h>
using namespace std;
const int mod = 1e9+7;
typedef long long ll;
ll fast_pow(ll a,ll b){
    ll res = 1;
    while(b>0){
        if(b&1) res = res*a%mod;
        a = a*a%mod;
        b >>= 1;
    }
    return res;
}

//get phi(n)
int phi(int n){
    int res = n;
    for(int i=2;(long long)i*i<=n;++i){
        if(n%i==0){
            res = res / i * (i-1);
            while(n%i==0){
                n /= i;
            }
        }
    }
    if(n>1) res = res/n*(n-1);
    return res;
}

int main(){
    int T,n,m,kase=0;
    long long inv2 = fast_pow(2,mod-2); //Çó2Ä£modµÄÄæÔª
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&m,&n);
        ll ans = 0;
        int tmp = n;
        for(int i=1;(long long)i*i<=n;++i){
            if(i*i==n){
                ans = (ans+phi(i)*fast_pow(m,i))%mod;
            }
            else if(n%i==0){
                ans = (ans+phi(n/i)*fast_pow(m,i)%mod+phi(i)*fast_pow(m,n/i)%mod)%mod;
            }
        }
        ans = ans*fast_pow(tmp,mod-2)%mod; //	ans / n
        ans = (ans+inv2*(fast_pow(m,(tmp+1)/2)+fast_pow(m,(tmp/2)+1)%mod)%mod)%mod;
        ans = ans*inv2%mod;
        printf("Case #%d: %lld\n",++kase,ans);
    }
    return 0;
}