LGTB 与序列

LGTB 有一个长度为N 的序列A，现在他想构造一个新的长度为N 的序列B，使得B 中的任意两个数都
互质。
并且他要使ai与bi对应项之差最小
请输出最小值
输入
第一行包含一个数N 代表序列初始长度
接下来一行包含N 个数A1, A2, …, AN，代表序列A
对于40% 的数据，1 N 10
对于100% 的数据，1 N 100, 1 ai 30
输出
输出包含一行，代表最小值
样例
样例输入样例输出
51
6 4 2 8
3
样例说明
样例中B 数组可以为1 5 3 1 8

考虑到bi最多变成58，如果变成更大的数还不如变成1，而且58之内只有16个素数
并且可以证明对于两个数a > b，变成的数A >= B
所以最多只有16个最大的数变成素数，其他的数都会变成1

所以直接对于前16个数dp，dp[i][j]代表到第i个数，哪些素因子被用过了，花费最少为多少。枚举一个数来转移即可

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<cmath>
 6 using namespace std;
 7 int f[101][(1<<16)],a[101],n,tot,pri[18],inf,ans,zt[101];
 8 bool vis[61];
 9 int abs(int a)
10 {
11     if (a<0) return -a;
12     return a;
13 }
14 bool cmp(int a,int b)
15 {
16     return a>b;
17 }
18 int main()
19 {
20     int i,j,k,p,l,x,flag;
21     cin>>n;
22     for (i=1; i<=n; i++)
23     {
24         scanf("%d",&a[i]);
25     }
26     sort(a+1,a+n+1,cmp);
27     for (i=2; i<=60; i++)
28         if (vis[i]==0)
29         {
30             for (j=2*i; j<=60; j+=i)
31                 vis[j]=1;
32         }
33     for (i=2; i<=58; i++)
34         if (vis[i]==0)
35         {
36             tot++;
37             pri[tot]=i;
38         }
39         for (i=2;i<=60;i++)
40         {
41             zt[i]=0;
42             for (j=1;j<=tot;j++)
43             if (i%pri[j]==0) 
44             {
45                 zt[i]|=(1<<j-1);
46             }
47         }
48     memset(f,127/3,sizeof(f));
49     inf=f[0][0];
50     f[0][0]=0;
51     for (i=1; i<=min(n,16); i++)
52     {
53         for (j=0; j<=(1<<16)-1; j++)
54             if(f[i-1][j]!=inf)
55             {
56                 f[i][j]=min(f[i][j],f[i-1][j]+abs(a[i]-1));
57                 for (k=2;k<=60;k++)
58                 if ((j&zt[k])==0)
59                 {
60                     f[i][j|zt[k]]=min(f[i][j|zt[k]],f[i-1][j]+abs(a[i]-k));
61                 }
62                 
63             }
64     }
65     ans=inf;
66     for (i=0; i<=(1<<16)-1; i++)
67             ans=min(ans,f[min(n,16)][i]);
68     for (i=min(n,16)+1;i<=n;i++)
69             ans+=abs(a[i]-1);
70     cout<<ans;
71 }

 

转载于:https://www.cnblogs.com/Y-E-T-I/p/7723652.html