LCA+线段树 NOIP2016 天天爱跑步

天天爱跑步

题目描述

小c同学认为跑步非常有趣,于是决定制作一款叫做《天天爱跑步》的游戏。«天天爱跑步»是一个养成类游戏,需要玩家每天按时上线,完成打卡任务。

这个游戏的地图可以看作一一棵包含

现在有 $SiS_iSi,终点为 TiT_iTi 。每天打卡任务开始时,所有玩家在第000秒同时从自己的起点出发, 以每秒跑一条边的速度, 不间断地沿着最短路径向着自己的终点跑去, 跑到终点后该玩家就算完成了打卡任务。 (由于地图是一棵树, 所以每个人的路径是唯一的)$

小C想知道游戏的活跃度, 所以在每个结点上都放置了一个观察员。在结点 $jjj的观察员会选择在第WjW_jWj秒观察玩家, 一个玩家能被这个观察员观察到当且仅当该玩家在第WjW_jWj秒也理到达了结点 jjj 。小C想知道每个观察员会观察到多少人?$

注意: 我们认为一个玩家到达自己的终点后该玩家就会结束游戏, 他不能等待一段时间后再被观察员观察到。即对于把结点 $若他在第WjW_jWj秒重到达终点,则在结点jjj的观察员不能观察到该玩家;若他正好在第WjW_jWj秒到达终点,则在结点jjj的观察员可以观察到这个玩家。$

输入输出格式

输入格式：

第一行有两个整数

接下来

接下来一行 $nnn个整数,其中第jjj个整数为WjW_jWj , 表示结点jjj出现观察员的时间。$

接下来 $mmm行,每行两个整数SiS_iSi,和TiT_iTi,表示一个玩家的起点和终点。$

对于所有的数据,保证 $S_i,T_ileq n, 0leq W_jleq n1≤Si,Ti≤n,0≤Wj≤n 。$

输出格式：

输出1行

输入输出样例

输入样例#1：

6 3
2 3
1 2 
1 4 
4 5 
4 6 
0 2 5 1 2 3 
1 5 
1 3 
2 6

输出样例#1：

2 0 0 1 1 1

输入样例#2：

输出样例#2：

1 2 1 0 1

说明

【样例1说明】

对于1号点， $Wi=0W_i=0Wi=0，故只有起点为1号点的玩家才会被观察到，所以玩家1和玩家2被观察到，共有2人被观察到。$

对于2号点，没有玩家在第2秒时在此结点，共0人被观察到。

对于3号点，没有玩家在第5秒时在此结点，共0人被观察到。

对于4号点，玩家1被观察到，共1人被观察到。

对于5号点，玩家1被观察到，共1人被观察到。

对于6号点，玩家3被观察到，共1人被观察到。

【子任务】

每个测试点的数据规模及特点如下表所示。提示: 数据范围的个位上的数字可以帮助判断是哪一种数据类型。

【提示】

如果你的程序需要用到较大的栈空问 (这通常意味着需要较深层数的递归), 请务必仔细阅读选手日录下的文本当rumung:/stact.p″, 以了解在最终评测时栈空问的限制与在当前工作环境下调整栈空问限制的方法。

在最终评测时，调用栈占用的空间大小不会有单独的限制，但在我们的工作

环境中默认会有 8 MB 的限制。这可能会引起函数调用层数较多时，程序发生

栈溢出崩溃。

我们可以使用一些方法修改调用栈的大小限制。例如，在终端中输入下列命

令 ulimit -s 1048576

此命令的意义是，将调用栈的大小限制修改为 1 GB。

例如，在选手目录建立如下 sample.cpp 或 sample.pas

将上述源代码编译为可执行文件 sample 后，可以在终端中运行如下命令运

行该程序

./sample

如果在没有使用命令“ ulimit -s 1048576”的情况下运行该程序， sample

会因为栈溢出而崩溃；如果使用了上述命令后运行该程序，该程序则不会崩溃。

特别地，当你打开多个终端时，它们并不会共享该命令，你需要分别对它们

运行该命令。

请注意，调用栈占用的空间会计入总空间占用中，和程序其他部分占用的内

存共同受到内存限制。

此题我在csdn上写过了，还是比较详细，贴链接

这里宰贴一份别人的tarjan求LCA的代码。

 1 #include <iostream>
 2 #include <cstdlib>
 3 #include <cstring>
 4 #include <cstdio>
 5 #include <cmath>
 6 #include <algorithm>
 7 #include <ctime>
 8 #include <vector>
 9 #include <queue>
10 #include <map>
11 #include <set>
12 #include <string>
13 using namespace std;
14 typedef long long LL;
15 const int MAXN = 300011;
16 const int MAXM = 600011; 
17 int n,m,ecnt,first[MAXN],next[MAXM],to[MAXM],f[MAXN][20],deep[MAXN],ans[MAXN],val[MAXN],tong[MAXN],MAXD,w[MAXN],num[1000011];
18 int head[MAXN],tt[MAXM],nn[MAXM],father[MAXN],vis[MAXN];
19 vector<int>ljh[MAXN],ljh2[MAXN],ljh3[MAXN];
20 struct node{ int s,t,lca,len;}a[MAXN];
21 inline int getint(){
22     int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
23     if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
24 }
25 inline void link(int x,int y){ next[++ecnt]=first[x]; first[x]=ecnt; to[ecnt]=y; }
26 inline void LINK(int x,int y){ nn[++ecnt]=head[x]; head[x]=ecnt; tt[ecnt]=y; }
27 inline int find(int x){ if(father[x]!=x) father[x]=find(father[x]); return father[x]; }
28 inline void init(int x,int fa){ 
29     father[x]=x; vis[x]=1;
30     for(int i=head[x];i;i=nn[i]) {
31         int v=tt[i];
32         if(x==a[v].s&&vis[a[v].t]) a[v].lca=find(a[v].t);
33         if(x==a[v].t&&vis[a[v].s]) a[v].lca=find(a[v].s);
34     }
35     for(int i=first[x];i;i=next[i]) {
36         int v=to[i]; if(v==fa) continue; 
37         deep[v]=deep[x]+1; init(v,x); father[v]=x;
38         f[v][0]=x; 
39     }    
40 }
41  
42 inline int lca(int x,int y){
43     if(deep[x]<deep[y]) swap(x,y); int t=0; while((1<<t)<=deep[x]) t++; t--;
44     for(int i=t;i>=0;i--) if(deep[x]-(1<<i)>=deep[y]) x=f[x][i]; if(x==y) return y;
45     for(int i=t;i>=0;i--) if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i]; return f[x][0];
46 }
47  
48 inline void dfs(int x,int fa){
49     int now=w[x]+deep[x],cun; if(now<=MAXD) cun=tong[now];
50     for(int i=first[x];i;i=next[i]) {
51         int v=to[i]; if(v==fa) continue;
52         dfs(v,x);
53     }
54     tong[deep[x]]+=val[x]; if(now<=MAXD) ans[x]=tong[now]-cun;
55     for(int i=0,ss=ljh[x].size();i<ss;i++) tong[deep[ljh[x][i]]]--;
56 }
57  
58 inline void DFS(int x,int fa){
59     int now=deep[x]-w[x],cun; now+=300000; cun=num[now];
60     for(int i=first[x];i;i=next[i]) {
61         int v=to[i]; if(v==fa) continue;
62         DFS(v,x);
63     }
64     for(int i=0,ss=ljh2[x].size();i<ss;i++) num[300000+ljh2[x][i]]++;
65     ans[x]+=num[now]-cun;
66     for(int i=0,ss=ljh3[x].size();i<ss;i++) num[300000+ljh3[x][i]]--;
67 }
68  
69 inline void work(){  
70     n=getint(); m=getint(); int x,y; for(int i=1;i<n;i++) { x=getint(); y=getint(); link(x,y); link(y,x); }
71     for(int i=1;i<=n;i++) w[i]=getint(); ecnt=0;
72     for(int i=1;i<=m;i++) { a[i].s=getint(),a[i].t=getint(),val[a[i].s]++; LINK(a[i].s,i); LINK(a[i].t,i);}
73     deep[1]=1; init(1,0); for(int i=1;i<=n;i++) MAXD=max(MAXD,deep[i]);
74     for(int j=1;j<=19;j++) for(int i=1;i<=n;i++) f[i][j]=f[f[i][j-1]][j-1]; 
75     for(int i=1;i<=m;i++) {
76         a[i].len=deep[a[i].s]+deep[a[i].t]-deep[a[i].lca]*2;
77         ljh[a[i].lca].push_back(a[i].s);
78     }
79     dfs(1,0); 
80     for(int i=1;i<=m;i++) {
81         ljh2[a[i].t].push_back(deep[a[i].t]-a[i].len);
82         ljh3[a[i].lca].push_back(deep[a[i].t]-a[i].len);
83     }
84     DFS(1,0);
85     for(int i=1;i<=m;i++) if(deep[a[i].s]-deep[a[i].lca]==w[a[i].lca]) ans[a[i].lca]--;
86     for(int i=1;i<=n;i++) { printf("%d",ans[i]); if(i<n) printf(" "); }
87 }
88  
89 int main()
90 {
91     work();
92     return 0;
93 }