【广东工业大学程序设计竞赛决赛】tmk找三角

个人感觉这道题还是很有营养的，YY能力得到大幅度提高。

题目链接：http://gdutcode.sinaapp.com/problem.php?cid=1057&pid=5

楼主弱鸡是当时跟着一起搞这场比赛的，1个半钟边刷榜边过了4题。。他们服务器太弱，跑题跑太慢了。。就没心情搞下去了。。

然后开始YY这道题，用了各种YY功力都没有yy过他的巨大数据应该如何处理。。直到如今才脑洞大开。

---

题意：题目上是说，在一个图上，某东西需要走最短路，然后问他要从s走到t，问走的这条路线上的所有边能否有三条边构成三角形。其实这个人走的路都是套路。因为图是一个生成树，所以，无论怎么走都是最短路。。。

然后这题套路深就深在，他的点很多1e5，边值看起来并没有什么用1e9，其实。。这东西才是关键，对于1e9范围内，如果算去边，最多能剩下多少条边使他们不能构成三角形。这道题看起来就像ccpc长春现场赛的一道题。。当时有脑洞。。现在没脑洞了。。想要构成最多条不能构成三角形的一个区间，只需要区间内所有数都满足斐波那契序列即可，因为三角形两边之和一定大于第三边，所以满足斐波那契序列的一定构不成三角形。

所以，对于1e9以内，可以构成多少不能构成三角形的。。好像50条是最多了。。所以。只要超过50条边，他们一定能构成三角形（我列个擦哦）。所以这个东西可以完美解决常数的访问量1e5，每次访问我们最多只用走50条边。

然后这个考虑完了之后，我们还需要考虑一个黑科技，就是如何才能最快的在一棵树内找到两点的所有边。这个对于一个节点无规律的树来说是很有难度的。但是对于一棵满足最大堆的树来说，是很简单的，只要每一层的数都比上一层的数大，我们访问的时候就能够通过比较两个节点当前点的权值，来判断他们谁应该向上移动，这样的访问速度为他们之间的边的条数，而他们的边又只有50条，所以更加完美解决这道题的数据。

然而如何才能构造出这样子的树呢，其实只需要bfs一次，从1这个结点往下一层一层给id，保证上一层的id小于当前层任意一个点的id即可（由于我比较弱，只能想到这种写法，若有大佬有其他写法请指教）

本人的构图能力比较差，可能写出来有垃圾代码。。各位多见谅，一定虚心学习构图能力。。毕竟边点的构建我删了又写了好几次。。。

附上弱弱的代码：

#include <bits/stdc++.h>

#define MAXN 100005

using namespace std;

struct Edge{
    int from,to,weight;
    Edge(int _from = 0,int _to = 0,int _weight = 0) : from(_from),to(_to),weight(_weight){}
};

struct Point{
    int pre,index,u;
    vector<Edge> v;
}p[MAXN];

int dis[MAXN];

void addEdge(int u,int v,int w){
    p[u].v.push_back(Edge(u,v,w));
}

void bfs(){
    queue<Point> q;
    q.push(p[1]);
    int idx = 1;
    p[1].index = idx;
    while(!q.empty()){
        Point tmp = q.front();
        q.pop();
        for(int i = 0;i < tmp.v.size();i++){
            if(tmp.v[i].to == tmp.pre){
                continue;
            }
            p[tmp.v[i].to].pre = tmp.u;
            p[tmp.v[i].to].index = ++idx;
            dis[tmp.v[i].to] = tmp.v[i].weight;
            q.push(p[tmp.v[i].to]);
        }
    }
}

void solve(int a,int b){
    vector<int> vs;
    int x = p[a].index,y = p[b].index;
    while(vs.size() < 50 && x != y){
        if(x > y){
            x = p[p[a].pre].index;
            vs.push_back(dis[a]);
            a = p[a].pre;
        } else {
            y = p[p[b].pre].index;
            vs.push_back(dis[b]);
            b = p[b].pre;
        }
    }
    bool flag = false;
    if(vs.size() == 50){
        printf("Yes\n");
        flag = true;
    } else {
        sort(vs.begin(),vs.end());
        for(int i = 2;i < vs.size();i++){
            if(vs[i - 2] + vs[i - 1] > vs[i]){
                printf("Yes\n");
                flag = true;
                break;
            }
        }
    }
    if(!flag){
        printf("No\n");
    }
}

void init(int n){
    for(int i = 0;i <= n;i++){
        p[i].v.clear();
    }
}

int main(){
    int T,a,b,u,v,w,n,q;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        init(n);
        for(int i = 1;i <= n;i++){
            p[i].u = i;
        }
        for(int i = 0;i < n - 1;i++){
            scanf("%d %d %d",&u,&v,&w);
            addEdge(u,v,w);
            addEdge(v,u,w);
        }
        bfs();
        scanf("%d",&q);
        while(q--){
            scanf("%d %d",&a,&b);
            solve(a,b);
        }
    }
}