CF453B Little Pony and Harmony Chest（数学+状压dp）

给定一个序列a，每个数不超过30，让你找出一个合法的序列b使得
$\sum\limits_{i=1}^n|a_i-b_i|$最小。
合法的序列b定义为b中的数两两互质。

数列中的数两两互质这个要求怎么转化呢？其实就是要求每个质因数只能出现在一个数中。我们首先观察到b中的数最大为2*max{a}-1.（否则变成1肯定更优）

所以b每个位置可填的数也就只有60种。可能的质因数只有17个。于是我们可以状压dp来做。

dp[i][S]->前i位,已经存在S状态的质因数了的最小答案
frac[i]->i的质因数，枚举b来转移即可。
记录路径即可。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 110
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,a[N],prime[N],tot=0,bin[20],frac[N],pre[N][132000];//frac[i]->i的质因数
bool notprime[N];int dp[N][132000];//dp[i][S]->前i位,已经存在S状态的质因数了的最小答案
inline void getprime(int nn){
    notprime[1]=1;
    for(int i=2;i<=nn;++i){
        if(!notprime[i]) prime[++tot]=i;
        for(int j=1;prime[j]*i<=nn;++j){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}
inline void print(int i,int S){
    if(!i) return;
    print(i-1,S^frac[pre[i][S]]);printf("%d ",pre[i][S]);
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();for(int i=1;i<=n;++i) a[i]=read();
    getprime(60);bin[0]=1;
    for(int i=1;i<=tot;++i) bin[i]=bin[i-1]<<1;
    for(int i=1;i<=59;++i)
        for(int j=1;j<=tot;++j)
            if(i%prime[j]==0) frac[i]|=bin[j-1];
    memset(dp,inf,sizeof(dp));dp[0][0]=0;
    for(int i=1;i<=n;++i)
        for(int k=1;k<2*a[i];++k){
            int S=(~frac[k])&(bin[tot]-1);
            for(int s=S;;s=S&(s-1)){
                if(dp[i-1][s]+abs(k-a[i])<dp[i][s|frac[k]])
                    dp[i][s|frac[k]]=dp[i-1][s]+abs(k-a[i]),pre[i][s|frac[k]]=k;
                if(!s) break;
            }
        }int ans=0;
    for(int i=1;i<=bin[tot]-1;++i) if(dp[n][i]<dp[n][ans]) ans=i;
    print(n,ans);
    return 0;
}