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. js动画之无缝滚动

    效果图如下: HTML代码如下: <!DOCTYPE html> <html lang="en"> <head> <meta charse ...

  2. linux文件共享之samba,nfs的搭建

    Samba server简介 Samba是在Linux和UNIX系统上实现SMB协议的一个免费软件,由服务器及客户端程序构成.SMB(Server Messages Block,信息服务块)是一种在局 ...

  3. ethereum(以太坊)(三)--合约单继承与多继承

    pragma solidity ^0.4.0; // priveta public internal contract Test{ //defualt internal uint8 internal ...

  4. PHP开发环境搭建一:PHP集成环境XAMPP 的安装与配置

    1. XMAPP简介 XAMPP(Apache+MySQL/MariaDB+PHP+Perl)开头的X代表X-OS,代表可以在任何常见操作系统下使用,包括Windows.Mac.Linux,开源平台. ...

  5. 学习Pytbon第九天,函数1 过程和参数

    函数def func1():定义函数 '''testing1'''#函数的说明 print("in the func1")#定义过程 return 0 #得到函数的执行结果.还是程 ...

  6. linux最大进程数

    使用 ulimit -a 命令,查看 max user processes 的输出,就是系统最大进程数 core file size (blocks, -c) unlimited data seg s ...

  7. Android MultiType第三方库的基本使用和案例+DiffUtil的简单用法

    1.MultiType简单介绍 1.1.MultiType用于比较复杂的页面. 如下图,今日头条用到了MultiType处理各种复杂的页面.    这种还是比较简单的类型.因为一个页面也就这种类型. ...

  8. [bzoj1999][noip2007]Core树网的核

    好久没写题解了.这题不算太水就写一下题解. 话说回来,虽然不水但是挺裸.可以说题意即一半题解了. 我猜粘了题面也没有人去看的,所以直接人话题意了. 给一棵树,点数1e6,(当年noip的n当然是只有3 ...

  9. 多个".h"文件中声明及定义 全局变量和函数

    一.".h"文件必须以如下格式书写 例:文件<CZ_efg_hi.h"> ------------文件内容----------- #ifndef CZ_Efg ...

  10. 利用split方法计算字符串中出现字母最多的次数

    最近练习一些简单的算法题,知道自己很不聪明,但却没想到用了这么久,划算不划算是个需要考虑的问题, 其中有个算法是:统计一个字符串出现最多的字母,网上很多自己的见解,但是才疏学浅,有些地方看的有点困难, ...