当且仅当盘子的路线是鲜果的路子的子路径(例如,个盘子就是顶点aiai到顶点bibi的不二法门(由于是树澳门金冠网站主页

4009: [Hnoi2015]接水果

Time Limit: 60 Sec  Memory Limit: 512 MB

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 的路径。 
菲菲每一回须要选取一个盘子去接当前的瓜果:一个行情能接住一个水果,当且仅当盘子的路径是鲜果的路径的子路径(例如
图1中从 3到7 的门径是从1到8的门路的子路径)。 
此地规定:从a 到b的门径与从b到
a的门路是一样条路子。当然为了增强难度,对于第 i
个水果,你要求采用能接住它的富有盘子中,权值第 kiki 小的那些盘子,每个盘子可重复使用(没有动用次数的上限:一个盘子接完一个水果后,前边还可继承接其余水果,只要它是鲜果路径的子路径)。幽香认为这些游乐很难,你能自在解决给他看呢? 
澳门金冠网站主页 1

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 的路子。 
清香每回需求选择一个行情去接当前的果品:一个盘子能接住一个水果,当且仅当盘子的路子是水果的路子的子路径(例如
图1中从 3到7 的途径是从1到8的途径的子路径)。 
那边规定:从a 到b的不二法门与从b到
a的途径是一模一样条途径。当然为了压实难度,对于第 i
个瓜果,你须求选用能接住它的装有盘子中,权值第 kiki 小的那几个盘子,每个盘子可重复使用(没有运用次数的上限:一个行情接完一个水果后,后边还可连续接其余水果,只要它是水果路径的子路径)。幽香认为这几个游戏很难,你能自在解决给她看吗? 
澳门金冠网站主页 2

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
小一定存在。

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

对于每个果子,输出一行表示选择的行情的权值。

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 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

Sample Output

442139372 
333042360 
442139372 
283254485 
283254485 
217394434 
217394434 
217394434 
217394434 
217394434

Hint

澳门金冠网站主页 3

  我尚未权力号…因为时限大就变成权限题真的好呢?

  好像我也讲不出什么倒霉的道理。

  那题挺有趣的,对于自己那种码废+代码能力差+思维混乱+弱的人的话,那题挺杀时间的…一个早晨就给它跪了。

  啊反正你就是能把一个盘子变成一个或者八个矩形?把一个水果变成一个点?

  然后就是统计覆盖点的权值第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;
}

  

 

Hint

澳门金冠网站主页 4

  我并未权限号…因为时限大就改成权限题真的好吧?

  好像自己也讲不出什么糟糕的道理。

  那题挺好玩的,对于自己那种码废+代码能力差+思维混乱+弱的人来说,那题挺杀时间的…一个中午就给它跪了。

  啊反正你就是能把一个盘子变成一个要么多个矩形?把一个水果变成一个点?

  然后就是统计覆盖点的权值第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;
}

  

 

相关文章