【BZOJ5100】[POI2018]Plan metra 构造
构造 plan
2023-09-11 14:15:24 时间
【BZOJ5100】[POI2018]Plan metra
Description
有一棵n个点的无根树，每条边有一个正整数权值，表示长度，定义两点距离为在树上的最短路径的长度。
已知2到n-1每个点在树上与1和n的距离，请根据这些信息还原出这棵树。
Input
第一行包含一个正整数n(2<=n<=500000)，表示点数。
第二行包含n-2个正整数d(1,2),d(1,3),...,d(1,n-1)，分别表示每个点到1的距离。
第三行包含n-2个正整数d(n,2),d(n,3),...,d(n,n-1)，分别表示每个点到n的距离。
输入数据保证1<=d<=1000000。
Output
若无解，输出NIE。
否则第一行输出TAK，接下来n-1行每行三个正整数u,v,c(1<=u,v<=n,1<=c<=1000000)
表示存在一条长度为c的连接u和v两点的树边。
若有多组解，输出任意一组。
Sample Input
7
6 6 2 2 1
5 3 5 1 4
6 6 2 2 1
5 3 5 1 4
Sample Output
TAK
1 5 2
5 7 1
5 2 4
7 3 3
1 4 2
1 6 1
1 5 2
5 7 1
5 2 4
7 3 3
1 4 2
1 6 1
题解：如果我们已经确定了点1到点n的距离，那么我们就可以将1到n这条路拿出来，然后其他点都想挂链一样挂到这条路径上即可。即：
如果1到n的长度是m，我们将所有d1-dn相同的点放到一起，那么对于每组d1-dn相同的点，要么|d1-dn|=m，要么d1+dn的最小值=m。所以用数组记录一下即可。
输出方案的时候将所有d1+dn=m的点排个序即可。
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int N=500010;
const int M=1000010;
int n,m,tot,cnt;
int d1[N],d2[N],s1[M<<1],s2[M<<1],p[M],mn[M<<1],pa[N],pb[N],pc[N];
inline char nc()
{
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int rd()
{
int ret=0,f=1; char gc=nc();
while(!isdigit(gc)) {if(gc=='-') f=-f; gc=nc();}
while(isdigit(gc)) ret=ret*10+(gc^'0'),gc=nc();
return ret*f;
}
inline int abs(int x) {return x>0?x:-x;}
inline bool check(int x)
{
if(((x>=M||x<=-M)?0:(s1[M+x]+s1[M-x]))+s2[x]==n-2)
{
m=x;
return 1;
}
return 0;
}
inline void add(int a,int b,int c)
{
pa[++cnt]=a,pb[cnt]=b,pc[cnt]=c;
}
int main()
{
n=rd();
if(n==2)
{
printf("TAK\n1 2 1");
return 0;
}
register int i,last;
for(i=2;i<n;i++) d1[i]=rd();
for(i=2;i<n;i++)
{
d2[i]=rd();
int &t=mn[d2[i]-d1[i]+M];
if(!t||d1[t]+d2[t]>d1[i]+d2[i]) t=i;
s1[d2[i]-d1[i]+M]++;
}
for(i=2;i<n;i++) if(i==mn[d2[i]-d1[i]+M]) s2[d1[i]+d2[i]]+=s1[d2[i]-d1[i]+M];
for(i=2;i<n;i++) if(check(abs(d2[i]-d1[i]))||check(d1[i]+d2[i])) break;
if(!m)
{
puts("NIE");
return 0;
}
for(i=2;i<n;i++)
{
if(d2[i]-d1[i]==m) add(i,1,d1[i]);
else if(d1[i]-d2[i]==m) add(i,n,d2[i]);
else if(i==mn[d2[i]-d1[i]+M]) p[d1[i]]=i;
else
{
int t=mn[d2[i]-d1[i]+M];
add(i,t,(d1[i]+d2[i]-d1[t]-d2[t])>>1);
}
}
p[m]=n,d1[n]=m;
for(last=i=1;i<=m;i++) if(p[i]) add(last,p[i],d1[p[i]]-d1[last]),last=p[i];
for(i=1;i<n;i++) if(pc[i]<=0)
{
puts("NIE");
return 0;
}
puts("TAK");
for(i=1;i<n;i++) printf("%d %d %d\n",pa[i],pb[i],pc[i]);
return 0;
}//6 2 5 4 8 4 1 6 4