POJ - 1087 (A Plug for UNIX)

题意：一个房间里有 n 个插座，现在有 m 个插头送过来匹配，还有 k 个转换器，一种插头匹配一种插座，转换器的格式 A B 表示把 A类型的插头变成B类型的插头，问 m 个插头最少剩几个不能匹配成功；

 

分析：网络流建图，我们用map给不同类型的插头插座编号，并设最大编号为res，设一种插头插座的编号为 u，则 2*u-1 代表 插头，2*u 代表插座，

① 0 为 源点，由 0 向 m个插头连一条容量为1的边；

② 2*res+1为汇点，由 n 个插座向汇点连接一条容量为1的边；

③ 对于每种插头插座编号 u ，由 2*u-1 向 2*u 连一条容量正无穷的边；

④对于 K 个转换器，由 A类型的插座向B类型的插头连一条容量为1的边；

建图完成就可以直接套最大流模板了；

 

 

代码：

1) EK

#include<map> 
#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

const int N = 1000;
const int M = 205*205;
const int INF = 1E9+7;

int n;

struct Edge{int to;int next;int value;}e[M<<1];
struct Pre{int node;int id;}pre[M<<1];//pre[i].node表示编号为i的点最短路的上一个点,pre[i].id表示最短路上连接i点的边的编号
int head[N],cnt=-1;//编号从0开始,原因见下
bool vis[N];
queue<int> q;
void add(int from,int to,int value)
{
    cnt++;
    e[cnt].to=to;
    e[cnt].value=value;
    e[cnt].next=head[from];
    head[from]=cnt;
}
bool bfs(int s,int t)//用来寻找s,t的最短路并记录,如果s,t不连通则返回0
{
    q=queue<int>();//清空队列
    memset(vis,0,sizeof(vis));
    memset(pre,-1,sizeof(pre));
    pre[s].node=s;
    vis[s]=1;
    q.push(s);
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        for(int i=head[x];i>-1;i=e[i].next)
        {
            int now=e[i].to;
            if(!vis[now]&&e[i].value)//忽略流量为0的边
            {
                pre[now].node=x;//用pre记录最短路
                pre[now].id=i;
                vis[now]=1;
                if(now==t)return 1;//找到
                q.push(now);
            }
        }
    }
    return 0;
}

int EK(int s,int t)
{
    int ans=0;
    while(bfs(s,t))
    {
        int minv=INF;
        for(int i=t;i!=s;i=pre[i].node)
            minv=min(minv,e[pre[i].id].value);
        for(int i=t;i!=s;i=pre[i].node)
        {
            e[pre[i].id].value-=minv;
            e[pre[i].id^1].value+=minv;//x^1表示x边的反向边,此方法仅在边的编号从0开始时有效
        }
        ans+=minv;
    }
    return ans;
}

string A[105],B[105],K[105][2];
int a,b,k,res;
map<string,int>mp;
int sum[10000000];

void input(){
	res=0;
	cin>>a;
    for(int i=1;i<=a;i++){
    	cin>>A[i];
    	if(mp[A[i]]==0) mp[A[i]]=++res;
		sum[mp[A[i]]*2]+=1; 
	}
    cin>>b;
    for(int i=1;i<=b;i++){
    	cin>>B[i]>>B[i];
    	if(mp[B[i]]==0) mp[B[i]]=++res;
    	sum[mp[B[i]]*2-1]+=1; 
	}
    cin>>k;
    for(int i=1;i<=k;i++){
    	cin>>K[i][0]>>K[i][1];
    	if(mp[K[i][0]]==0) mp[K[i][0]]=++res;
    	if(mp[K[i][1]]==0) mp[K[i][1]]=++res;
    	sum[(mp[K[i][0]]*2)*10000+(mp[K[i][1]])*2-1]+=1; 
	}
    
}

void work(){
	for(int i=1;i<=a;i++){
    	int num=mp[A[i]]*2;
    	if(sum[num]>0){
    		add(num,n,sum[num]);
    		add(n,num,0);
    	    //cout<<num<<' '<<n<<' '<<sum[num]<<endl;
    	    sum[num]=0;
		}
	}
	for(int i=1;i<=b;i++){
    	int num=mp[B[i]]*2-1;
    	if(sum[num]>0){
    		add(0,num,sum[num]);
    		add(num,0,0);
    	    //cout<<0<<' '<<num<<' '<<sum[num]<<endl;
    	    sum[num]=0;
		}
	}
    for(int i=1;i<=k;i++){
    	int Y=mp[K[i][1]],X=mp[K[i][0]];
    	int num=(X*2)*10000+(Y*2-1);
    	if(sum[num]>0){
    		add(X*2,Y*2-1,num);
    		add(Y*2-1,X*2,0);
    	    //cout<<X*2<<' '<<Y*2-1<<' '<<sum[num]<<endl;
    	    sum[num]=0;
		}
	}
    for(int i=1;i<=res;i++){
    	add(i*2-1,i*2,INF);
    	add(i*2,i*2-1,0);
	}
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	memset(head,-1,sizeof(head));
    input();
    n=2*res+1;
    work();
    cout<<b-EK(0,n);
}

 

2) Dinic

#include<map> 
#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

const int N = 1000;
const int M = 205*205;
const int INF = 1E9+7;

int n;

int dis[N];
queue<int> q;

struct Edge{int to;int value;int next;}e[M<<1];
int head[N],cnt=-1;
void add(int from,int to,int value)
{
    cnt++;
    e[cnt].to=to;
    e[cnt].value=value;
    e[cnt].next=head[from];
    head[from]=cnt;
}

bool bfs(int s,int t)//bfs功能和EK算法的相似,不同的是Dinic中的bfs要求出所有点到源点s的最短路dis[i]
{
    q=queue<int>();//清空队列
    memset(dis,-1,sizeof(dis));
    dis[s]=0;
    q.push(s);
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        for(int i=head[x];i>-1;i=e[i].next)
        {
            int now=e[i].to;
            if(dis[now]==-1&&e[i].value!=0)
            {
                dis[now]=dis[x]+1;
                q.push(now);
            }
        }
    }
    return dis[t]!=-1;
}
int cur[N];

int dfs(int x,int t,int maxflow)
{
    if(x==t)return maxflow;
    int ans=0;
    for(int i=cur[x];i>-1;i=e[i].next)
    {
        int now=e[i].to;
        if(dis[now]!=dis[x]+1||e[i].value==0||ans>=maxflow)continue;
        cur[x]=i;//此路可行，记录此路
        int f=dfs(now,t,min(e[i].value,maxflow-ans));
        e[i].value-=f;
        e[i^1].value+=f;
        ans+=f;
    }
    return ans;
}

int Dinic(int s,int t)
{
    int ans=0;
    while(bfs(s,t))
    {
        for(int i=0;i<=n;i++) cur[i]=head[i];
		//memcpy(cur,head,sizeof(head)); //初始化，用这个会T
        ans+=dfs(s,t,INF);
    }
    return ans;
}



string A[105],B[105],K[105][2];
int a,b,k,res;
map<string,int>mp;
int sum[10000000];

void input(){
	res=0;
	cin>>a;
    for(int i=1;i<=a;i++){
    	cin>>A[i];
    	if(mp[A[i]]==0) mp[A[i]]=++res;
		sum[mp[A[i]]*2]+=1; 
	}
    cin>>b;
    for(int i=1;i<=b;i++){
    	cin>>B[i]>>B[i];
    	if(mp[B[i]]==0) mp[B[i]]=++res;
    	sum[mp[B[i]]*2-1]+=1; 
	}
    cin>>k;
    for(int i=1;i<=k;i++){
    	cin>>K[i][0]>>K[i][1];
    	if(mp[K[i][0]]==0) mp[K[i][0]]=++res;
    	if(mp[K[i][1]]==0) mp[K[i][1]]=++res;
    	sum[(mp[K[i][0]]*2)*10000+(mp[K[i][1]])*2-1]+=1; 
	}
    
}

void work(){
	for(int i=1;i<=a;i++){
    	int num=mp[A[i]]*2;
    	if(sum[num]>0){
    		add(num,n,sum[num]);
    		add(n,num,0);
    	    //cout<<num<<' '<<n<<' '<<sum[num]<<endl;
    	    sum[num]=0;
		}
	}
	for(int i=1;i<=b;i++){
    	int num=mp[B[i]]*2-1;
    	if(sum[num]>0){
    		add(0,num,sum[num]);
    		add(num,0,0);
    	    //cout<<0<<' '<<num<<' '<<sum[num]<<endl;
    	    sum[num]=0;
		}
	}
    for(int i=1;i<=k;i++){
    	int Y=mp[K[i][1]],X=mp[K[i][0]];
    	int num=(X*2)*10000+(Y*2-1);
    	if(sum[num]>0){
    		add(X*2,Y*2-1,num);
    		add(Y*2-1,X*2,0);
    	    //cout<<X*2<<' '<<Y*2-1<<' '<<sum[num]<<endl;
    	    sum[num]=0;
		}
	}
    for(int i=1;i<=res;i++){
    	add(i*2-1,i*2,INF);
    	add(i*2,i*2-1,0);
	}
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	memset(head,-1,sizeof(head));
    input();
    n=2*res+1;
    work();
    cout<<b-Dinic(0,n);
}