给定一个序列a,每个数不超过30,让你找出一个合法的序列b使得
∑i=1n|ai−bi|
最小。
合法的序列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;
}