（5943 Kingdom of Obsession）简单的二分图匹配问题

题目链接:http://acm.split.hdu.edu.cn/showproblem.php?pid=5943

题意:有n个人，每个人的标号以此是s+1到s+n，要求将所有人重新排序之后满足每个人的位置y能够保证被他的标号整除，就是数组重排之后满足每一位的a[i]%i==0。

解题思路：由于每个人都要满足a[i]%i==0,所以a[i]的值为质数时只能被1整除，所以a[i]为质数超过1个时，便返回No。所以可知匹配的数区域不大，可以直接通过匈牙利算法进行完全二分图匹配。注意:存在n和s交叉重叠的情况，重叠部分的数不用管，直接取n为min(n,s)，s为max(n,s)进行匹配。记住所有的数组都要进行合理的初始化，我在这里WA很多。

解题代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
using namespace std;
bool line[1000][1000];
int peo[1000];
int used[1000];
int prime(int n)
{
    for(int i=2;i<=sqrt(n);i++)
        if(n%i==0)
           return 0;
    return 1;
}
int find(int x,int n)
{
    int i,j;
    for(j=1;j<=n;j++)
        if(line[x][j]==true&&peo[j]==0)
    {
        peo[j]=1;
        if(used[j]==0||find(used[j],n))
        {
            used[j]=x;
            return 1;
        }
    }
    return 0;
}
int main()
{
    int T,n,s,ans,summ,t;
    scanf("%d",&T);
    for(int i=0;i<T;i++)
    {
        memset(used,0,sizeof(used));
        memset(line,false,sizeof(line));
        ans=0;
        summ=0;
        scanf("%d%d",&n,&s);
        if(n>=s)
        {
            t=s;
            s=n;
            n=t;
        }
        for(int j=1;j<=n;j++)
        {
            if(prime(s+j)==1)
                ans++;
            if(ans>=2)
            {
                printf("Case #%d: No\n",i+1);
                break;
            }
        }
        if(ans==2)
            continue;
        for(int j=1;j<=n;j++)
            for(int k=1;k<=n;k++)
                if((j+s)%k==0)
                    line[k][j]=true;
        for(int j=1;j<=n;j++)
        {
            memset(peo,0,sizeof(peo));
            summ=summ+find(j,n);
        }
        if(summ==n)
            printf("Case #%d: Yes\n",i+1);
        else
            printf("Case #%d: No\n",i+1);
    }
    return 0;
}