（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;}