題目鏈接：http://acm.nyist.net/JudgeOnline/problem.php?pid=688

此題題解 不太懂，因為對這些概念，定理太模糊，理解起來比較困難，不過想想還是應該把代碼寫出來；

題意：給你一個數 n ，代表 n 種顏色和n個珠子，問你可以組合多少種長度為n的項鏈；不需要用掉n種顏色，項鏈的旋轉和翻轉都是為同一條

題解：?http://pan.baidu.com/s/1ntwPVkd#dir/path=%2FACM%E5%A5%97%E9%A2%98%2F09%E5%B9%B4%20-%203rd%20Central%20South%20China%20Programming%20Contest%2Fcsc2009%2FC%2Fsolution

點進去之后點擊C.ppt

公式：

#include<stdio.h> #include<math.h> #include<string.h> #include<algorithm> using namespace std; typedef long long LL; const LL mod = 20090531; LL x,y; LL euler(LL p)//歐拉定理 {LL res = p;for(LL i = 2;i <= sqrt(p);i++){if(p%i==0){while(p%i==0) p/=i;res = res/i*(i-1);}}if(p > 1) res = res / p * (p-1);return res%mod; } LL fast_pow(LL a,LL b)//快速冪 {LL res = 1;while(b){if(b&1) res = (res * a) % mod;a = (a * a) % mod;b >>= 1;}return res; } void edgcd(LL a,LL b)//求逆元 {if(b == 0){x = 1;y = 0;return;}edgcd(b,a%b);LL temp = x;x = y;y = temp - a/b*y; } void solve(LL m) {LL res = 0;for(LL i = 1;i <= sqrt(m);i++){if(m%i==0){if(m/i != i) res = (res + euler(i)*fast_pow(m,m/i-1))%mod;res = (res + euler(m/i)*fast_pow(m,i-1))%mod;}}if(m&1){res = (fast_pow(m,(m+1)/2) + res) % mod;edgcd(2,mod);}else{res = (fast_pow(m,m/2+1) + fast_pow(m,m/2) + 2 * res) % mod;edgcd(4,mod);}while(x < 0) x = (x + mod) % mod;printf("%lld\n",x * res % mod); } int main() {int n;LL m;scanf("%d",&n);while(n--){scanf("%lld",&m);solve(m);} }

總結

以上是生活随笔為你收集整理的nyoj-Color the necklace（Ploya定理 + 欧拉函数 + 扩展欧几里得（求逆元））的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。