问题 B: Fiolki
时间限制: 3 Sec 内存限制: 128 MB
题目描述
化学家吉丽想要配置一种神奇的药水来拯救世界。
吉丽有n种不同的液体物质，和n个药瓶（均从1到n编号）。初始时，第i个瓶内装着g[i]克的第i种物质。吉丽需要执行一定的步骤来配置药水，第i个步骤是将第a[i]个瓶子内的所有液体倒入第b[i]个瓶子，此后第a[i]个瓶子不会再被用到。瓶子的容量可以视作是无限的。
吉丽知道某几对液体物质在一起时会发生反应产生沉淀，具体反应是1克c[i]物质和1克d[i]物质生成2克沉淀，一直进行直到某一反应物耗尽。生成的沉淀不会和任何物质反应。当有多于一对可以发生反应的物质在一起时，吉丽知道它们的反应顺序。每次倾倒完后，吉丽会等到反应结束后再执行下一步骤。
吉丽想知道配置过程中总共产生多少沉淀。
输入
第一行三个整数n,m,k(0<=m<n<=200000,0<=k<=500000)，分别表示药瓶的个数（即物质的种数），操作步数，可以发生的反应数量。
第二行有n个整数g[1],g[2],…,g[n]（1<=g[i]<=10^9)，表示初始时每个瓶内物质的质量。
接下来m行，每行两个整数a[i],b[i](1<=a[i],b[i]<=n,a[i]≠b[i])，表示第i个步骤。保证a[i]在以后的步骤中不再出现。
接下来k行，每行是一对可以发生反应的物质c[i],d[i](1<=c[i],d[i]<=n,c[i]≠d[i])，按照反应的优先顺序给出。同一个反应不会重复出现。
输出
配置过程中总共产生多少沉淀。
样例输入
3 2 1 2 3 4 1 2 3 2 2 3
样例输出
6 乱搞压正解，暴力出奇迹。 数据水到一定地步，当时本来期望打50分的暴力A了……
1 #include<iostream>
2 #include<cstdlib>
3 #include<cstdio>
4 #include<cstring>
5 #include<queue>
6 #include<algorithm>
7 #include<cmath>
8 #define N 200005
9 using namespace std;
10 int n,m,t,zz2;
11 int g[N],a[N],b[N];
12 int f[N][2],zz;
13 struct ro{
14 int to,l,from;
15 int next;
16 bool friend operator > (ro a,ro b)
17 {
18 return a.l>b.l;
19 }
20 }road[1000005],road2[1000006];
21 void build(int x,int y,int z){
22 zz++;
23 road[zz].to=y;
24 road[zz].from=x;
25 road[zz].l=z;
26 road[zz].next=a[x];
27 a[x]=zz;
28 }
29 void build2(int x,int y){
30 zz2++;
31 road2[zz2].from=x;
32 road2[zz2].to=y;
33 road2[zz2].next=b[x];
34 b[x]=zz2;
35 }
36 int fa[N];
37 int find(int x){
38 if(fa[x]==x)return x;
39 fa[x]=find(fa[x]);
40 return fa[x];
41 }
42 long long hav[N];
43 void hb(int a,int b){
44 int x=find(a);
45 int y=find(b);
46 fa[x]=y;
47 hav[y]+=hav[x];
48 build2(y,x);
49 }
50 bool cle[N],cle2[N];
51 priority_queue<ro,vector<ro>,greater<ro > > q1;
52 void dfs(int x,int tt){
53 //cout<<x<<endl
54 if(!cle2[x])
55 {
56 for(int i=a[x];i>0;i=road[i].next)
57 {
58 int y=road[i].to;
59 if(find(y)==tt)
60 {
61 q1.push(road[i]);
62 }
63 }
64 }
65 bool yx=1;
66 for(int i=b[x];i>0;i=road2[i].next)
67 {
68 int y=road2[i].to;
69 if(!cle[y])
70 {
71 //yx=0;
72 dfs(y,tt);
73 }
74 if(!cle[y])
75 {
76 yx=0;
77 }
78 }
79 cle[x]=yx&cle2[x];
80 }
81 int main(){
82 scanf("%d%d%d",&n,&m,&t);
83 for(int i=1;i<=n;i++)
84 scanf("%d",&g[i]);
85 for(int i=1;i<=m;i++)
86 scanf("%d%d",&f[i][0],&f[i][1]);
87 for(int i=1;i<=t;i++)
88 {
89 int x,y;
90 scanf("%d%d",&x,&y);
91 build(x,y,i);
92 build(y,x,i);
93 }
94 for(int i=1;i<=n;i++)
95 fa[i]=i;
96 for(int i=1;i<=m;i++)
97 {
98 int from=f[i][0],to=f[i][1];
99 dfs(from,to);
100 while(!q1.empty())
101 {
102 ro tt=q1.top();
103 q1.pop();
104 int x=tt.from,y=tt.to;
105 int p=min(g[x],g[y]);
106 g[x]-=p;
107 g[y]-=p;
108 hav[to]+=p*2;
109 if(!g[x]) cle2[x]=1;
110 if(!g[y]) cle2[y]=1;
111 }
112 hb(from,to);
113 }
114 long long sum=0;
115 for(int i=1;i<=n;i++)
116 {
117 if(find(i)==i)
118 {
119 sum+=hav[i];
120 }
121 }
122 printf("%lld\n",sum);
123 //while(1);
124 return 0;
125 }
时间： 2024-11-09 03:21:39