bzoj4514: [Sdoi2016]数字配对（费用流）

题目传送门
。

解法：
哇跟网上做法完全不一样诶。
我本来是这样想的。
分成两排点，都是1~n。
能连的话第一排连第二排。
这样的话出现了一个问题。如果第二个数字只能用一次。
第一个数字去跟第二个数字匹配。第二个数字去跟第三个数字匹配。
这种情况我无法限制。

然后。我就想到一个比较讨巧的方法。
还是两排。
能连的话第一排连第二排。然后反过来的也要连。
这样的话就不会出现上一种情况了。
如果第一个数字去跟第二个数字匹配。表示这个费用是比较大的。
那么第二个数字去跟第一个数字匹配其实是一样的。所以会再加多一次答案。
最后答案除以2就行。。

怎么判断两个数是否满足条件匹配。
先把每个数质因数次数的和算出来。
如果能整除然后次数和相差1那么肯定相差一个质数没乘。然后就做出来了。。

翻了下跟网上的做法完全不一样嘿嘿嘿很开心啊。。

代码实现：

#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
struct node {int x,y,c;ll d,next,other;}a[110000];int len,last[1100];
void ins(int x,int y,int c,ll d) {
    int k1,k2;
    len++;k1=len;a[len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;
    a[len].next=last[x];last[x]=len;
    len++;k2=len;a[len].x=y;a[len].y=x;a[len].c=0;a[len].d=-d;
    a[len].next=last[y];last[y]=len;
    a[k1].other=k2;a[k2].other=k1;
}
int list[1100],head,tail,fbian[1100],st,ed;bool v[1100];ll d[1100];
bool spfa() {
    head=1;tail=2;list[1]=st;memset(v,false,sizeof(v));v[st]=true;
    for(int i=1;i<=ed;i++)d[i]=-99999999999;d[st]=0;
    while(head!=tail) {
        int x=list[head];
        for(int k=last[x];k;k=a[k].next) {
            int y=a[k].y;
            if(a[k].c>0&&d[y]<d[x]+a[k].d) {
                d[y]=d[x]+a[k].d;fbian[y]=k;
                if(v[y]==false) {v[y]=true;list[tail++]=y;if(tail==ed+1)tail=1;}
            }
        }head++;v[x]=false;if(head==ed+1)head=1;
    }return d[ed]==-99999999999?false:true;
}
int ans=0;ll sum=0;bool bk=false;
void wh() {
    int x=ed,flow=999999999;
    while(x!=st) {int k=fbian[x];flow=min(flow,a[k].c);x=a[k].x;}
    if(sum+flow*d[ed]>=0) {sum+=flow*d[ed];ans+=flow;}
    else {ans+=sum/(-d[ed]);bk=true;}
    x=ed;while(x!=st) {int k=fbian[x];a[k].c-=flow;a[a[k].other].c+=flow;x=a[k].x;}
    return;
}
int n;
int A[210],B[210];ll C[210];
int f[210];
int main() {
    scanf("%d",&n);st=2*n+1;ed=st+1;
    for(int i=1;i<=n;i++) {
        scanf("%d",&A[i]);
        int t=int(sqrt(double(A[i]+1)));int x=A[i];
        for(int j=2;j<=t;j++) {
            while(x%j==0) {x/=j;f[i]++;}
        }if(x!=1)f[i]++;
    }len=0;memset(last,0,sizeof(last));
    for(int i=1;i<=n;i++) {scanf("%d",&B[i]);ins(st,i,B[i],0);ins(i+n,ed,B[i],0);}
    for(int i=1;i<=n;i++) scanf("%lld",&C[i]);
    for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)if(i!=j) {
        if(A[i]%A[j]==0&&f[i]-f[j]==1) {
            ins(i,j+n,999999999,C[i]*C[j]);
            ins(j,i+n,999999999,C[i]*C[j]);
        }
    }bk=false;
    while(spfa()==true) {wh();if(bk==true)break;}
    printf("%d\n",ans/2);
    return 0;
}