Uva10601 立方体 置换

这题的置换集合有很多，得慢慢找。。。

代码：

#include<iostream>
#include<cstdio>
#include<vector>
#include<string>
#include<queue>
#include<cstring>
#define maxn 15
#define INF 0xfffffff
#define mem(a,b) memset(a,b,sizeof(a))
#define FOR(i,s,t) for(int i=s;i<=t;i++)
#define ull unsigned long long
#define ll long long
using namespace std;
inline int rd(int &x)
{
    char c=getchar();
    while(!isdigit(c))c=getchar();
    x=0;
    while(isdigit(c))
    {
        x=x*10+c-'0';
        c=getchar();
    }
    return x;
}
int buf[10];
inline void ot(int i)
{
    int p=0;
    if(i==0) p++;
    else while(i)
        {
            buf[p++]=i%10;
            i/=10;
        }
    for(int j=p-1; j>=0; j--) putchar('0'+buf[j]);
}
int a[maxn],b[maxn];
ll c[20][20];
ll cal(int k)//循环节长度为k，把12条边分成12/k组
{
    int n=0;//总共n组，可能减掉一部分边
    for(int i=0;i<6;i++)
    {
        if(b[i]%k==0)
        {
            b[i]=b[i]/k;//b[i]能染几组
            n+=b[i];
        }
        else
        return 0;//如果有一种颜色不是k的倍数，则不能染整数组
    }
    ll ret=1;
//    cout<<"n:"<<n<<endl;
    for(int i=0;i<6;i++)
    {
        ret*=c[n][b[i]];
        n-=b[i];
    }
    return ret;
}
ll fun1()//不动
{
    memcpy(b,a,sizeof(a));
    return cal(1);//循环节长度为1
}
ll fun2()//对面旋转
{
    ll ret=0;
    memcpy(b,a,sizeof(a));
    ret+=3*2*cal(4);//90 270
    memcpy(b,a,sizeof(a));
    ret+=3*cal(2);//180
    return ret;
}
ll fun3()//对点旋转
{
    memcpy(b,a,sizeof(a));
    ll ret=0;
    ret+=4*2*cal(3);
    return ret;
}
ll fun4()//对棱旋转
{
    ll ret=0;
    memcpy(b,a,sizeof(a));
    //有两条棱是不变的,随意连，可以相同也可以不相同
    for(int i=0;i<6;i++)
    {
        for(int j=0;j<6;j++)
        {
            memcpy(b,a,sizeof(a));
            b[i]--,b[j]--;//先把不是不变的两个染色
            if(b[i]>=0&&b[j]>=0)
            {
                ret+=6*cal(2);//算一次都会改变b[i]
            }
        }
    }
    return ret;
}
ll solve()
{
    ll ret=0;
    //不动
//    cout<<ret<<endl;
    ret+=fun1();
//    cout<<ret<<endl;
    //对面旋转
    ret+=fun2();
//    cout<<ret<<endl;
    //对点旋转
    ret+=fun3();
//    cout<<ret<<endl;
    //对棱旋转
    ret+=fun4();
//    cout<<ret<<endl;
    return ret/24;//置换集合为24种
}
void init()
{
    for(int i=0;i<13;i++)
    {
        c[i][0]=1;
        c[i][i]=1;
        for(int j=1;j<i;j++)
        {
            c[i][j]=c[i-1][j-1]+c[i-1][j];
        }
    }
}
void output()
{
    for(int i=0;i<12;i++)
    {
        for(int j=0;j<=i;j++)
        {
            cout<<c[i][j]<<' ';
        }
        cout<<endl;
    }
}
int main()
{
    int t,num;
    init();
//    output();
    scanf("%d",&t);
    while(t--)
    {
        mem(a,0);
        for(int i=0;i<12;i++)
        {
            scanf("%d",&num);
            a[num-1]++;
        }
        printf("%lld\n",solve());
    }
	return 0;
}