POJ 3270 Cow Sorting 置换群

Cow Sorting

接触到置换群的概念,对于一个(1~n)的一个排列a1, a2, a3...an

    1   2   3   4   5 ... n

   a1  a2 a3 a4 a5...an

看作是一个置换,那么将其写为若干个不相交的循环的乘积形式(A1, A2, ... Ap1)(B1, B2, ... Bp2)... ...,例如

   1 2 3 4 5 6

   3 5 6 4 2 1

可以写为(1 3 6)(2 5)(4)的形式,这里可以看成有3个循环(可以理解为上面的1 3 6这个集合,对应在下面也是 1 3 6)

 

那么对应在每一个循环中,可以通过最多Ki - 1此交换使得这个循环有序,这里的Ki表示第i个循环的集合的大小。也就是说每交换一次可以使得一个元素回到自己的位置,由于交换的代价是这两个元素值的和,那么根据贪心的思想,每次选择用这个循环中的最小值将其他元素交换到对应的位置上去。那么总代价就是

SUMi + (Ki - 2) * Mi

这里Mi表示第i个循环的最小值,SUMi表示第i个循环的和,Ki表示这个循环的元素的个数。

 

另外,从另一个角度出发,如果整个序列的最小值非常小,那么我们判断是不是可以先用整个序列的最小值把第i个循环的最小值Mi换出来,然后将这个集合(循环)交换到有序后,再将Mi交换回来,这样总代价就是

SUMi + (Ki + 1) * MIN + Mi

其中MIN表示整个序列的最小值。

 

那么,最后的答案就是  ans = SIGMA{ SUMi + min{(Ki - 2) * Mi,  (Ki + 1) * MIN + Mi} }

 

 1 #pragma comment(linker, "/STACK:1677721600")
 2 #include <map>
 3 #include <set>
 4 #include <stack>
 5 #include <queue>
 6 #include <cmath>
 7 #include <ctime>
 8 #include <bitset>
 9 #include <vector>
10 #include <cstdio>
11 #include <cctype>
12 #include <cstdarg>
13 #include <cstring>
14 #include <cstdlib>
15 #include <iostream>
16 #include <algorithm>
17 using namespace std;
18 #define INF 0x3f3f3f3f
19 #define inf (-((LL)1<<40))
20 #define root 1, 1, n
21 #define middle ((L + R) >> 1)
22 #define lson k<<1, L, (L + R)>>1
23 #define rson k<<1|1,  ((L + R)>>1) + 1, R
24 #define mem0(a) memset(a,0,sizeof(a))
25 #define mem1(a) memset(a,-1,sizeof(a))
26 #define mem(a, b) memset(a, b, sizeof(a))
27 #define FIN freopen("in.txt", "r", stdin)
28 #define FOUT freopen("out.txt", "w", stdout)
29 #define rep(i, a, b) for(int i = a; i <= b; i ++)
30 #define dec(i, a, b) for(int i = a; i >= b; i --)
31 
32 //typedef __int64 LL;
33 typedef long long LL;
34 typedef pair<int, int> Pair;
35 const int MAXN = 300000 + 10;
36 const int MAXM = 1100000;
37 const double eps = 1e-12;
38 LL MOD = 1000000007;
39 
40 int n;
41 struct Node {
42     int a, id;
43     bool operator < (const Node &A) const {
44         return a < A.a;
45     }
46 }node[MAXN];
47 bool vis[MAXN];
48 
49 int find_ans(int p, int &mi, int &cnt) {
50     int sum = 0;
51     mi = INF;
52     cnt = 0;
53     do {
54         sum += node[p].a;
55         vis[p] = 1;
56         cnt ++;
57         mi = min(mi, node[p].a);
58         p = node[p].id;
59     } while(!vis[p]);
60     return sum;
61 }
62 
63 int main()
64 {
65 #ifndef ONLINE_JUDGE
66     FIN;// FOUT;
67 #endif
68     while(~scanf("%d", &n)) {
69         int mi = INF;
70         rep (i, 1, n) {
71             scanf("%d", &node[i].a);
72             node[i].id = i;
73             mi = min(mi, node[i].a);
74         }
75         sort(node + 1, node + n + 1);
76         mem0(vis);
77         int cur_min, cur_cnt, ans = 0;
78         rep (i, 1, n) if(!vis[i]) {
79             int sum = find_ans(i, cur_min, cur_cnt); // 找到当前循环的元素个数,最小值
80             if(cur_cnt > 1) {// mi表示所有数字中的最小值
81                 ans += sum + min( (cur_cnt - 2) * cur_min, (cur_cnt + 1) * mi + cur_min );
82             }
83         }
84         cout << ans << endl;
85     }
86     return 0;
87 }

 

posted @ 2015-08-08 16:03  再见~雨泉  阅读(214)  评论(0编辑  收藏  举报