国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 學院 > 開發設計 > 正文

[POJ2114]Boatherds(點分治+二分)

2019-11-08 19:44:14
字體:
來源:轉載
供稿:網友

題目描述

傳送門 題意:一棵樹,每一次詢問樹上是否存在兩個點距離為x

題解

裸的點分治 暴力就行

代碼

#include<algorithm>#include<iostream>#include<cstring>#include<cstdio>#include<cmath>using namespace std;#define N 10005int n,x,y,m,sum,root;int q[N],ans[N];int tot,point[N],nxt[N*2],v[N*2],c[N*2];int size[N],big[N],d[N],deep[N];bool vis[N];void clear(){ x=y=m=sum=root=0; memset(ans,0,sizeof(ans)); tot=0;memset(point,0,sizeof(point)); memset(vis,0,sizeof(vis));}void add(int x,int y,int z){ ++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y; c[tot]=z;}void getroot(int x,int fa){ size[x]=1;big[x]=0; for (int i=point[x];i;i=nxt[i]) if (v[i]!=fa&&!vis[v[i]]) { getroot(v[i],x); size[x]+=size[v[i]]; big[x]=max(big[x],size[v[i]]); } big[x]=max(big[x],sum-size[x]); if (big[x]<big[root]) root=x;}void getdeep(int x,int fa){ deep[++deep[0]]=d[x]; for (int i=point[x];i;i=nxt[i]) if (v[i]!=fa&&!vis[v[i]]) { d[v[i]]=d[x]+c[i]; getdeep(v[i],x); }}int findl(int l,int r,int k){ int ans=0; while (l<=r) { int mid=(l+r)>>1; if (deep[mid]==k) ans=mid,r=mid-1; else if (deep[mid]<k) l=mid+1; else r=mid-1; } return ans;}int findr(int l,int r,int k){ int ans=-1; while (l<=r) { int mid=(l+r)>>1; if (deep[mid]==k) ans=mid,l=mid+1; else if (deep[mid]<k) l=mid+1; else r=mid-1; } return ans;}int calc(int x,int now,int k){ d[x]=now;deep[0]=0; getdeep(x,0); sort(deep+1,deep+deep[0]+1); int t=0; for (int i=1;i<=deep[0];++i) { if (deep[i]+deep[i]>k) break; int l=findl(i+1,deep[0],k-deep[i]); int r=findr(i+1,deep[0],k-deep[i]); t+=r-l+1; } return t;}void dfs(int x){ for (int i=1;i<=m;++i) ans[i]+=calc(x,0,q[i]); vis[x]=1; for (int i=point[x];i;i=nxt[i]) if (!vis[v[i]]) { for (int j=1;j<=m;++j) ans[j]-=calc(v[i],c[i],q[j]); sum=size[v[i]];root=0; getroot(v[i],0); dfs(root); }}int main(){ while (~scanf("%d",&n)) { if (!n) break; clear(); for (int i=1;i<=n;++i) { while (~scanf("%d",&x)) { if (!x) break; scanf("%d",&y); add(i,x,y);add(x,i,y); } } while (~scanf("%d",&x)) { if (!x) break; q[++m]=x; } big[0]=N; sum=n;root=0; getroot(1,0); dfs(root); for (int i=1;i<=m;++i) if (ans[i]) puts("AYE"); else puts("NAY"); puts("."); }}
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 英德市| 新宾| 霸州市| 民权县| 孟津县| 罗山县| 从江县| 湾仔区| 泗洪县| 来宾市| 囊谦县| 密云县| 邢台市| 阜平县| 宁晋县| 栾城县| 湟源县| 南木林县| 鄂伦春自治旗| 湟中县| 伊金霍洛旗| 卓尼县| 贵阳市| 定襄县| 邳州市| 大邑县| 石林| 台前县| 保康县| 潼南县| 彭山县| 梓潼县| 新和县| 万山特区| 泰安市| 平阳县| 迁安市| 聂拉木县| 陈巴尔虎旗| 贵阳市| 始兴县|