poj 2409 Polya

最新推荐文章于 2021-02-18 10:52:23 发布

caoyu01

最新推荐文章于 2021-02-18 10:52:23 发布

阅读量422

点赞数

分类专栏： Polya 文章标签：数论

本文链接：https://blog.csdn.net/u013113958/article/details/38678697

版权

Polya 专栏收录该内容

0 篇文章 0 订阅

订阅专栏

题意:

m中颜色的珠子，个数不限，求组成长度为n的项链，求不重复的个数？（重复：翻转或者旋转能重合在一起）。

解法：

这个还真是没学的很明白，网上的资料也少，硬看书，还是很费劲。

Polya.

//稍后做笔记，先等着多校被虐吧。~~~~~

/*************************************************************************
    > File Name: poj2409.cpp
    > Author: cy
    > Mail: 1002@qq.com 
    > Created Time: 2014/8/19 10:00:14
 ************************************************************************/

#include<iostream>
#include<cstring>
#include <algorithm>
#include<cstdlib>
#include<vector>
#include<cmath>
#include<stdlib.h>
#include<iomanip>
#include<list>
#include<deque>
#include<map>
#include <stdio.h>
#include <queue>

#define maxn 50000+5

#define inf 0x3f3f3f3f
  #define INF 0x3FFFFFFFFFFFFFFFLL
#define rep(i,n) for(i=0;i<n;i++)
 #define reP(i,n) for(i=1;i<=n;i++)

#define ull unsigned long long
 #define ll long long

#define cle(a) memset(a,0,sizeof(a))

using namespace std;
int gcd(int x,int y)
{
	if(y==0)return x;
	else return gcd(y,x%y);
}
int main()
{
#ifndef ONLINE_JUDGE
     freopen("in.txt","r",stdin);
     //freopen("out.txt","w",stdout);
#endif
	 int n,m;
	while(cin>>m>>n)
	{
		if(n==0&&m==0)break;
		int sum=0;
		int i,j,k;
		for(i=1;i<=n;i++)
		{
			int temp=gcd(n,i);
            sum+=(int)pow(m*1.0,gcd(n,i)*1.0);
		}
		if(n%2==1)sum+=(int)(n*(int)pow(m*1.0,(n+1)/2.0));
		else{
			sum+=(int)(n/2*pow(m*1.0,(n+2)/2.0));
			sum+=(int)(n/2*pow(m*1.0,(n)/2.0));
		}
		sum/=(n*2);
		cout<<sum<<endl;
	}
    return 0;
}