[NOIP2010] 关押罪犯

 

【问题描述】

S城现有两座监狱,一共关押着N名罪犯,编号分别为1~N。他们之间的关系自然也极不和谐。很多罪犯之间甚至积怨已久,如果客观条件具备则随时可能爆发冲 突。我们用“怨气值”(一个正整数值)来表示某两名罪犯之间的仇恨程度,怨气值越大,则这两名罪犯之间的积怨越多。如果两名怨气值为c的罪犯被关押在同一 监狱,他们俩之间会发生摩擦,并造成影响力为c的冲突事件。每年年末,警察局会将本年内监狱中的所有冲突事件按影响力从大到小排成一个列表,然后上报到S 城Z市长那里。公务繁忙的Z市长只会去看列表中的第一个事件的影响力,如果影响很坏,他就会考虑撤换警察局长。

在详细考察了N名罪犯间的矛盾关系后,警察局长觉得压力巨大。他准备将罪犯们在两座监狱内重新分配,以求产生的冲突事件影响力都较小,从而保住自己的乌纱 帽。假设只要处于同一监狱内的某两个罪犯间有仇恨,那么他们一定会在每年的某个时候发生摩擦。那么,应如何分配罪犯,才能使Z市长看到的那个冲突事件的影 响力最小?这个最小值是多少?

【输入】

输入文件名为prison.in。输入文件的每行中两个数之间用一个空格隔开。

第一行为两个正整数N和M,分别表示罪犯的数目以及存在仇恨的罪犯对数。

接下来的M行每行为三个正整数aj,bj,cj,表示aj号和bj号罪犯之间存在仇恨,其怨气值为cj。数据保证1<=aj<bj<N,0<cj<=1,000,000,000,且每对罪犯组合只出现一次。

【输出】

输出文件prison.out共1行,为Z市长看到的那个冲突事件的影响力。如果本年内监狱中未发生任何冲突事件,请输出0。


【输入样例】


4 6
1 4 2534
2 3 3512
1 2 28351
1 3 6618
2 4 1805
3 4 12884
【输出样例】


3512
【输入输出样例说明】
罪犯之间的怨气值如下面左图所示,右图所示为罪犯的分配方法,市长看到的冲突事件影响力是3512(由2 号和3 号罪犯引发)。其他任何分法都不会比这个分法更优。


【数据范围】

对于30%的数据有N≤15。

对于70%的数据有N≤2000,M≤50000。

对于100%的数据有N≤20000,M≤100000。

 【分析】

这是一道经典题目,解法有很多:并查集,二分答案+二分图染色,甚至还有人直接说可以直接跑二分图匹配……

 考虑到编写代码的复杂程度,我选用了拆点后用利用并查集维护关系然后贪心的办法。

 将每个“罪犯”拆成两个点:它本身所在的监狱和另一个监狱。每次找到原图中权值最大的一条边,尽可能地将这两个罪犯放在不同的监狱——将A,B分别与B,A的对称点合并。如果这条边对应的两个罪犯已经确定处于同一个集合,这时的边权就是答案,因为这时如果我们要将这两个罪犯放在不同的集合,就必然会与之前考虑过的关系冲突,产生更大的矛盾值,不符合题意。

 1 #include <algorithm>
 2 #include <cstdio>
 3 #include <cctype>
 4 #include <cmath>
 5 #include <vector>
 6 #include <cstring>
 7 #include <queue>
 8 
 9 #if defined DEBUG
10 FILE *in = fopen("test""r");
11 #define out stdout
12 #else
13 FILE *in = fopen("prison1.in""r");
14 FILE *out = fopen("prison1.out""w");
15 #endif
16 using namespace std;
17 inline void getint(int &x){
18     char c = fgetc(in);
19     while(!isdigit(c)) c = fgetc(in);
20     x = c - '0';
21     while(isdigit(c = fgetc(in)))x = x * 10 - '0' + c;
22 }
23 
24 /*===========================================*/
25 #define maxn (20002)
26 #define maxm (100005)
27 namespace UFS{
28     int pre[maxn<<1];
29     inline void init(int size){
30         int i;
31         for(i = 0;i <= size; ++i)
32             pre[i] = i;
33     }
34     int find(int k){
35         if(pre[k] == k)return k;
36         return pre[k] = find(pre[k]);
37     }
38 }
39 struct edge{
40     int x, y, w;
41 }E[maxm];
42 inline bool cmp(const edge&a, const edge&b){return a.w > b.w;}
43 int n, m;
44 using namespace UFS;
45 int main(){
46     int i;
47     getint(n), getint(m);
48     init(n << 1);
49     for(i = 0;i < m;++i)
50         getint(E[i].x),getint(E[i].y),getint(E[i].w);
51     sort(E, E + m, cmp);
52     for(i = 0;i < m;++i){
53         int x = find(E[i].x), y = find(E[i].y);
54         if(x == y){
55             fprintf(out"%d\n", E[i].w);
56             return 0;
57         }
58         pre[x] = find((y - 1 + n) % (n << 1) + 1);
59         pre[y] = find((x - 1 + n) % (n << 1) + 1);
60     }
61     fprintf(out"0\n");
62     return 0;
63 }
并查集 贪心
posted @ 2014-08-03 13:48  Asm.Definer  阅读(364)  评论(0编辑  收藏  举报