POJ 1486 - Sorting Slides 二分匹配 必须边

1. 枚举每一个边，看看最大匹配数是否会减少即可
2. 
   
3. 
   
4. #include<cstdio>  
5. #include<cstring>  
6. using namespace std;  
7. const int maxn=26+5;  
8.   
9. struct Max_Match  
10. {  
11.     int n,m;  
12.     bool g[maxn][maxn];  
13.     bool vis[maxn];  
14.     int left[maxn];  
15.   
16.     void init(int n)  
17.     {  
18.         this->n=n;  
19.         memset(g,0,sizeof(g));  
20.         memset(left,-1,sizeof(left));  
21.     }  
22.   
23.     bool match(int u)  
24.     {  
25.         for(int v=1;v<=n;v++)if(g[u][v] && !vis[v])  
26.         {  
27.             vis[v]=true;  
28.             if(left[v]==-1 || match(left[v]))  
29.             {  
30.                 left[v]=u;  
31.                 return true;  
32.             }  
33.         }  
34.         return false;  
35.     }  
36.   
37.     int solve()  
38.     {  
39.         int ans=0;  
40.         for(int i=1;i<=n;i++)  
41.         {  
42.             memset(vis,0,sizeof(vis));  
43.             if(match(i)) ans++;  
44.         }  
45.         return ans;  
46.     }  
47. }MM;  
48.   
49. int xmin[maxn],ymin[maxn],xmax[maxn],ymax[maxn];  
50. struct  
51. {  
52.     int x;  // edge[i].x=x 表示第i个矩形配对的数字 是x;  
53.     bool ok;//标记该边是否是 必需边  
54. }edge[maxn];  
55.   
56. int main()  
57. {  
58.     int n,kase=0;  
59.     while(scanf("%d",&n)==1&&n)  
60.     {  
61.         MM.init(n);  
62.         memset(edge,0,sizeof(edge));  
63.         for(int i=1;i<=n;i++)  
64.         {  
65.             scanf("%d%d%d%d",&xmin[i],&xmax[i],&ymin[i],&ymax[i]);  
66.         }  
67.         for(int i=1;i<=n;i++)  
68.         {  
69.             int x,y;  
70.             scanf("%d%d",&x,&y);  
71.             for(int j=1;j<=n;j++)  
72.             {  
73.                 if(xmin[j]<=x&&x<=xmax[j]&&ymin[j]<=y&&y<=ymax[j])  
74.                     MM.g[j][i]=true;  
75.             }  
76.         }  
77.         MM.solve();  
78.         int edge_num=n;  
79.         for(int i=1;i<=n;i++)  
80.         {  
81.             edge[MM.left[i]].x=i;  
82.             edge[MM.left[i]].ok=true;  
83.         }  
84.         for(int i=1;i<=n;i++)//尝试删除第i条匹配边  
85.         {  
86.             int j=edge[i].x;  
87.             MM.g[i][j]=false;//删除此边  
88.             memset(MM.left,-1,sizeof(MM.left));  
89.             int num = MM.solve();  
90.             if(num == n)//删除边后,匹配数不变  
91.             {  
92.                 edge[i].ok=false;  
93.                 edge_num--;  
94.             }  
95.             MM.g[i][j]=true;//还原此边  
96.         }  
97.         printf("Heap %d\n",++kase);  
98.         if(edge_num==0) printf("none\n");  
99.         else  
100.         {  
101.             for(int i=1;i<=n;i++)if(edge[i].ok)  
102.                 printf("(%c,%d) ",i-1+'A',edge[i].x);  
103.             printf("\n");  
104.         }  
105.         printf("\n");  
106.     }  
107.     return 0;  
108. }