题意
传送门 POJ 1286
题解
组合问题中,要求把旋转和翻转之后相同的状态看成是本质相同的状态,考虑 P o ˊ l y a Pólya Poˊlya 计数定理。设颜色有 m m m 种。对于旋转,有旋转 [ 0 , n ) [0,n) [0,n) 个 2 π / n 2\pi /n 2π/n 的角度的置换,即 n n n 个置换群。考虑旋转 k k k 个位置后和原来相同的染色方案个数,将颜色与第 i i i 个石头相同的石头的集合叫做 i i i 的轨迹,个数为 n / g c d ( k , n ) n/gcd(k,n) n/gcd(k,n)。因此,旋转 k k k 个位置后和原来相同的染色方案数就是 m g c d ( k , n ) m^{gcd(k,n)} mgcd(k,n)。
翻转同理,若 n n n 为奇数,对称轴必经过某一个点,共 n n n 种置换;若 n n n 为偶数,对称轴经过两点或两点连线的中点,各 n / 2 n/2 n/2 种置换。同时考虑旋转与翻转,共 2 n 2n 2n 个置换群。
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
int gcd(int a, int b)
{
return b == 0 ? a : gcd(b, a % b);
}
ll _pow(ll x, ll n)
{
ll res = 1;
while (n)
{
if (n & 1)
res *= x;
x = x * x;
n >>= 1;
}
return res;
}
int main()
{
int n, m = 3;
while (~scanf("%d", &n) && n != -1)
{
if (n == 0)
{
puts("0");
continue;
}
ll res = 0;
for (int i = 0; i < n; ++i)
res += _pow(m, gcd(i, n));
if (n & 1)
res += n * _pow(m, (n + 1) / 2);
else
res += n / 2 * (_pow(m, (n + 2) / 2) + _pow(m, n / 2));
res /= 2 * n;
printf("%lld\n", res);
}
return 0;
}