POJ 2114 树的分治

这题跟1741的不同之处在于查找==k的是否存在 修改一下统计函数即可

#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <cstdio>
#include <cmath>
#include <queue>
#include <map>
#include <set>
#define eps 1e-5
#define MAXN 11111
#define MAXM 55555
#define INF 1000000000
using namespace std;
struct EDGE
{
    int v, next, w;
}edge[MAXM];
int head[MAXN], e;
int n, k, vis[MAXN], ans, root, num;
void init()
{
    memset(head, -1, sizeof(head));
    e = 0;
}
void add(int u, int v, int w)
{
    edge[e].v = v;
    edge[e].w = w;
    edge[e].next = head[u];
    head[u] = e++;
}
int mx[MAXN], size[MAXN], mi, dis[MAXN];
void dfssize(int u, int fa) //处理子树的大小
{
    size[u] = 1;
    mx[u] = 0;
    for(int i = head[u]; i != -1; i = edge[i].next)
    {
        int v = edge[i].v;
        if(v != fa && !vis[v])
        {
            dfssize(v, u);
            size[u] += size[v];
            if(size[v] > mx[u]) mx[u] = size[v];
        }
    }
}
void dfsroot(int r, int u, int fa) //求重心
{
    if(size[r] - size[u] > mx[u]) mx[u] = size[r] - size[u];
    if(mx[u] < mi) mi = mx[u], root = u;
    for(int i = head[u]; i != -1; i = edge[i].next)
    {
        int v = edge[i].v;
        if(v != fa && !vis[v]) dfsroot(r, v, u);
    }
}
void dfsdis(int u, int d, int fa) //求距离
{
    dis[num++] = d;
    for(int i = head[u]; i != -1; i = edge[i].next)
    {
        int v = edge[i].v;
        if(v != fa && !vis[v]) dfsdis(v, d + edge[i].w, u);
    }
}
int calc(int u, int d)
{
    int ret = 0;
    num = 0;
    dfsdis(u, d, 0);
    sort(dis, dis + num);
    int i = 0, j = num - 1;
    while(i < j) //经典
    {
        if(dis[i] + dis[j] < k) i++;
        else if(dis[i] + dis[j] > k) j--;
        else
        {
            if(dis[i] == dis[j])
            {
                ret += (j - i) * (j - i + 1) / 2;
                break;
            }
            int st = i, ed = j;
            while(dis[st] == dis[i]) st++;
            while(dis[ed] == dis[j]) ed--;
            ret += (st - i) * (j - ed);
            i = st, j = ed;
        }
    }
    return ret;
}
void dfs(int u)
{
    mi = n;
    dfssize(u, 0);
    dfsroot(u, u, 0);
    ans += calc(root, 0);
    vis[root] = 1;
    for(int i = head[root]; i != -1; i = edge[i].next)
    {
        int v = edge[i].v;
        if(!vis[v])
        {
            ans -= calc(v, edge[i].w);
            dfs(v);
        }
    }
}
void solve()
{
    memset(vis, 0, sizeof(vis));
    ans = 0;
    dfs(1);
    if(ans > 0) printf("AYE\n");
    else printf("NAY\n");
}
int main()
{
    while(scanf("%d", &n) != EOF && n)
    {
        int v, w;
        init();
        for(int i = 1; i <= n; i++)
        {
            while(scanf("%d", &v) != EOF && v)
            {
                scanf("%d", &w);
                add(i, v, w);
                add(v, i, w);
            }
        }
        while(scanf("%d", &k) != EOF && k) solve();
        printf(".\n");
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值