B. Little Pony and Harmony Chest(状压dp+数论)

https://codeforces.com/problemset/problem/453/B

---

思路：

题意要求b序列中两两互质，也就是说b序列中每个数的质因子分解不能有相同。而ai<=30,若b取到59时，|1-30|=|30-59|=29，其实bi最大范围在59就可以满足要求了，而这里面所有数质因子分解最多16个。所以状压考虑这个质因子集合出现的状态。然后最后记录前驱输出方案数。

#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<cstdio>
#include<algorithm>
#define debug(a) cout<<#a<<"="<<a<<endl;
using namespace std;
const int maxn=1<<17;
typedef long long LL;
inline LL read(){LL x=0,f=1;char ch=getchar();	while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;}
LL a[110];
LL dp[105][maxn];
LL vis[105][maxn];
LL s[70];
int primes[20]={2,3 ,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59};
void dfs(LL pos,LL S){
     if(pos==0) return;
     ///LL laststate=state^s[dp[pos][state]];
     LL u=vis[pos][S];
     dfs(pos-1,S^s[u]);
     cout<<u<<" ";
}
int main(void){
   cin.tie(0);std::ios::sync_with_stdio(false);
   LL n;cin>>n;
   for(LL i=1;i<=n;i++) cin>>a[i];
   ///处理每个数的分解质因数集合
   for(LL i=1;i<=60;i++){
       for(LL j=0;j<16;j++){
           LL a=i;
           if(a%primes[j]==0){
              a/=primes[j];
              s[i]|=(1<<j);
           }
       }
   }
   memset(dp,0x3f,sizeof(dp));
   memset(vis,-1,sizeof(vis));
   vis[0][0]=0;
   for(LL j=0;j<(1<<16);j++) dp[0][j]=0;
   for(LL i=1;i<=n;i++){///前i个数
       for(LL j=0;j<(1<<16);j++){///当前拥有的质因数分解状态
           for(LL k=1;k<60;k++){
               if((s[k]&j)!=0) continue;///该数填进去和已经有的质因子集合碰撞导致存在两两不互质
               LL lastval=dp[i][j];
               dp[i][j]=min(dp[i][j],dp[i-1][j^s[k]]+abs(a[i]-k));
               if(lastval>dp[i][j]){///状态被更新了
                  vis[i][j]=k;///当前状态填k这个数字
               }
           }
       }
   }
   LL ans=0x3f3f3f3f;LL id=0;///最后的质因子集合状态
   for(LL j=0;j<(1<<16);j++){
       if(ans>dp[n][j]&&vis[n][j]!=-1){
          ans=dp[n][j];
          id=j;
       }
   }
   dfs(n,id);
   return 0;
}