【bzoj1119】[POI2009]SLO

Description

对于一个1-N的排列(ai)，每次你可以交换两个数ax与ay(x<>y)，代价为W(ax)+W(ay) 若干次交换的代价为每次交换的代价之和。请问将(ai)变为(bi)所需的最小代价是多少。

Input

第一行N。第二行N个数表示wi。第三行N个数表示ai。第四行N个数表示bi。 2<=n<=1000000 100<=wi<=6500 1<=ai,bi<=n ai各不相等，bi各不相等 (ai)<>(bi) 样例中依次交换数字(2,5)(3,4)(1,5)

Output

一个数，最小代价。

Sample Input

6

2400 2000 1200 2400 1600 4000

1 4 5 3 6 2

5 3 2 4 6 1
Sample Output

11200
HINT

感谢MT大牛贡献译文.

题解

转自hzwer

1.找出初始状态和目标状态。明显，目标状态就是排序后的状态。
2.画出置换群，在里面找循环。例如，数字是8 4 5 3 2 7
明显，目标状态是2 3 4 5 7 8，能写为两个循环：
(8 2 7)(4 3 5)。
3.观察其中一个循环，明显地，要使交换代价最小，应该用循环里面最小的数字2，去与另外的两个数字，7与8交换。这样交换的代价是：
sum – min + (len – 1) * min
化简后为：
sum + (len – 2) * min
其中，sum为这个循环所有数字的和，len为长度，min为这个环里面最小的数字。

4.考虑到另外一种情况，我们可以从别的循环里面调一个数字，进入这个循环之中，使交换代价更小。例如初始状态：
1 8 9 7 6
可分解为两个循环：
(1)(8 6 9 7)，明显，第二个循环为(8 6 9 7)，最小的数字为6。我们可以抽调整个数列最小的数字1进入这个循环。使第二个循环变为：(8 1 9 7)。让这个1完成任务后，再和6交换，让6重新回到循环之后。这样做的代价明显是：
sum + min + (len + 1) * smallest
其中，sum为这个循环所有数字的和，len为长度，min为这个环里面最小的数字，smallest是整个数列最小的数字。

5.因此，对一个循环的排序，其代价是sum – min + (len – 1) * min和sum + min + (len + 1) * smallest之中小的那个数字。

6.我们在计算循环的时候，不需要记录这个循环的所有元素，只需要记录这个循环的最小的数及其和。

代码

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<set>
#include<complex>
#define ll long long
#define mod 10000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int N=1000005;
int n,w[N],a[N],b[N],pos[N],flag[N],MIN=1e9;
ll ans;
inline void solve(int x)
{
    flag[x]=1;int t=x,len=1,mn=w[x];ll sum=w[x];
    while (a[pos[t]]!=x)
    {
        t=a[pos[t]];
        flag[t]=1;
        sum+=w[t];
        mn=min(w[t],mn);
        len++;
    }
    ans+=min(sum+(ll)(len-2)*mn,sum+mn+(ll)(len+1)*MIN);
}
int main()
{
    n=read();
    for (int i=1;i<=n;i++)
    {
        w[i]=read();
        MIN=min(w[i],MIN);
    }
    for (int i=1;i<=n;i++) a[i]=read();
    for (int i=1;i<=n;i++) b[i]=read(),pos[b[i]]=i;
    for (int i=1;i<=n;i++)if (!flag[i]) solve(i);
    cout<<ans;
    return 0;
}