jzoj 5898. 【NOIP2018模拟10.6】距离统计 动态树分治

Description在这里插入图片描述
Input
在这里插入图片描述
Output
在这里插入图片描述
Sample Input

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

Sample Output

3
6
2

Data Constraint
jzoj 5898. 【NOIP2018模拟10.6】距离统计 动态树分治_第1张图片

分析:
我们可以先二分一个答案,那么问题就变成了从一个点 x x x开始的路径长度 ≤ m i d ≤mid mid的有多少条。
考虑用动态树分治。
每个点维护一个vector表示以这个点开始的路径的长度,并排序。还维护一些vector表示他的每一棵子树的到他路径上的长度,并排序。因为从 x x x开始,那么 x x x所在子树是没有贡献的,可以用总贡献减去这棵子树的贡献。然后再维护一个map表示每个点属于第几棵子树。
距离可以用倍增lca求。
总复杂度是 O ( n l o g 3 n ) O(nlog^3n) O(nlog3n)的。
好像数据有锅,所有人都是85分。

代码:

#include 
#include 
#include 
#include 
#include 
#include 
#include 

const int maxn=5e4+7;

using namespace std;

int n,m,x,y,w,cnt,root,sum,k;
int ls[maxn];
int size[maxn],height[maxn],p[maxn],vis[maxn],dis[maxn];
int f[maxn][16],dep[maxn],d[maxn];

vector  sondis[maxn*20];

struct node{
	vector  dis;
	map  num;
	int fa;
	int find(int x,int len)
	{
		int k1=upper_bound(dis.begin(),dis.end(),len)-dis.begin();
		int k2=upper_bound(sondis[x].begin(),sondis[x].end(),len)-sondis[x].begin();
		return k1-k2;
	}
}t[maxn];

struct edge{
	int y,w,next;
}g[maxn*2];

queue  q;

void add(int x,int y,int w)
{
	g[++cnt]=(edge){y,w,ls[x]};
	ls[x]=cnt;
}

void solve(int x,int fa)
{
	dep[x]=dep[fa]+1;
	f[x][0]=fa;
	for (int i=ls[x];i>0;i=g[i].next)
	{
		int y=g[i].y;
		if (y==fa) continue;
		d[y]=d[x]+g[i].w;
		solve(y,x);
	}
}

int getlca(int x,int y)
{
    if (dep[x]>dep[y]) swap(x,y);
    int d=dep[y]-dep[x],k=15,t=1<=t) d-=t,y=f[y][k];
        t/=2,k--;
    }
    if (x==y) return x;
    k=15;
    while (k>=0)
    {
        if (f[x][k]!=f[y][k])
        {
            x=f[x][k];
            y=f[y][k];
        }
        k--;
    }
    return f[x][0];
}

int getdis(int x,int y)
{
	int lca=getlca(x,y);
	return d[x]+d[y]-d[lca]*2;
} 

void findroot(int x,int fa)
{
    size[x]=1;
    height[x]=0;
    for (int i=ls[x];i>0;i=g[i].next)
    {
        int y=g[i].y;
        if ((y==fa) || (vis[y])) continue;
        findroot(y,x);
        size[x]+=size[y];
        height[x]=max(height[x],size[y]);
    }
    height[x]=max(height[x],sum-size[x]);
    if ((height[x]0;i=g[i].next)
	{
		int y=g[i].y;
		if ((y==fa) || (vis[y])) continue;
		dis[y]=dis[x]+g[i].w;
		t[root].num[y]=cnt;
		t[root].dis.push_back(dis[y]);
	    sondis[cnt].push_back(dis[y]);
		dfs(y,x);
		size[x]+=size[y];
		if (x==root)
		{
			p[y]=x;
			q.push(y);
			sort(sondis[cnt].begin(),sondis[cnt].end());
			cnt++;
		}
	}
}

int check(int x,int len)
{
	int sum=0;
	for (int i=x;i>0;i=t[i].fa)
	{
		int d=t[i].num[x];
		sum+=t[i].find(d,len-getdis(x,i));
	}
	return sum-1;
}

int main()
{
	freopen("tree.in","r",stdin);
	freopen("tree.out","w",stdout);
	scanf("%d%d",&n,&m);
	for (int i=1;i=k) ans=mid,r=mid-1;
			                else l=mid+1;
		}
		printf("%d\n",ans);
	}
}

你可能感兴趣的:(点分治)