HDU4623 CRIME 【状压DP】【同类项合并】

题目大意:

求相邻元素互质的排列个数。

题目分析:

由于互质只与质因数有关，所以我们对于质因数种类相同的数合并为一类，特殊的，1,17,19,23是一类，因为没有数与他们不互质。

那么我们做各个位进制不同的状压DP。转移就是在末尾添加哪个数。

代码:

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<algorithm>
 #include<vector>
 #include<iostream>
 using namespace std;

 #define P __gcd

 int n,mod;
 int tot,res,maxx;

 int GCD[][];
 int a[];
 //1,2,3,5,2*3,7,2*5,11,13,2*7,3*5,3*7,2*11,2*13
 int bel[]={,,,,,,,,,,,,,,,,,,,,,,,,,,,,};
 int data[]={,,,,,,,,,,,,,};
 int sum[];
 int fac[]={,,,,,};

 int f[][];

 void GetSubSet(){
     memset(a,,sizeof(a));maxx=;
     for(int i=;i<=n;i++) a[bel[i]]++,maxx=max(maxx,bel[i]);
     sum[] = ;
     for(int i=;i<=maxx;i++) sum[i] = sum[i-]*a[i-]+sum[i-];
 }

 vector<int> v[];
 void dfs(int now,int dd,int num){
     if(now > maxx){
     v[num].push_back(dd);
     }else{
     for(int i=;i<=a[now];i++){
         dfs(now+,dd+i*sum[now],num+i);
     }
     }
 }

 int nw[];
 void work(){
     for(int i=;i<=maxx;i++) f[sum[i]][i] = ;
     for(int i=;i<=n;i++) v[i].clear();
     dfs(,,);
     for(int i=;i<n;i++){
     for(int j=;j<v[i].size();j++){
         int p=v[i][j];
         for(int k=maxx;k>=;k--){nw[k] = p/sum[k];p%=sum[k];}
         p = v[i][j];
         for(int k=;k<=maxx;k++){
         if(nw[k] ==  || f[p][k]==) continue;
         for(int l=;l<=maxx;l++){
             if(a[l]-nw[l] == ||GCD[l][k]!=) continue;
             f[p+sum[l]][l] += f[p][k];
             f[p+sum[l]][l] %= mod;
         }
         }
     }
     }
     int ans = ;
     for(int i=;i<=maxx;i++) ans = (ans+f[v[n][]][i])%mod;
     for(int i=;i<=maxx;i++)
         ans=(ans*fac[a[i]])%mod;
     printf("%d\n",ans);
 }

 int main(){
     int t; scanf("%d",&t);
     for(int i=;i<;i++)for(int j=;j<;j++)GCD[i][j]=P(data[i],data[j]);
     while(t--){
     scanf("%d%d",&n,&mod);
     memset(f,,sizeof(f));
     GetSubSet();
     work();
     }
     return ;
 }