poj2154 Color ——Polya定理

最新推荐文章于 2021-02-17 13:26:15 发布

weixin_30563319

最新推荐文章于 2021-02-17 13:26:15 发布

阅读量116

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/Zinn/p/9251212.html

本文介绍Polya定理的应用案例，并提供了一个具体的代码实现。通过解决一个关于组合数学的问题，展示了如何使用Polya定理进行计算。文章包含了完整的C++代码，可供读者学习参考。

题目：http://poj.org/problem?id=2154

今天学了个高端的东西，Polya定理...

此题就是模板，然而还是写了好久好久...

具体看这个博客吧：https://blog.youkuaiyun.com/wsniyufang/article/details/6671122

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int const maxn=1e6+5;
int X,n,p,pri[maxn],cnt,ans;
bool vis[maxn];
int phi(int n)
{
    int ret=n;
    for(int i=2;i*i<=n;i++)
    {
        if(n%i)continue;
//        (ret-=ret/i)%=p;// 有除法，不可随时取模 
        ret-=ret/i;
        while(n%i==0)n/=i;
    }
//    if(n!=1)(ret-=ret/n)%=p;// 有除法，不可随时取模 
    if(n!=1)ret-=ret/n;
    return ret%p;
} 
int pw(int a,int b)
{
    int ret=1; a%=p;//!!!
    for(;b;b>>=1,(a*=a)%=p)
        if(b&1)(ret*=a)%=p;
    return ret;
}
int main()
{
    scanf("%d",&X);
    while(X--)
    {
        scanf("%d%d",&n,&p);
        ans=0;
        for(int i=1;i*i<=n;i++)
        {
            if(n%i)continue;
            (ans+=phi(n/i)*pw(n,i-1))%=p;
            if(i*i!=n)(ans+=phi(i)*pw(n,n/i-1))%=p;
        }
        printf("%d\n",ans);
    }
}

转载于:https://www.cnblogs.com/Zinn/p/9251212.html