AC日记——[Sdoi2016]数字配对 bzoj 4514
思路：
很受伤现在，，测了那么多次不过的原因就是因为INF不够大；
解法有两种：
解法1：
把n个点按照质因数个数为奇或偶分为两个点集（很容易就可以想到）；
然后，按照题目连边跑最大费用流；
当累计的能量马上就要小于0时，退出循环，输出答案；
解法2：
把n个点拆成2*n个点，也是两个集合；
如果ai[i]到ai[j]可以连边，则i连j+n，同时j连i+n；
当累计的能量马上就要小于0时，退出循环，输出答案/2；
来，上代码：
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; #define ll long long #define INF 0x7ffffffffff ll n,prime[1000050],num_prime=0,ai[205],bi[205],ci[205]; ll s,t,head[205],E[100005],V[100005],W[100005],F[100005]; ll cnt=1,pre[205],que[200005],dis[205]; bool if_prime[1000050],d1iv[205],if_[205]; inline void in(ll &now) { ll if_z=1;now=0; char Cget=getchar(); while(Cget>'9'||Cget<'0') { if(Cget=='-') if_z=-1; Cget=getchar(); } while(Cget>='0'&&Cget<='9') { now=now*10+Cget-'0'; Cget=getchar(); } now*=if_z; } void euler(ll limit) { for(ll i=2;i<=limit;i++) { if(!if_prime[i]) prime[++num_prime]=i; for(ll j=1;prime[j]*i<=limit&&j<=num_prime;j++) { if_prime[i*prime[j]]=true; if(i%prime[j]==0) break; } } } inline bool cut(ll op) { ll res=0,now=0; while(op!=1&&now<10500) { now++; while(!(op%prime[now])) op/=prime[now],res++; } return res&1; } inline void edge_add(ll u,ll v,ll w,ll f) { E[++cnt]=head[u],head[u]=cnt,W[cnt]=w,F[cnt]=f,V[cnt]=v; E[++cnt]=head[v],head[v]=cnt,W[cnt]=-w,F[cnt]=0,V[cnt]=u; } inline bool spfa() { for(ll i=s;i<=t;i++) dis[i]=-INF,if_[i]=false,pre[i]=-1; ll h=0,tail=1;que[0]=s,if_[s]=true,dis[s]=0; while(h<tail) { ll now=que[h++];if_[now]=false; for(ll i=head[now];i;i=E[i]) { if(F[i]>0&&dis[V[i]]<dis[now]+W[i]) { pre[V[i]]=i; dis[V[i]]=dis[now]+W[i]; if(!if_[V[i]]) { if_[V[i]]=true; que[tail++]=V[i]; } } } } return dis[t]!=-INF; } int main() { in(n);euler(1000040),t=n+1; for(ll i=1;i<=n;i++) in(ai[i]); for(ll i=1;i<=n;i++) in(bi[i]); for(ll i=1;i<=n;i++) in(ci[i]); for(ll i=1;i<=n;i++) { if(cut(ai[i])) { d1iv[i]=true; edge_add(i,t,0,bi[i]); } else edge_add(s,i,0,bi[i]); } for(ll i=1;i<=n;i++) { if(!d1iv[i]) { for(ll j=1;j<=n;j++) { if(d1iv[j]) { ll a=max(ai[i],ai[j]),b=min(ai[i],ai[j]); if(a==b||a%b) continue; if(!if_prime[a/b]) edge_add(i,j,ci[i]*ci[j],INF); } } } } ll ans=0,cii=0; while(spfa()) { ll now=t,pos=INF; while(pre[now]!=-1) { pos=min(pos,F[pre[now]]); now=V[pre[now]^1]; } now=t; while(pre[now]!=-1) { F[pre[now]]-=pos; F[pre[now]^1]+=pos; now=V[pre[now]^1]; } if(cii+dis[t]*pos<0) { ans+=cii/(-dis[t]); break; } ans+=pos,cii+=dis[t]*pos; } cout<<ans; return 0; }
#include <cmath> #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; #define maxn 205 #define maxm 1000005 #define ll long long #define INF 0x7ffffffffff ll que[maxm],E[maxm],V[maxm],F[maxm],W[maxm],cnt=1,bi[maxn]; ll n,head[maxn<<1],s,t,dis[maxn<<1],pre[maxn<<1],ai[maxn],ci[maxn]; bool if_[maxn<<1]; inline void in(ll &now) { ll if_z=1;now=0; char Cget=getchar(); while(Cget>'9'||Cget<'0') { if(Cget=='-') if_z=-1; Cget=getchar(); } while(Cget>='0'&&Cget<='9') { now=now*10+Cget-'0'; Cget=getchar(); } now*=if_z; } inline void edge_add(ll u,ll v,ll w,ll f) { E[++cnt]=head[u],V[cnt]=v,W[cnt]=w,F[cnt]=f,head[u]=cnt; E[++cnt]=head[v],V[cnt]=u,W[cnt]=-w,F[cnt]=0,head[v]=cnt; } inline bool judge(ll op) { ll lit=sqrt(op); for(ll i=2;i<=lit;i++) if(!(op%i)) return false; return true; } inline bool spfa() { for(ll i=s;i<=t;i++) dis[i]=-INF,pre[i]=-1,if_[i]=false; ll h=0,tail=1;dis[s]=0,que[0]=s,if_[s]=true; while(h<tail) { ll now=que[h++];if_[now]=false; for(ll i=head[now];i;i=E[i]) { if(F[i]>0&&dis[V[i]]<dis[now]+W[i]) { pre[V[i]]=i,dis[V[i]]=dis[now]+W[i]; if(!if_[V[i]]) if_[V[i]]=true,que[tail++]=V[i]; } } } return dis[t]!=-INF; } int main() { freopen("menci_pair.in","r",stdin); freopen("menci_pair.out","w",stdout); in(n),t=n<<1|1; for(ll i=1;i<=n;i++) in(ai[i]); for(ll i=1;i<=n;i++) { in(bi[i]); edge_add(s,i,0,bi[i]); edge_add(i+n,t,0,bi[i]); } for(ll i=1;i<=n;i++) in(ci[i]); for(ll i=1;i<=n;i++) { for(ll j=1;j<=n;j++) { if(ai[i]<=ai[j]||ai[i]%ai[j]) continue; if(judge(ai[i]/ai[j])) { edge_add(i,j+n,ci[i]*ci[j],INF); edge_add(j,i+n,ci[i]*ci[j],INF); } } } ll ans=0,cii=0; while(spfa()) { ll now=t,pos=INF; while(pre[now]!=-1) { pos=min(pos,F[pre[now]]); now=V[pre[now]^1]; } now=t; while(pre[now]!=-1) { F[pre[now]]-=pos; F[pre[now]^1]+=pos; now=V[pre[now]^1]; } if(cii+dis[t]*pos<0) { ans+=cii/abs(dis[t]); break; } ans+=pos,cii+=dis[t]*pos; } cout<<ans/2; return 0; }
