【分析】
首先,第一步tarjan缩点还是很容易想到的。
一个点至少要经过一次,所以拆点,连两条边,第一条流量为1,费用为1,第二条流量为INF,费用为0。
连通块之间连边,流量为INF,费用为0。 若块x可作为起点,则由S–>x’连边,流量为INF,费用为0。若可作为终点,则由x”–>T连边,流量为INF,费用为0。 增广次数即为答案。
(又是一道把我心态做崩的题。。一开始打成了每条边走一次╮(╯▽╰)╭。。三个小时才发现错了。。蠢得可以。)
【代码】
#include <cstdio>
#include <iostream>
#include <queue>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#define N 2005
#define M 100005
#define INF 1000000000
using namespace std;
typedef long long ll;
typedef pair<ll,ll> pa;
int read()
{
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
int n,m,cnt=1,S,T,ans,CaseT,Cnt,num_s,num_t,ind,scc,Ans,sum;
int b[M],p[N],nextedge[M],w[M],c[M];
int B[20005],P[N],Next[20005];
int Dis[N],Pre[N],Dfn[N],Low[N];
int aa[N],bb[N];
bool Flag[N],Instack[N],Mark[N];
int Belong[N];
void Add(int x,int y,int z,int cost)
{
cnt++;
b[cnt]=y;
nextedge[cnt]=p[x];
p[x]=cnt;
w[cnt]=z;
c[cnt]=cost;
}
void Anode(int x,int y,int z,int cost){
Add(x,y,z,cost);Add(y,x,0,-cost);
}
void Ins(int x,int y)
{
Cnt++;
B[Cnt]=y;
Next[Cnt]=P[x];
P[x]=Cnt;
}
void Input_Init()
{
n=read();m=read();num_s=read(),num_t=read();
Cnt=ind=scc=ans=Ans=sum=0;
for(int i=1;i<=n;i++) P[i]=0;
for(int i=1;i<=num_s;i++) aa[i]=read();
for(int i=1;i<=num_t;i++) bb[i]=read();
for(int i=1;i<=m;i++)
{
static int x,y,z;
x=read(),y=read();
Ins(x,y);
}
}
stack<int>st;
void Tarjan(int x)
{
Dfn[x]=Low[x]=++ind;
st.push(x);
Instack[x]=1;
for(int i=P[x];i;i=Next[i])
{
int v=B[i];
if(!Dfn[v])
{
Tarjan(v);
Low[x]=min(Low[x],Low[v]);
}
else if(Instack[v]) Low[x]=min(Low[x],Dfn[v]);
}
if(Low[x]==Dfn[x])
{
int now=-1;scc++;
while(now!=x)
{
now=st.top();st.pop();
Instack[now]=0;
Belong[now]=scc;
}
}
}
void Get_Scc()
{
for(int i=1;i<=n;i++) Dfn[i]=0;
for(int i=1;i<=n;i++) if(!Dfn[i])
Tarjan(i);
}
void Build_Graph()
{
T=scc<<1|1;
cnt=1;
for(int i=0;i<=T;i++) p[i]=Pre[i]=0;
for(int i=1;i<=n;i++) Mark[i]=0;
for(int i=1;i<=num_s;i++) Anode(0,Belong[aa[i]],INF,0);
for(int i=1;i<=num_t;i++) Anode(Belong[bb[i]]+scc,T,INF,0);
for(int i=1;i<=n;i++)
for(int j=P[i];j;j=Next[j])
if(Belong[i]!=Belong[B[j]])
Anode(Belong[i]+scc,Belong[B[j]],INF,0);
for(int i=1;i<=scc;i++)
{
Anode(i,i+scc,INF,0);
Anode(i,i+scc,1,-1);
}
}
bool Bfs()
{
queue<int>q;
q.push(S);
for(int i=0;i<=T;i++) Dis[i]=INF;
Dis[S]=0;
while(!q.empty())
{
int k=q.front();q.pop();Flag[k]=0;
for(int i=p[k];i;i=nextedge[i])
{
int v=b[i],f=w[i];
if(Dis[v]>Dis[k]+c[i]&&f)
{
Dis[v]=Dis[k]+c[i];
Pre[v]=i;
if(!Flag[v])
{
Flag[v]=1;
q.push(v);
}
}
}
}
return Dis[T]!=INF;
}
bool Mcf()
{
int Maxf=INF;
for(int i=Pre[T];i;i=Pre[b[i^1]])
Maxf=min(Maxf,w[i]);
for(int i=Pre[T];i;i=Pre[b[i^1]])
{
w[i]-=Maxf,w[i^1]+=Maxf;
ans-=c[i]*Maxf;
}
}
void MCF()
{
while(Bfs())
{
int pre=ans;
Mcf();
if(pre==ans) break;
else Ans++;
}
if(ans!=scc) printf("no solution\n");
else printf("%d\n",Ans);
}
int main()
{
CaseT=read();
while(CaseT--)
{
Input_Init();
Get_Scc();
Build_Graph();
MCF();
}
return 0;
}