Time Limit: 5 Sec  Memory Limit: 256 MBSec  Special Judge
Submit: 293  Solved: 144
[
Submit][Status][Discuss]

Description

给定一棵有n个节点的树,相邻两点之间的距离为1。

请找到一个点x,使其满足所有m条限制,其中第i条限制为dist(x,a[i])+dist(x,b[i])<=d[i]。

Input

第一行包含一个正整数t(1<=t<=1000),表示数据组数。

对于每组数据,第一行包含两个正整数n,m(2<=n,m<=300000),表示点数、限制数。

接下来n-1行,每行两个正整数x,y(1<=x,y<=n),表示树上的一条边。

接下来m行,每行三个正整数a[i],b[i],d[i](1<=a[i],b[i]<=n,1<=d[i]<=600000),描述一条限制。

输入数据保证所有n之和不超过300000,所有m之和也不超过300000。

Output

输出t行。第i行输出第i组数据的答案,如果无解输出NIE,否则输出TAK,

然后输出x,如有多组解,输出任意一组。

Sample Input

2
5 3
1 2
2 3
2 4
3 5
1 4 2
5 5 5
3 2 1
3 2
1 2
2 3
1 1 2
3 3 1

Sample Output

TAK 2
NIE

HINT

Source

鸣谢Claris上传

题解:

        ①首先呢要搞懂:(设1为根)

               满足第i个限制条件:dist(x,a[i])+dist(x,b[i])<=d[i] 的所有点中,

               离根节点最近的点的到达根的距离D是有一个公式的: 

                D=max(   0   ,   (dis(a[i],1) + dis(b[i],1) - d[i]) / 2    )

                [温馨提示:画个图走一波就可以懂]

        ②然后呢,根据上文结论呢,设 D 值 最大的点时 P。

               有一个更高级的结论:如果P都不能合法,那么没有什么点会合法了。

        ③插几句,怎么找到P点:

                      根据P点必须满足最苛刻的条件,可以从对应的a[i],b[i]为根遍历树得到dis[],

                      然后呢找到满足最苛刻条件,且到1最近的那个点就是我们D对应的P点了。

        ④我是这样理解这个结论的:

                      (1)P点的位置是由最苛刻的限制条件得出的(即离根最远)

                      (2)对于和它统一棵子树的限制条件,肯定他可以满足所有(因为他是最苛刻的)

                      (3)对于和它不在同一子树的限制条件,也必须想办法满足它才可以。

#include<cstdio>
#include<cstring>
#define go(i,a,b) for(int i=a;i<=b;i++)
#define mem(a,b) memset(a,b,sizeof(a))
#define fo(i,a,x) for(int i=a[x],v=e[i].v;i;i=e[i].next,v=e[i].v)
const int N=400010;
struct E{int v,next;}e[N<<1];
int T,n,m,i,x,y,t,ed,a[N],b[N],d[N],D[3][N],k,head[N];
void ADD(int u,int v){e[k]=(E){v,head[u]};head[u]=k++;}
void dfs(int u,int fa,int z,int p)
{
D[p][x]=z++;
fo(i,head,u)if(v!=fa)dfs(v,u,z,p);
}
int main()
{
scanf("%d",&T);
while(T--&&scanf("%d%d",&n,&m))
{
go(i,1,n)head[i]=0;k=1;
go(i,2,n){int u,v;scanf("%d%d",&u,&v),ADD(u,v),ADD(v,u);}
go(i,1,m)scanf("%d%d%d",a+i,b+i,d+i);
dfs(1,0,0,0);x=0;
go(i,1,m)
{
t=D[0][a[i]]+D[0][b[i]]-d[i];
if(!x||t>y)x=i,y=t;
}
dfs(a[x],0,0,1),dfs(b[x],0,0,2);t=0;
go(i,1,n)if(D[1][i]+D[2][i]<=d[x])if(!t||D[0][i]<y)t=i,y=D[0][i];
if(!t){puts("NIE");continue;}dfs(t,0,0,0);
go(i,1,m)if(D[0][a[i]]+D[0][b[i]]>d[i]){t=0;break;}
if(t)printf("TAK %d\n",t);else puts("NIE");
}
return 0;
}//Paul_Guderian

外面下起了小雨,雨滴轻飘飘得像我的年轻岁月,

我脸上蒙着雨水就像蒙着幸福……——————————汪峰《青春》

【BZOJ 4151 The Cave】的更多相关文章

  1. 【Bzoj 1835 基站选址】

    基站选址的区间里隐藏着DP优化的机密…… 分析:       不论是做过乘积最大还是石子合并,或者是其他的入门级别的区间DP题目的人呐,大米并认为读题后就能够轻松得出一个简洁明了的Dp转移方程.    ...

  2. 【BZOJ 2744 朋友圈】

    Time Limit: 30 Sec  Memory Limit: 128 MBSubmit: 1570  Solved: 532[Submit][Status][Discuss] Descripti ...

  3. 【BZOJ 5038 不打兔子】

    Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 22  Solved: 8[Submit][Status][Discuss] Description 勤 ...

  4. 【BZOJ 1088 扫雷Mine】模拟

    http://www.lydsy.com/JudgeOnline/problem.php?id=1088 2*N的扫雷棋盘,第二列的值a[i]记录第 i 个格子和它8连通的格子里面雷的数目. 第一列的 ...

  5. 【BZOJ做题记录】07.07~?

    在NOI一周前重开一个坑 最后更新时间:7.08 07:38 7.06 下午做的几道CQOI题: BZOJ1257: [CQOI2007]余数之和sum:把k mod i写成k-k/i*i然后分段求后 ...

  6. 【bzoj5050】【bzoj九月月赛H】建造摩天楼

    讲个笑话,这个题很休闲的. 大概是这样的,昨天看到这个题,第一眼星际把题目看反了然后感觉这是个傻逼题. 后来发现不对,这个修改一次的影响是很多的,可能导致一个数突然可以被改,也可能导致一个数不能被改. ...

  7. 【BZOJ 2458 最小三角形】

    Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1551  Solved: 549[Submit][Status][Discuss] Descripti ...

  8. 【BZOJ 5000 OI树】

    Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 107  Solved: 64[Submit][Status][Discuss] Description ...

  9. 【BZOJ 5047 空间传送装置】

    Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 282  Solved: 121[Submit][Status][Discuss] Descriptio ...

随机推荐

  1. 【例题收藏】◇例题·III◇ 木と整数 / Integers on a Tree

    ◇例题·III◇ 木と整数 / Integers on a Tree 只需要一个美妙的转换,这道题就会变得无比美妙…… 来源:+AtCoder 2148(ARC-063 E)+ ◆ 题目大意 给定一棵 ...

  2. matlab2018a安装后帮助文档打不开解决方法

    安装matlab2018a破解版后,帮助文档提示需要许可证问题(破解版没有可用许可证): 解决方法是把文档设置为离线即可(预设---->帮助---->安装在本地---->小窗口)

  3. selenium webdriver 移动到某个位置,滚动到某个位置

    https://blog.csdn.net/xiaosongbk/article/details/70231564

  4. AWS CentOS7 实例修改主机名

    问题描述: AWS EC2 实例在升级到CentOS7以后,我们发现主机名的修改不再像之前的版本(CentOS 5/6)一样简单. 每次新建实例之后,修改好主机名,重启或者克隆之后的机器,主机名还是会 ...

  5. B1086 就不告诉你 (15分)

    B1086 就不告诉你 (15分) 做作业的时候,邻座的小盆友问你:"五乘以七等于多少?"你应该不失礼貌地围笑着告诉他:"五十三."本题就要求你,对任何一对给定 ...

  6. POJ:2674-Linear world(名字交换碰撞)

    Linear world Time Limit: 3000MS Memory Limit: 65536K Total Submissions: 4514 Accepted: 1025 Descript ...

  7. [学习笔记]CSS选择器

    CSS语法结构 selector {     property1 : value;     property2 : value; } 如果包含多个属性,那么属性每个属性的结尾需要有一个分号.如果属性的 ...

  8. Python locale 多语言模块和我遇到的坑

    Table of Contents 1. locale遇到的问题 1.1. locale 简介 1.1.1. 什么是locale 1.1.2. locale 相关命令 1.2. Python loca ...

  9. 减少Android staido 占用C 盘

    1.gradle 更换文件夹: 设置GRADLE_USER_HOME环境变量 在/etc/profile或~/.bash_profile增加如下: export GRADLE_USER_HOME=D: ...

  10. 手把手教你如何逐步安装OpenStack

    [TechTarget中国原创] 尽管OpenStack官方提供的在线安装教程和分步向导能够为管理员提供很大帮助,但是依然存在很多不尽如人意的地方.因此在Ubuntu上安装OpenStack的过程当中 ...