hdu 3547 (polya定理 + 小高精）

DIY Cube
Time Limit: 2000/2000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)
Total Submission(s): 207    Accepted Submission(s): 111

Problem Description

Mr. D is interesting in combinatorial enumeration. Now he want to find out the number of ways on painting the vertexes of a cube. Suppose there are C different colors and two paintings are considered the same if they can transform from one to another by rotation.

Input
There are multiple test cases in the input, the first line of input contains an integer denoting the number of test cases.
For each test case, there are only one integer C, denoting the number of colors. (1 <= C <= 1000000000)
 
Output
For each test case, output the the number of painting ways. And if the number is equal or larger than 1015, output the last 15 digits.
 
Sample Input

3  1  2  112

Sample Output

Case 1: 1  Case 2: 23  Case 3: 031651434916928

Author
HyperHexagon
 
Source
HyperHexagon's Summer Gift (Original tasks)
 
Recommend
zhengfeng

---

polya定理的应用，需要加高精最后输出。如果想不到，你可以用置换群乘法让计算机代替你去算每种操作的循环节，我这里是已经在纸上算好的了。

假设有x种颜色。

对于一个cube有四种大置换：

1：固定对立的面旋转，共有3对对立面：

　　　　可得旋转90°与旋转270°互为逆操作，都有两个循环节，共有 3*2*x^2个不动点;

　　　　旋转180°有四个循环节，共有3*1*x^4个不动点；

　　共有3*2+3*1=9种置换。

2：固定对立的边旋转，共有6对对立边：

　　　　只可旋转180°，有四个循环节，共有6*1*x^4个不动点；

　　共有6种置换。

3：固定对立的角旋转，共有4对对立的角：

　　　　旋转120°与旋转270°互为逆操作，都有四个循环节，共有4*2*x^4 个不动点。

　　共4*2=8种置换。

4：不动：

　　　　有8个循环节，有X^8个不动点。

　　共1种置换。

所以共有24种置换。

共有x^8+17*x^4+6*x^2个不动点。

由L=1/|G| *Σ（D（ai））得：

等价类L=1/24*（x^8+17*x^4+6*x^2）；

由于保留后15位，数字也过大，用大数与小数的高精度去处理。

下面给出代码：

 #include<cstdio>
 #include<iostream>
 #include<cstring>
 #define clr(x) memset(x,0,sizeof(x))
 #define cop(x,y) memcpy(x,y,sizeof(y))
 #define LL long long
 #define lim 10
 using namespace std;
 LL ans[];
 LL xpow[];//x的幂
 LL dv[];//x^(i/2)的系数
 int dd[]={,,,,};
 void add(LL *a,LL *b);//高精加
 void mul(LL *a,int b);//高精乘
 void div(LL *a,int b);//高精除
 int main()
 {
     int T,x;
     scanf("%d",&T);
     for(int t=;t<=T;t++)
     {
         scanf("%d",&x);
         printf("Case %d: ",t);
         clr(ans);
         clr(xpow);
         xpow[]=;
         for(int i=;i<=;i++)
         {
             mul(xpow,x);
             if(i%==)
             {
                 cop(dv,xpow);
                 mul(dv,dd[i/]);
                 add(ans,dv);
             }
         }
         div(ans,);
         int v=;
         while(ans[v]== && v>=)
             v--;
         if(v>) v=;
         for(int i=v;i>=;i--)
             printf("%d",ans[i]);
         printf("\n");
     }
     return ;
 }
 void mul(LL *a,int b)
 {
     LL ret=;
     int v=;
     while(a[v]== && v>=)
         v--;
     for(int i=;i<=v || ret!=;i++)
     {
         ret=ret/lim+a[i]*(LL)b;
         a[i]=ret%lim;
     }
     return ;
 }
 void add(LL *a,LL *b)
 {
     LL ret=;
     int v=;
     while(a[v]== && b[v]== && v>=)
         v--;
     for(int i=;i<=v || ret!=;i++)
     {
         ret=ret/lim+a[i]+b[i];
         a[i]=ret%lim;
     }
     return ;

 }
 void div(LL *a,int b)
 {
     LL ret=;
     int v=;
     while(a[v]== && v>=)
         v--;
     for(int i=v;i>=;i--)
     {
         ret=ret*lim+a[i];
         a[i]=ret/(LL)b;
         ret%=(LL)b;
     }
     return ;
 }