hdu6675 度度熊与排列(二分图最小字典序匹配)

题目

度熊有一个机器，这个机器有一个1∼M 的排列 p[1..M] 当作参数，

若丢进一个长度为 M 的字符串，此机器会将此字符串重新排列后再输出，

重新排列的方式为：原本第 i 个位置的字符会变到第 p[i] 个位置。

举例来说，当 M = 3，p[1]=3,p[2]=1,p[3]=2，

那么丢 "abc" 进入这个机器后，机器会输出"bca";若丢进的是 "ded"，那么机器会输出 "edd"。

某天，度熊不小心忘记这个机器的参数了，只记得参数的长度是 M，

于是他丢了 N 个长度为M 的字符串进去，并记录下对于每个字符串机器的输出结果，

请你根据这些结果，帮度熊找回这个机器的参数。

若有多组参数都满足度熊的记录，请输出字典序最小的排列作为参数。

若并不存在任何参数满足度熊的记录，请输出 -1。

T<=100，N<=20，M<=50，字符只由小写字母构成

思路来源

https://blog.csdn.net/weixin_43851525/article/details/99709594

题解

先对n个可能的集合取交集，然后考虑二分图匹配

第一行i的所有可能位置记为pos[i]，那么i就应该向pos[i]连边，

最后，要让求得的左边的匹配集cx[i]的字典序最小

方法：增序连边，倒序匹配

对于每个i，pos[i]按字典序增序连，这样对于i，保证i搜到的第一个是能匹配的字典序最小的

对于n个位置，应从n到1倒序连，这样每次后来的若能匹配到字典序最小的，会把先来的挤到一个字典序更大的上

从而保证了，后来的字典序不会比挤掉的大，而后来的实际位于序列前面，所以保证了开头的字典序小

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=25;
const int M=55;
int t,n,m;
ll dp[N][M],ans[M],now;
char a[M],b[M];
bool ok;
vector<int>E[M]; 
int cx[M*2],cy[M*2],res;
bool vis[M*2];
struct edge{int to,next,w;}e[M*M*2];
void init()
{
    for(int i=1;i<M;++i)
    E[i].clear();
}
void add(int u,int v)
{
    E[u].push_back(v);
}
bool dfs(int u)
{
    for(int i=0;i<E[u].size();++i)
    {
        int v=E[u][i];
        if(!vis[v])
        {
            vis[v]=1;
            if(cy[v]==-1||dfs(cy[v]))
            {
                cx[u]=v;
                cy[v]=u;
                return 1;
            }
        }
    }
    return 0;
}
int hungary(int n)
{
    int res=0;
    memset(cx,-1,sizeof cx);
    memset(cy,-1,sizeof cy);
    for(int u=n;u>=1;--u)
    {
        if(cx[u]==-1)
        {
            memset(vis,0,sizeof vis);
            res+=dfs(u);
        }
    }
    return res;
}
int main()
{
    scanf("%d",&t);
    while(t--)
    {
        ok=1;
        memset(dp,0,sizeof dp);
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;++i)
        {
            scanf("%s%s",a+1,b+1);
            for(int j=1;j<=m;++j)
            {
                now=0;
                for(int k=1;k<=m;++k)
                {
                    if(a[j]==b[k])
                    dp[i][j]|=(1ll<<k);
                }
            }
        }
        for(int i=1;i<=n;++i)
        {
            for(int j=1;j<=m;++j)
            {
                if(i==1)ans[j]=dp[i][j];
                else ans[j]&=dp[i][j];
                if(!ans[j])ok=0;
            }
        }
        if(!ok){puts("-1");continue;}
        init();
        for(int i=1;i<=m;++i)
        {
            for(int j=1;j<=m;++j)
            if(ans[i]>>j&1)add(i,j);
        }
        res=hungary(m);
        if(res!=m){puts("-1");continue;}
        for(int i=1;i<=m;++i)
        printf("%d%c",cx[i],i==m?'\n':' ');
    }
    return 0;
}