POJ 1486 Sorting Slides

http://poj.org/problem?id=1486

题意：给n个矩形的4个边界的坐标（左上和右下），分别是：xmin、xmax、ymin、ymax。又给出四个数字的坐标。每个数字只能属于一个矩形。求出每个数字的从属关系。

题解：二分图最大匹配问题：数字和矩形的匹配。要求出每一条必须边。先求出最大匹配ans，然后删除每一条边，再进行匹配，看最大匹配是否等于ans：如果相等，则这条边不是必须边；如果  小于ans，则这条边是必须边。

注意点：删除边之后要恢复这条边；输出格式……

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <map>
#include <queue>
#include <stack>
#include <bitset>
#include <algorithm>
#include <numeric>
#include <functional>
#include <set>
#include <fstream>

using namespace std;

const int maxn=110;
int V;
int G[maxn][maxn];
int match[maxn];
int used[maxn];
int n;
struct cos{
    int xmin,xmax,ymin,ymax;
}s[maxn];
struct numb{
    int x,y;
}num[maxn];

bool dfs(int v)
{
    for (int i=0; i<n; i++) {
        if (!used[i]&&G[v][i]) {
            used[i]=1;
            if (match[i]==-1||dfs(match[i])) {
                match[i]=v;
                return true;
            }
        }
    }
    return false;
}

int bipartite_matching()
{
    int res=0;
    memset(match, -1, sizeof(match));
    for (int v=0; v<V; v++) {
       // if (match[v]<0) {
            memset(used, 0, sizeof(used));
            if (dfs(v)) {
                res++;
            }
        //}
    }
    return res;
}
int main()
{
    //freopen("/Users/apple/Desktop/POJ 1486/POJ 1486/in", "r", stdin);
    // freopen("/Users/apple/Desktop/POJ 1486/POJ 1486/out", "w", stdout);
    int mycase=0;
    while ((scanf("%d",&n))!=EOF&&n!=0) {
        V=n;
        mycase++;
        memset(G, 0, sizeof(G));
        for (int i=0; i<n; i++) {
            scanf("%d%d%d%d",&s[i].xmin,&s[i].xmax,&s[i].ymin,&s[i].ymax);
        }
        for (int i=0; i<n; i++) {
            scanf("%d%d",&num[i].x,&num[i].y);
        }
        for (int i=0; i<n; i++) {
            for (int j=0; j<n; j++) {
                if ((num[j].x>=s[i].xmin)&&(num[j].x<=s[i].xmax)&&(num[j].y>=s[i].ymin)&&(num[j].y<=s[i].ymax)) {
                    G[i][j]=1;
                }
            }
        }
        printf("Heap %d
",mycase);
        int ans=bipartite_matching();
        bool flag=false;
        for (int i=0; i<n; i++) {
            for (int j=0; j<n; j++) {
                if (G[i][j]==0) {
                    continue;
                }
                G[i][j]=0;
                if (bipartite_matching()<ans) {
                    flag=true;
                    printf("(%c,%d) ",'A'+i,j+1);
                }
                G[i][j]=1;
            }
        }
        if (flag==0) {
            printf("none

");
        }
        else printf("

");
    }
    return 0;
}