# BZOJ 4514: [Sdoi2016]数字配对

## 4514: [Sdoi2016]数字配对

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1606  Solved: 608
[Submit][Status][Discuss]

一行一个数，最多进行多少次配对

3
2 4 8
2 200 7
-1 -2 1

4

## HINT

n≤200，ai≤10^9，bi≤10^5，∣ci∣≤10^5

## 代码：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#include<map>
//by NeighThorn
#define inf 0x3f3f3f3f3f3f3f
using namespace std;

const int maxn=200+5,maxm=32000+5,maxe=100000+5;

int n,a[maxn],b[maxn],c[maxn];
int cnt,no[maxn],vis[maxm],pri[maxm];
int S,T,hd[maxn],fl[maxe],to[maxe],nxt[maxe],Min[maxn],from[maxn];
long long w[maxe],dis[maxn];

map<int,int> mp;

inline void prework(void){
for(int i=2;i<=32000;i++){
if(!vis[i])
vis[i]=1,pri[++cnt]=i,mp[i]=1;
for(int j=1;j<=cnt&&pri[j]*i<=32000;j++){
vis[i*pri[j]]=1;
if(i%pri[j]==0) break;
}
}
}

inline void add(int x,int y,int s,long long l){
w[cnt]= l;fl[cnt]=s;to[cnt]=y;nxt[cnt]=hd[x];hd[x]=cnt++;
w[cnt]=-l;fl[cnt]=0;to[cnt]=x;nxt[cnt]=hd[y];hd[y]=cnt++;
}

inline bool spfa(void){
for(int i=S;i<=T;i++) dis[i]=-inf,Min[i]=0x3f3f3f3f;
queue<int> q;q.push(S),dis[S]=0;vis[S]=1;
while(!q.empty()){
int top=q.front();q.pop();vis[top]=0;
for(int i=hd[top];i!=-1;i=nxt[i])
if(fl[i]&&dis[to[i]]<dis[top]+w[i]){
from[to[i]]=i;
dis[to[i]]=dis[top]+w[i];
Min[to[i]]=min(Min[top],fl[i]);
if(!vis[to[i]])
vis[to[i]]=1,q.push(to[i]);
}
}
return dis[T]!=-inf;
}

inline long long find(void){
for(int i=T;i!=S;i=to[from[i]^1])
fl[from[i]]-=Min[T],fl[from[i]^1]+=Min[T];
return dis[T]*Min[T];
}

inline int mcmf(void){
long long t,mincost=0,maxflow=0;
while(spfa()){
t=find();
if(mincost+t>=0) mincost+=t,maxflow+=Min[T];
else{
maxflow+=mincost/abs(dis[T]);
return maxflow;
}
}
return maxflow;
}

signed main(void){
scanf("%d",&n);prework();S=0;
memset(hd,-1,sizeof(hd));T=n+1;
memset(vis,0,sizeof(vis));
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
for(int i=1;i<=n;i++) scanf("%d",&b[i]);
for(int i=1;i<=n;i++) scanf("%d",&c[i]);
for(int i=1,tmp;i<=n;i++){
tmp=a[i];
for(int j=1;j<=cnt;j++)
while(tmp%pri[j]==0)
no[i]++,tmp/=pri[j];
if(tmp>1) no[i]++,mp[tmp]=1;
}
for(int i=1;i<=n;i++)
if(no[i]&1)
else
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(a[i]%a[j]==0&&mp.find(a[i]/a[j])!=mp.end()){