POJ 2409 ploya定理

//题意：给定颜色种数和环上的珠子总数，问有多少种染色方案（通过旋转和翻转相同的算同一种）。
//
//分析：polya定理。在这里只谈一下polya定理是如何应用的。对于排成一排的带编号的小球，按照某一种方案改变其中一些球的放置顺
//序，可以称之为置换。每一种置换方法可以用两排数字來表示，第一排数字和第二排数字一一对应，第一排数字表示小球的原来位置（1～n），
//第二排数字表示小球交换后的位置。现在我们有n个小球，m种颜色。有k种置换方法，我们认为能通过置换方法交换位置后变成同一种染色情况
//（颜色的排列状况相同，忽略小球编号），则我们认为这些互相通过置换能达到的状态为同一种染色方法。我们现在要求总共有多少种染色方法。
//要计算方法数，我们先要计算k种置换方法中每种置换方法中含有的环数，即建立一个图，有n个点，把每个置换方法两排数字中的上下一一对应的数字
//对看成边的起点和终点，计算这个图中有几个环。我们设环数分别为c1~ck。那么染色方法数为(m^c1+m^c2+...+m^ck)/k。以上就是polya定理，这里要
//注意的是置换方法集合必须是群，需要满足封闭性，即如果把通过该集合中的若干个方法连续进行置换压缩成一个置换方法（用两排数子表示），那么
//这种新的置换方法也必须属于该集合。
//
//这道题可以这样以来就是一道赤裸裸的polya定理题了。
//
//旋转：
//n种旋转方法每种旋转i个格（1<=i<=n）循环结有gcd（i，n）个
//翻转：
//
//（1）这种是经过某个顶点i与中心的连线为轴的翻转，由于n为偶数，有对称性，所以此种共n/2种翻转：
//
//（2）这种是以顶点i和i+1的连线的中点与中心的连线为轴的翻转，同样，根据对称性，也有n/2种翻转：
//
//所以给定长度n，共有2n种置换。
//
//题意：用k种颜色对n个珠子构成的环上色，旋转翻转后相同的只算一种，求不等价的着色方案数。
//
//Burnside定理的应用：
//
//当n为奇数时，有n种翻转，每种翻转都是以一个顶点和该顶点对边的中点对称。有k^(n/2+1)*n种。
//
//当n为偶数时，有n种翻转，其中一半是以两个对应顶点，另一半是以两条对边对称。有k^(n/2+1)*n/2+k^(n/2)*n/2种。
//
//考虑旋转：枚举旋转角度360/n*i，(0<i<=n)，也就是一个置换。经过该置换，颜色仍保持不变的着色方案有k^GCD(n,i)种。
//
//
//
//一个长度为n的环，每i个上同一种颜色，可以上多少种颜色。
//
//假设起点在x，则x，x+i，x+2*i，……，x+k*i，……
//
//假设在第t次，第一次回到起点，则x=(x+t*i)%n => t*i%n=0 => t=LCM(i,n)/i=n*i/GCD(n,i)/i=n/GCD(n,i)。
//
//那么可以上n/t种颜色，即n/(n/GCD(n,i))种，所以旋转的着色方案有k^GCD(n,i)种
//
//
//大概学学

//(LL) m
//(ll)n;
//
//k为颜色数
//不同的方法数=几个循环群数*pow（k,循环节数）之和  然后再除以所有的置换数之和
//  ans +=      n pow(k,(n+1)/2)                                      n+n  奇数
//                                                                    n+n/2+n/2 偶数
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<string>
#define ll long long
using namespace std;

ll gcd(ll a,ll b)
{
    if(b==0 ) return a;
    return gcd(b,a%b);
}

ll pow(int a,int b)
{
    return b==0?1:pow(a,b-1)*a;
}

int main()
{
    int n,k;
    while(scanf("%d%d",&k,&n)!=EOF)
    {
        if(n==0 && k==0)
            break;
        ll ans=0;
        if(n%2)
            ans+=n*pow(k,((n/2)+1));
        else
            ans+=n/2*pow(k,((n/2)+1))+n/2*pow(k,(n/2));


        for(int i=0;i<n;i++)
            ans+=pow(k,gcd(i,n));
        ans=ans/2/n;
    cout<<ans<<endl;
    }
}