POJ 1486（二分图匹配）二分图的完全匹配的必须边

题意：给你n个幻灯片，每个幻灯片有个数字编号1~n,现在给每个幻灯片用A~Z进行编号，那么问有多少个幻灯片的数字和字母是唯一对应的

二分图的完全匹配的必须边，裸的。

这道题目可以用

http://www.cnblogs.com/proverbs/archive/2012/08/29/2662638.html

方法来解。。但是我的第一反应就是删边。于是乎，用删边的方法做的

 

先做一遍匹配，匈牙利。然后对于匹配中的每条边都删除一遍，如果不能达到完全匹配，则这条边是必须边。

大致思想就是这样了

 

代码里几个我犯错的地方写了注释

 

View Code 
 1 #include <cstdio>
 2 #include <cstdlib>
 3 #include <cstring>
 4 #define N 40
 5 #define M 1600
 6 using namespace std;
 7 int n,x1[N],x2[N],y1[N],y2[N],x,y,head[N],next[M],to[M],bh[N][N],ban,cnt,linky[N],first_ans,cas;
 8 bool vis[N];
 9 void add(int u,int v)
10 {
11     to[cnt]=v; next[cnt]=head[u]; head[u]=cnt++;
12 }
13 void read()
14 {
15     memset(head,-1,sizeof head);cnt=0;
16     for(int i=1;i<=n;i++) scanf("%d%d%d%d",&x1[i],&x2[i],&y1[i],&y2[i]);
17     for(int i=1;i<=n;i++)
18     {
19         scanf("%d%d",&x,&y);
20         for(int j=1;j<=n;j++)
21             if(x>x1[j]&&x<x2[j]&&y>y1[j]&&y<y2[j])
22             {
23                 bh[i][j]=cnt;//两点连线的编号 
24                 add(i,j);
25             }
26     }
27 }
28 bool dfs(int u)
29 {
30     for(int i=head[u];~i;i=next[i])
31         if(!vis[to[i]]&&ban!=i)
32         {
33             vis[to[i]]=true;
34             if(linky[to[i]]==-1||dfs(linky[to[i]]))
35             {
36                 linky[to[i]]=u;
37                 return true;
38             }
39         }
40     return false;
41 }
42 void first_match()
43 {
44     memset(linky,-1,sizeof linky);
45     for(int i=1;i<=n;i++)
46     {
47         memset(vis,0,sizeof vis);
48         if(dfs(i)) first_ans++;
49     }
50 }
51 void go()
52 {
53     ban=-1;
54     first_match();
55     bool bj=false;
56     printf("Heap %d\n",cas);
57     for(int i=1,tmp;i<=n;i++)
58     {
59         tmp=linky[i];
60         ban=bh[tmp][i];//ban是删的边，注意tmp是左边的点，i是右边的点 
61         linky[i]=-1;
62         memset(vis,0,sizeof vis);
63         if(!dfs(tmp))
64         {
65             bj=true,printf("(%c,%d) ",'A'+i-1,tmp);
66             linky[i]=tmp;//必须放在括号内，因为如果dfs返回true的话说明有新的匹配，linky不能更新 
67         }
68     }
69     if(bj) printf("\n\n");
70     else printf("none\n\n");
71 }
72 int main()
73 {
74     while(scanf("%d",&n),n)
75     {
76         cas++;
77         read();
78         go();
79     }
80     return 0;
81 }

 

话说，我这个蒟蒻POJ终于到了300题了。。

还是太弱，XLk神牛每天40题呢。。

OTZZZZZZZZZZZZZZZZZZZZZZZZZZZZZ

posted @ 2012-08-29 22:05  proverbs  阅读(552)  评论(0)  编辑  收藏  举报