/**
1.一个n环的全排列是 (n-1)! : (第一个位置是环头尾相连的地方.不用改变)
2.插入一个n环中 有n个位置可以插入.
3.匈牙利算法 :它是一种用增广路径求二分图最大匹配的算法。
http://blog.csdn.net/dark_scope/article/details/8880547
题意:有2*n个珠子,n个阴珠子,n个阳珠子.给你m个关系,i,j; 表示第i个阳珠子和第j个阴珠子如果相邻阳珠子就会变暗
且要求阳珠子和阴珠子交替出现.
思路:因为是环.我们先全排列出所有银珠子构成的环,在其中n个地方插入n个阳珠子.
对这n个 位置 进行关系建图
进行匈牙利算法 ,求出最大匹配sum.则最小变暗的阳珠子就是 n - sum个;
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 20;
///yin:枚举阴珠子的排列 dis:表示阳珠子和珠子的变暗关系,can:表示第i个位置可以接受的阳珠子(有关系)
int yin[maxn],dis[maxn][maxn],can[maxn][maxn];
int vis[maxn],yang[maxn]; ///这两个匈牙利算法的标记数组.
inline void read(int &a)
{
a = 0;
char ch = getchar();
while(ch==' '||ch == '\n') ch = getchar();
while(ch>='0'&&ch<='9'){ a=a*10+ch-'0';ch = getchar(); }
}
///匈牙利算法的递归
bool dfs(int k,int n)
{
for(int i=1;i<=n;i++)
{
if(!vis[i]&&can[k][i])
{
vis[i] = 1;
if(yang[i]==-1||dfs(yang[i],n))
{
yang[i] = k;
return true;
}
}
}
return false;
}
int solve(int n)
{
///这两个初始化位置放错了 wa 了 两次.
memset(yang,-1,sizeof(yang));
memset(can,0,sizeof(can));
for(int i=1;i<=n;i++)
{
int pre = yin[i],last = i==1? yin[n] : yin[i-1];
for(int j = 1;j<=n;j++)
{
if(!dis[j][pre] && !dis[j][last]){
can[i][j] = 1;
}
}
}
int ans = 0;
for(int i=1;i<=n;i++){
memset(vis,0,sizeof(vis));
if(dfs(i,n)) ans++;
}
return n-ans;
}
int main()
{
int n,m;
while(~scanf("%d%d",&n,&m))
{
int a,b;
memset(dis,0,sizeof(dis));
memset(yin,0,sizeof(yin));
for(int i=1;i<=n;i++) yin[i] = i;
for(int i=1;i<=m;i++){
read(a);read(b);
dis[a][b] = 1;
}
if(!n) {printf("0\n");continue;}
int ans = maxn;
do
{
ans=min(solve(n),ans);
}while(next_permutation(yin+2,yin+1+n));
printf("%d\n",ans);
}
return 0;
}
HDU-5727Necklace-项链 全排列环+匈牙利
最新推荐文章于 2024-09-03 22:33:53 发布