码迷,mamicode.com
首页 > 其他好文 > 详细

POJ1741:Tree——点分治

时间:2019-10-01 10:00:08      阅读:82      评论:0      收藏:0      [点我收藏+]

标签:isp   暴力   getch   个数   时间   bool   sed   bre   暴力枚举   

题面

  POJ1741

解析

   刚学了点分治,练一练模版题

  过程就不多说了,主要说说细节

  在每次查询下一棵子树时, 传进去的整棵子树大小是上一次的$siz$, 这个数据其实是错的, 但好像并不影响时间复杂度, 这样的话找重心就必须找最大子树最小的点了,否则会错。因此需要存一个当前最大子树最小的点的最大子树的大小, 以及当前的重心, 每次找重心之前,前者要赋为$inf$

  在每次统计以当前点为$lca$的链的答案时, 要先把这个重心加入数组中, 再按$dis$排序, 然后用双指针法或是其他可行的方法统计答案, 不能暴力枚举点对,不然复杂度不对

 代码:

技术图片
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn = 40005, inf = 0x3f3f3f3f;

inline int read()
{
    int ret, f=1;
    char c;
    while((c=getchar())&&(c<0||c>9))if(c==-)f=-1;
    ret=c-0;
    while((c=getchar())&&(c>=0&&c<=9))ret=(ret<<3)+(ret<<1)+c-0;
    return ret*f;
}

int n, k;

int head[maxn], tot;
struct edge{
    int nxt, to, d;
}e[maxn<<1];

void Addedge(int x, int y, int w)
{
    e[++tot] = (edge){head[x], y, w};
    head[x] = tot;
}

int cnt;
int dis[maxn], siz[maxn], bel[maxn], p[maxn], ans, val, rt;
bool vis[maxn];

void dfs(int x, int fa, int sum)
{
    int ret = 0;    
    siz[x] = 1;
    for(int i = head[x]; i; i = e[i].nxt)
    {
        int id = e[i].to;
        if(vis[id] || id == fa)    continue;
        dfs(id, x, sum);
        siz[x] += siz[id];
        ret = max(ret, siz[id]);
    }
    ret = max(ret, sum - siz[x]);
    if(ret < val)
    {
        val = ret;
        rt = x;
    }
}

bool cmp(int x, int y)
{
    return dis[x] < dis[y];
}

void get_dis(int x, int fa, int b)
{
    p[++cnt] = x;
    bel[x] = b;
    for(int i = head[x]; i; i = e[i].nxt)
    {
        int id = e[i].to;
        if(id == fa || vis[id])    continue;
        dis[id] = dis[x] + e[i].d;
        get_dis(id, x, b);
    }
}

int num[maxn];

void calc(int x)
{
    cnt = 0;    
    p[++cnt] = x;
    dis[x] = 0;
    bel[x] = x;
    for(int i = head[x]; i; i = e[i].nxt)
    {
        int id = e[i].to;
        if(vis[id])    continue;
        dis[id] = e[i].d;
        get_dis(id, x, id);
    }
    sort(p + 1, p + cnt + 1, cmp);    
    for(int i = 1; i <= cnt; ++i)
        ++ num[bel[p[i]]];
    int l = 0, r = cnt;
    while(l < r - 1)
    {
        ++ l;
        -- num[bel[p[l]]];
        while(dis[p[r]] + dis[p[l]] > k)
        {
            -- num[bel[p[r]]];
            -- r;
            if(r <= l)
                break;
        }
        if(r <= l)    break;
        ans += r - l - num[bel[p[l]]];    
    }
    if(r > l)
        -- num[bel[p[r]]];
}

void solve(int x)
{
    vis[x] = 1;
    calc(x);
    for(int i = head[x]; i; i = e[i].nxt)
    {
        int id = e[i].to;
        if(vis[id])    continue;
        val = inf;
        dfs(id, x, siz[id]);
        solve(rt);
    }
}

int main()
{
    while(scanf("%d%d", &n, &k), n != 0)
    {
        for(int i = 1; i < n; ++i)
        {
            int u = read(), v = read(), w = read();
            Addedge(u, v, w);
            Addedge(v, u, w);
        }
        val = inf;
        dfs(1, 0, n);
        solve(rt);
        printf("%d\n", ans);
        ans = tot = 0;
        for(int i = 1; i <= n; ++i)
            head[i] = vis[i] = 0;
    }
    return 0;
}
View Code

 

POJ1741:Tree——点分治

标签:isp   暴力   getch   个数   时间   bool   sed   bre   暴力枚举   

原文地址:https://www.cnblogs.com/Joker-Yza/p/11614532.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!