poj 2409 poj 1286【马克】

polya！！！~~~

好吧就是水题…………1A的说。。。！

然后，对于每次旋转，我是暴力去找的循环数……听说用gcd可以，，让我再去学习一下。然后对于翻转，分奇偶讨论。奇数的时候，翻转的幂数每次就是floor((n+1)/2)，有n次，对于偶数的翻转，对角线的时候是floor((n+1)/2)+1，有n/2次，对称轴的时候和奇数的是一样的，也有n/2次。

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
using namespace std;
int a[2][33];
int vis[33];
int main()
{
    int n,c;
    while(scanf("%d%d",&c,&n)!=EOF&&n+c)
    {
        int i,j,k;
        for(i=1;i<=n;i++) a[0][i]=i;
        int flag=1;
        double sum=0;
        sum+=pow((double)c,(double)n);
        for(i=1;i<n;i++)
        {
            int tmp=a[!flag][1];
            for(j=1;j<n;j++)
                a[flag][j]=a[!flag][j+1];
            a[flag][n]=tmp;

        //    for(j=1;j<=n;j++) cout<<a[flag][j]<<"  "; cout<<endl;
            int d=n;
            int len=0;
            memset(vis,0,sizeof(vis));
            while(d)
            {
                j=1;
                while(vis[j]) j++;
                k=j;
                while(1)
                {
                    d--;
                    vis[j]=1;
                    j=a[flag][j];
                    if(j==k) break;
                }
                len++;
            }
        //    cout<<"len="<<len<<endl;
            sum+=pow((double)c,len);
            flag=!flag;
        }
        double m=floor((n+1)/2.0);
        if(n%2)
        sum+=n*pow((double)c,m);
        else
        {
            sum+=(n/2)*pow((double)c,m)+(n/2)*pow((double)c,m+1);
        }
        cout<<sum/(2*n)<<endl;
    }
    return 0;
}

******************************************************************

那啥，，这两道一样样的嘛，就是把c改成3就好了，然后注意特判当n等于0的时候，输出0…

不对！WA了好多发QAQ！明明答案都一样了啊！不能理解！！

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
using namespace std;
int a[2][33];
int vis[33];
int main()
{
    int n,c;
    c=3;
    //freopen("in.txt","r",stdin);freopen("outt.txt","w",stdout);
    while(scanf("%d",&n)!=EOF)
    {
        if(n==-1) break;

        if(n==0) {printf("0\n");continue;}
        //int c=3;
        int i,j,k;
        for(i=1;i<=n;i++) a[0][i]=i;
        int flag=1;
        double sum=0;
        sum+=pow((double)c,(double)n);
        for(i=1;i<n;i++)
        {
            int tmp=a[!flag][1];
            for(j=1;j<n;j++)
                a[flag][j]=a[!flag][j+1];
            a[flag][n]=tmp;

        //    for(j=1;j<=n;j++) cout<<a[flag][j]<<"  "; cout<<endl;
            int d=n;
            int len=0;
            memset(vis,0,sizeof(vis));
            while(d)
            {
                j=1;
                while(vis[j]) j++;
                k=j;
                while(1)
                {
                    d--;
                    vis[j]=1;
                    j=a[flag][j];
                    if(j==k) break;
                }
                len++;
            }
        //    cout<<"len="<<len<<endl;
            sum+=pow((double)c,len);
            flag=!flag;
        }
        double m=floor((n+1)/2.0);
        if(n%2)
        sum+=n*pow((double)c,m);
        else
        {
            sum+=(n/2)*pow((double)c,m)+(n/2)*pow((double)c,m+1);
        }
        printf("%.lf\n",sum/(2*n));
        //cout<<sum/(2*n)<<endl;
    }
    return 0;
}

**********************************

看到说可以用欧拉函数优化，在mark一下。