【JZOJ 4675】 Double-row

Description

科学家温斯顿在一张超长的白纸上写下了两行数，每一行数有N个。
但他写完后觉得看起来有点不和谐。他希望重新编排，使得每一行数中没有相同的数。
他每次可以调换同一列的两个数。
请帮他找到操作次数最少的方案。
N<=50000
数值Xi满足1<=X<=100000

Analysis

刚开始看像路径覆盖模型，但是并不知道如何实现交换操作
事实上是转化为图论模型
**Ai向Bi连有向边
最终会形成若干联通块
联通块的形状只可能是环或链
交换同一列两个数相当于将一条边换向**
遍历一遍，可以O(n)计算答案
还有一种连边方法，将每一列看成一个点，读者可自行思考

Code

#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define efo(i,v) for(int i=last[v];i;i=next[i])
using namespace std;
typedef long long ll;
const int N=500005,T=100000;
int n,mx,m,all,a[N],b[N];
int tot,now,ans,to[N],next[N],last[N],out[N],wei[N];
bool bz[N];
void read(int &t)
{
    t=0;
    int p=1;char ch;
    for(ch=getchar();ch<'0' || ch>'9';ch=getchar())
        if(ch=='-') p=-1;
    for(;'0'<=ch && ch<='9';ch=getchar()) t=t*10+ch-'0';
    t*=p;
}
void link(int u,int v,int w)
{
    to[++tot]=v,next[tot]=last[u],wei[tot]=w,last[u]=tot;
    out[u]++;
}
void dfs(int v)
{
    bz[v]=1;
    efo(i,v)
        if(!bz[to[i]])
        {
            m+=wei[i],all++;
            dfs(to[i]);
        }
}
int main()
{
    read(n);
    fo(i,1,n) read(a[i]),mx=max(mx,a[i]);
    fo(i,1,n) read(b[i]),mx=max(mx,b[i]);
    tot=1;
    fo(i,1,n) link(a[i],b[i],0),link(b[i],a[i],1);
    fo(i,1,mx)
        if(out[i]==1 && !bz[i])
        {
            m=all=0;
            dfs(i);
            ans+=min(m,all-m);
        }
    fo(i,1,mx)
        if(!bz[i])
        {
            m=all=0;
            dfs(i);
            all++,m+=wei[next[last[i]]]^1;
            ans+=min(m,all-m);
        }
    printf("%d",ans);
    return 0;
}