【BZOJ】4205: 卡牌配对-最大流&线性筛

传送门:bzoj4205

---

题解

又读错题了，是至多一项互质，而非存在一项互质。。。

考虑枚举一定不互质的剩下两项的约数——只需要枚举质数，而$<200$的质数只有$num<50$个，考虑建$3nu{m}^2$个中间点分别表示$(A,B,C)$中$AB,BC,AC$不互质时的公约数情况。

$S$向X组每个点连边，X组每个点的$AB,BC,AC$组队分别向包含其质因数的中间点连边。中间点向Y组中组队包含其质因数的点连边，Y组点向T连边。容量均为1。
因为$200<2\times 3\times 5\times 7$，所以每个点最多连出去$3\times 3\times 3$条边。

跑最大流即可。

因为是分层图，所以跑dinic，匈牙利会TLE。

---

代码

一堆变量名打错，跪

#include<bits/stdc++.h>
using namespace std;
const int N=6e4+10,M=202,inf=0x7f7f7f7f;

int n,m,num,p[M];
struct P{int v[3];}a[N],b[N];
bool pri[M];

const int id[3][3]=
{
{0,1,0},
{0,2,1},
{1,2,2}
};

void pre()
{
	int i,j;
	for(i=2;i<=200;++i){
		if(!pri[i]) p[++num]=i;
		for(j=1;j<=num && i*p[j]<=200;++j){
			pri[i*p[j]]=true;
			if(i%p[j]==0) break;
		}
	}
}

namespace MF{
#define MX 4000005
int head[N<<1],to[MX],nxt[MX],w[MX],tot=1;
int S,T,dep[N<<1];

inline void lk(int u,int v,int vv)
{
    to[++tot]=v;nxt[tot]=head[u];head[u]=tot;w[tot]=vv;
    to[++tot]=u;nxt[tot]=head[v];head[v]=tot;w[tot]=0;
}

queue<int>que;
inline bool bfs()
{
	int i,j,x;
	memset(dep,0xff,(T+1)<<2);
	dep[S]=1;que.push(S);
	for(;!que.empty();){
		x=que.front();que.pop();
		for(i=head[x];i;i=nxt[i]){
			j=to[i];if(dep[j]!=-1 || (!w[i])) continue;
			dep[j]=dep[x]+1;que.push(j);
		}
	}
	return (dep[T]!=-1);
}

int dfs(int x,int f)
{
	if(x==T) return f;
	int i,j,res,ss=0;
	for(i=head[x];i;i=nxt[i]){
		j=to[i];if(dep[j]!=dep[x]+1 || (!w[i])) continue;
		res=dfs(j,min(f-ss,w[i]));
		if(!res) continue;
		w[i]-=res;w[i^1]+=res;ss+=res;
		if(ss==f) return ss;
	}
	if(!ss) dep[x]=-1;
	return ss;
}

int sol()
{
	int i,j,k,x,y,bs=3*num*num,re=0;S=0;T=bs+n+m+1;
	for(x=1;x<=n;++x) lk(S,bs+x,1);
	for(x=1;x<=m;++x) lk(bs+n+x,T,1);
	for(x=1;x<=n;++x)
	 for(y=0;y<3;++y)
	  for(i=1;i<=num;++i) if(a[x].v[id[y][0]]%p[i]==0)
	    for(j=1;j<=num;++j) if(a[x].v[id[y][1]]%p[j]==0)
	     lk(bs+x,id[y][2]*num*num+(i-1)*num+j,1);
	for(x=1;x<=m;++x)
	 for(y=0;y<3;++y)
	  for(i=1;i<=num;++i) if(b[x].v[id[y][0]]%p[i]==0)
	    for(j=1;j<=num;++j) if(b[x].v[id[y][1]]%p[j]==0)
	     lk(id[y][2]*num*num+(i-1)*num+j,bs+n+x,1);
	for(;bfs();) re+=dfs(S,inf);
	return re;
}
#undef MX
}

int main(){
	int i,j;pre();
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;++i) 
	 for(j=0;j<3;++j)
	  scanf("%d",&a[i].v[j]);
	for(i=1;i<=m;++i) 
	 for(j=0;j<3;++j)
	  scanf("%d",&b[i].v[j]);
	printf("%d\n",MF::sol());
	return 0;
}