JZOJ-senior-5916. 【NOIP2018模拟10.20】flow

Time Limits: 1000 ms Memory Limits: 262144 KB

Description

你是申国的一个地方长官，你手下有n个城市。

为了加强基础设施建设，在2020全面建成小康社会，统筹推进经济建设、政治建设、文化建设、社会建设、生态文明建设，坚定实施科教兴国战略、人才强国战略、创新驱动发展战略、乡村振兴战略、区域协调发展战略、可持续发展战略、军民融合发展战略，突出抓重点、补短板、强弱项，特别是要坚决打好防范化解重大风险、精准脱贫、污染防治的攻坚战，使全面建成小康社会得到人民认可、经得起历史检验。你认为本省的水利调配非常有问题，这导致部分地区出现严重的缺水，而部分地区却全年洪灾泛滥。

于是你打算将原有的但是已经废弃了的m条水管重新使用。第i条水管连接城市xi和yi。这些水管联通了所有城市。每座城市对水的需求不同设为ai，部分城市处于缺水状态，ai为正，缺水量刚好为ai mol。部分城市因为有水库，ai为负，它需要向外输送-ai mol的水才能不形成洪灾。对于每条水管，你需要决定它的输送量fi，若fi为正则表示从xi向yi输送fi mol的水，fi为负则表示从yi向xi输送-fi mol的水。

你需要做到每个城市都刚好满足它的需求，即缺ai mol水的城市需要刚好输入ai的水，而多出-ai mol水的城市需要刚好输出-ai mol水。

你需要判断能否满足要求，若满足，你还需要输出所有的f。

Input

第一行一个整数n。
第二行n个整数ai。
第三行一个整数m。
接下来m行每行两个整数，xi,yi。
意义如题中所述。

Output

如果不能满足要求，输出”Impossible”。
否则第一行输出”Possible”，接下来m行，第i+1行输出fi。
若有多组解，随意输出一组即可。

Sample Input

4
3 -10 6 1
5
1 2
3 2
2 4
3 4
3 1

Sample Output

Possible
4
-6
8
-7
7

Data Constraint

对于10%：n<=5，|ai|<=20
对于20%：n<=10，|ai|<=20
对于另30%：n<=1000，m=n，|ai|<=20
对于100%：n<=200000，m<=300000，你使用longint(pascal)/int(c++)就能够存储ai。
数据纯随机。

Hint

Tips
mol是一种单位。

Solution

显然仅 $a$ 的和为 $0$ 时有解
我们发现，在有解的情况下，只要这些水管将所有的城市连通了就一定就可以保证有方案满足题意，于是我们搞出一棵生成树，对于 $x$，连接 $x$ 和它父亲节点的边权就是 $x$ 的子树和，注意判断一下方向来确定正负

Code

#include<algorithm>
#include<cstdio>

#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 ll long long

using namespace std;

const int N=3e5+5;
int n,m,num;
int b[N],f[N],v[N],sz[N],last[N];
ll s,g[N],ans[N];

struct node
{
	int x,y;
}a[N];

struct edge
{
	int to,next,id,k;
}e[2*N];

int get(int x)
{
	return f[x]=f[x]==x?x:get(f[x]);
}

void merge(int x,int y)
{
	f[get(x)]=get(y);
}

void link(int x,int y,int z,int w)
{
	e[++num]=(edge){y,last[x],z,w},last[x]=num;
}

void dfs(int x)
{
	v[x]=1;
	for(int w=last[x];w;w=e[w].next)
	{
		int y=e[w].to;
		if(v[y]) continue;
		dfs(y);
		g[x]+=g[y];
		int id=e[w].id;
		ans[id]=e[w].k*g[y];
	}
}

int main()
{
	freopen("flow.in","r",stdin);
	freopen("flow.out","w",stdout);
	scanf("%d",&n);
	fo(i,1,n)
	{
		scanf("%d",&b[i]),s+=b[i];
		f[i]=i,sz[i]=1,g[i]=b[i];
	}
	if(s)
	{
		printf("Impossible\n");
		return 0;
	}
	else printf("Possible\n");
	scanf("%d",&m);
	fo(i,1,m)
	{
		scanf("%d%d",&a[i].x,&a[i].y);
		if(get(a[i].x)!=get(a[i].y))
		{
			merge(a[i].x,a[i].y);
			link(a[i].x,a[i].y,i,1);
			link(a[i].y,a[i].x,i,-1);
		}
	}
	dfs(1);
	fo(i,1,m) printf("%lld\n",ans[i]);
}