借使它是水果路线的子路线)

4009: [Hnoi2015]接水果

Time Limit: 60 Sec  Memory Limit: 512 MB

Description

风见芬芳特别喜欢玩贰个号称osu!的玩耍,此中他最爱怜玩的格局正是接水果。由于他早就DT FC 了The big
black,
她以为那一个游乐太轻巧了,于是发明了二个特别难的本子。首先有贰个地形图,是一棵由
n 个极端、n-1 条边组成的树(例如图 1交给的树包罗 8 个顶峰、7
条边)。那颗树上有 P 个盘子,每一种盘子实际上是一条路线(比如图 1 中顶点 6
到终端 8 的不二等秘书技),何况每一个盘子还应该有一个权值。第 i
个盘子便是顶点aiai到顶点bibi的门路(由于是树,所以从aiai到bibi 的路线是独一的), 权值为cici。接下来依次会有Q个瓜果掉下来,每一个水果本质上也是一条门路,第i
个水果是从顶点 uiui 到终端vivi 的路径。 
菲菲每趟须求选拔贰个市场价格去接当前的果品:二个盘子能接住三个水果,当且仅当盘子的门路是鲜果的不二等秘书技的子路线(比方图第11中学从 3到7 的路子是从1到8的门路的子路线)。 
此地规定:从a 到b的渠道与从b到
a的路线是一致条路线。当然为了升高难度,对于第 i
个瓜果,你须要选取能接住它的有所盘子中,权值第 kiki 小的那么些盘子,每一个盘子可重复使用(没有行使次数的上限:一个增势接完三个水果后,前边还可承袭袭其余水果,只要它是水果路线的子路线)。芳香以为那几个游乐很难,你能轻轻易松化解给她看吗? 
图片 1

Input

先是行两个数 n和P
和Q,表示树的分寸和物价指数的个数和瓜果的个数。 
接下去n-1 行,每行七个数
a、b,表示树上的a和b 之间有一条边。树中顶点按1到 n标号。 
接下去 P 行,每行四个数
a、b、c,表示路线为 a 到 b、权值为 c 的物价指数,在那之中0≤c≤109109,a不对等b。 
接下去Q行,每行五个数 u、v、k,表示路线为 u到 v的果品,其中u不等于v,你要求选择第 k小的生势,第k
小一定期存款在。

Output

对于每种果子,输出一行表示选拔的盘子的权值。

Sample Input

10 10 10 
1 2 
2 3 
3 4 
4 5 
5 6 
6 7 
7 8 
8 9 
9 10 
3 2 217394434 
10 7 13022269 
6 7 283254485 
6 8 333042360 
4 6 442139372 
8 3 225045590 
10 4 922205209 
10 8 808296330 
9 2 486331361 
4 9 551176338 
1 8 5 
3 8 3 
3 8 4 
1 8 3 
4 8 1 
2 3 1 
2 3 1 
2 3 1 
2 4 1 
1 4 1

Sample Output

442139372 
333042360 
442139372 
283254485 
283254485 
217394434 
217394434 
217394434 
217394434 
217394434

Hint

图片 2

  作者从不权力号…因为时间限制大就改为权限题真的好啊?

  好像自个儿也讲不出什么不佳的道理。

  那题挺有意思的,对于自个儿这种码废+代码才干差+思维混乱+弱的人的话,这题挺杀时间的…五在那之中午就给它跪了。

  啊反正你正是能把四个行情产生二个依然七个矩形?把贰个水果产生一个点?

  然后就是总计覆盖点的权值第k小的矩形是哪些?

  hin有道理啊作者干吗正是想不出去呢吗吗?

  然后如若你会扫描线的话
正是一道全体二分裸题了?

  假设不会的话,你势必会差分是吧…

  把五个矩形按x轴拆成左加右减,树状数组搞搞。

  然后就是一道全部二分裸题了?

  为啥二个个 一百行都不要的
小编却打了如此多吧?

  思维混乱+弱啊!

  注意数组要开两倍,因为盘子恐怕有多少个。

#include   <iostream>
#include   <cstdio>
#include    <cstdlib>
#include    <algorithm>
#include    <vector>
#include    <cstring>
#include    <queue>
#include    <complex>
#include    <stack>
#define LL long long int
#define dob double
using namespace std;

const int N = 40010;
struct Node{int to,next;}E[N<<1];
struct Plate{
  int xl,xr,yl,yr,val;
  bool operator <(const Plate &p)const{
    return val<p.val;
  }
}plate[N<<1];
struct Fruit{
  int x,y,k,id;
  bool operator <(const Fruit &f)const{
    return x<f.x;
  }
}fruit[N],que1[N],que2[N];
struct Data{
  int x,l,r,val;
  bool operator <(const Data &l)const{
    return x<l.x;
  }
}Line[N<<1];
int n,P,Q,Ans[N],cnt,head[N],tot;
int fa[21][N],dep[N],size[N],son[N],top[N],dfn[N],tim,last[N];
struct Bit{
  int A[N],vis[N];
  inline int lb(int x){
    return x&-x;
  }
  inline void upd(int x,int val){
    for(;x<=n;x+=lb(x)){
      if(vis[x]!=tim)A[x]=0,vis[x]=tim;
      A[x]+=val;
    }
  }
  inline void update(int l,int r,int val){
    upd(l,val);upd(r+1,-val);
  }
  inline int query(int x,int ans=0){
    for(;x;x-=lb(x)){
      if(vis[x]!=tim)A[x]=0,vis[x]=tim;
      ans+=A[x];
    }
    return ans;
  }
}T;

inline int gi(){
  int x=0,res=1;char ch=getchar();
  while(ch>'9'||ch<'0'){if(ch=='-')res*=-1;ch=getchar();}
  while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
  return x*res;
}

inline void link(int u,int v){
  E[++tot]=(Node){v,head[u]};
  head[u]=tot;
}

inline void join(int x){
  for(int i=1;i<=15;++i)
    fa[i][x]=fa[i-1][fa[i-1][x]];
}

inline void dfs1(int x,int fat){
  fa[0][x]=fat;dep[x]=dep[fat]+1;join(x);
  dfn[x]=++tim;
  for(int e=head[x];e;e=E[e].next){
    int y=E[e].to;if(y==fa[0][x])continue;
    dfs1(y,x);
  }
  last[x]=tim;
}

inline int lca(int x,int y){
  if(x==y)return x;
  if(dep[x]<dep[y])swap(x,y);
  for(int i=15;i>=0;--i)
    if(dep[fa[i][x]]>=dep[y])
      x=fa[i][x];
  if(x==y)return x;
  for(int i=15;i>=0;--i)
    if(fa[i][x]!=fa[i][y])
      x=fa[i][x],y=fa[i][y];
  return fa[0][x];
}

inline int jump(int x,int goal){
  for(int i=15;i>=0;--i)
    if(dep[fa[i][x]]>dep[goal])
      x=fa[i][x];
  return x;
}

/*
  把答案(盘子)二分。
  把左边的矩形加进去。
  然后check,calc一下盘子个数,和k比较一下,划分一下。
  递归处理。
*/

inline void solve(int optl,int optr,int l,int r){
  if(optl>optr || l>r)return;++tim;
  if(optl==optr){
    for(int i=l;i<=r;++i)
      Ans[fruit[i].id]=plate[optl].val;
    return;
  }
  int mid=(optl+optr)>>1,tot1=0,tot2=0,k=l,tmp=0,cnt1=1,cnt2=l;
  for(int i=optl;i<=mid;++i){
    Line[++tmp]=(Data){plate[i].xl,plate[i].yl,plate[i].yr,1};
    Line[++tmp]=(Data){plate[i].xr+1,plate[i].yl,plate[i].yr,-1};
  }
  sort(Line+1,Line+tmp+1);
  while(cnt1<=tmp && cnt2<=r){
    if(Line[cnt1].x<=fruit[cnt2].x){
      int xl=Line[cnt1].l,xr=Line[cnt1].r,val=Line[cnt1].val;
      T.update(xl,xr,val);cnt1++;
    }
    else{
      int sum=T.query(fruit[cnt2].y);
      if(sum>=fruit[cnt2].k)que1[++tot1]=fruit[cnt2++];
      else fruit[cnt2].k-=sum,que2[++tot2]=fruit[cnt2++];
    }
  }
  while(cnt2<=r){
    int sum=T.query(fruit[cnt2].y);
    if(sum>=fruit[cnt2].k)que1[++tot1]=fruit[cnt2++];
    else fruit[cnt2].k-=sum,que2[++tot2]=fruit[cnt2++];
  }
  for(int i=1;i<=tot1;++i)fruit[k++]=que1[i];
  for(int i=1;i<=tot2;++i)fruit[k++]=que2[i];
  solve(optl,mid,l,l+tot1-1);
  solve(mid+1,optr,l+tot1,r);
}

int main()
{
  freopen("fruit_hnoi2015.in","r",stdin);
  freopen("fruit_hnoi2015.out","w",stdout);
  n=gi();P=gi();Q=gi();
  for(int i=1;i<n;++i){
    int u=gi(),v=gi();
    link(u,v);link(v,u);
  }
  dfs1(1,1);
  for(int i=1;i<=P;++i){
    int a=gi(),b=gi(),c=gi(),u=lca(a,b);
    if(dfn[a]>dfn[b])swap(a,b);
    if(u!=a)plate[++cnt]=(Plate){dfn[a],last[a],dfn[b],last[b],c};
    else{
      int v=jump(b,a);
      if(dfn[v]>1)plate[++cnt]=(Plate){1,dfn[v]-1,dfn[b],last[b],c};
      if(last[v]<n)plate[++cnt]=(Plate){dfn[b],last[b],last[v]+1,n,c};
    }
  }
  for(int i=1;i<=Q;++i){
    int a=gi(),b=gi(),k=gi();
    if(dfn[a]>dfn[b])swap(a,b);
    fruit[i]=(Fruit){dfn[a],dfn[b],k,i};
  }
  sort(plate+1,plate+cnt+1);
  sort(fruit+1,fruit+Q+1);
  solve(1,cnt,1,Q);
  for(int i=1;i<=Q;++i)printf("%d\n",Ans[i]);
  return 0;
}

  

 

相关文章