原題連結
題意簡述
在序列中選擇若干個數,使得其 \(gcd=1\) 且對應代價最小
實施
假設答案裡,\(a_i\) 是最後一個選的,代表 \(i\) 前面存在某些數的組合的 \(gcd\) 與 \(a_i\) 互質
揹包+狀壓
再遍歷前面的數 \(j\) 和狀態,代表選 \(j\) 時,數 \(i\) 的質因子集合的狀態
code
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll N=114514;
ll dp[1024];
ll yz[305][14];
ll cnt[306]={0};
void solve()
{
int n;
cin>>n;
vector<int> c(n+5);
for(ll i=1;i<=n;i++)
{
ll x;
cin>>x;
for(ll j=2;j*j<=x;j++)
{
if(x%j==0)
{
yz[i][++cnt[i]]=j;
while(x%j==0) x/=j;
}
}
if(x>1) yz[i][++cnt[i]]=x;
}
for(int i=1;i<=n;i++) cin>>c[i];
ll ans=2e9;
for(int i=1;i<=n;i++)
{
memset(dp,0x3f,sizeof dp);
int bit=(1<<cnt[i])-1;
dp[bit]=c[i];
for(int j=1;j<i;j++)
{
ll tem=0;
for(int x=1;x<=cnt[i];x++)
{
bool flag=0;
for(int y=1;y<=cnt[j];y++)
{
if(yz[i][x]==yz[j][y])
{
flag=1;
break;
}
}
if(flag) tem|=(1<<(x-1));
}
for(int k=0;k<=bit;k++) dp[k&tem]=min(dp[k&tem],dp[k]+c[j]);
}
ans=min(ans,dp[0]);
}
if(ans==2e9) cout<<"-1\n";
else cout<<ans<<'\n';
}
int main()
{
ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
int t=1;
//cin>>t;
while(t--) solve();
return 0;
}