https://vjudge.net/problem/UVA-10294
题目
某人看到商店里面有项链和手镯，是由一圈珠子构成的。有k种颜色的珠子，每种珠子随便用，问n个珠子能制作成多少种不同的项链和手镯。
项链或手镯可以旋转，如果旋转后能够重合的两个项链或者手镯看作一种项链或手镯。
手镯除了可以旋转以外还可以翻转。如果翻转后能够重合的两个手镯看作一种手镯。
输入n,k，0<n<51，0<t<11
输出项链和手镯的种数，保证答案不超过11位数。
题解
用Burnside引理 https://oi-wiki.org/math/permutation-group/
$leftlvert X/G ight vert = frac{1}{leftlvert G ight vert} sum_{gin G} leftlvert X^g ight vert$
$ X^g={x|xin X,g(x)=x} $
那么对于项链，有旋转0次、旋转1次、旋转2次……旋转n-1次，如果要求旋转后不变（整个项链看成一个元素，旋转看成置换），那么旋转i次就要求第0,i,2i,3i……个珠子颜色相同，可以得到有$lcm(n,i)/i=n/gcd(n,i)$个珠子颜色相同，这就把项链分成了$gcd(n,i)$份，每一份有t种涂色方式，因此得到项链种数为
[frac{1}{n}sum_{i=0}^{n-1} {t^{gcd(n,i)}}]
剩下的还是直接套公式就可以了……
#include<cstdio>
using namespace std;
typedef long long ll;
int gcd(int a, int b) {
return b==0? a:gcd(b,a%b);
}
ll qpow(unsigned long long a, int b) {
ll ans=1;
for(;b;b>>=1) {
if(b&1) ans*=a;
a*=a;
}
return ans;
}
int main() {
int n,t;
while(~scanf("%d%d", &n, &t)) {
ll ans=0;
for(int i=0; i<n; i++) {
ans += qpow(t, gcd(n,i));
}
printf("%lld ", ans/n);
if(n&1) {
ans+=n*qpow(t,(n+1)/2);
printf("%lld
", ans/2/n);
} else {
ans+=n/2*(qpow(t,(n+2)/2)+qpow(t,n/2));
printf("%lld
", ans/2/n);
}
}
}