标签：return cst silver mat 材质 ems algo swap 接下来
第一行一个整数n，代表岛屿数量。
接下来n-1行，每行三个整数u,v,w，代表u号岛屿和v号岛屿由一条代价为c的桥梁直接相连，保证1<=u,v<=n且1<=c<=100000。
第n+1行，一个整数m，代表敌方机器能使用的次数。
接下来m行，每行一个整数ki，代表第i次后，有ki个岛屿资源丰富，接下来k个整数h1,h2,…hk，表示资源丰富岛屿的编号。
输出有m行，分别代表每次任务的最小代价。
对于100%的数据，2<=n<=250000,m>=1,sigma(ki)<=500000,1<=ki<=n-1
虚树
传说中的虚树，原来就是利用LCA和栈重建树啊……
码题过程十分曲折。先是一个a[i]打成i调半天，之后交上去又TTT，看了别的题解才注意到初始化邻接表不需要memset(hd)，只要标号置0，dp完一个点就删一个点的hd就可以了，交上去答案爆炸，发现没有开LL……开了LL程序突然爆了，发现是因为手写的min函数形参是int……删了min继续交，INF又开大了……
对自己的智商感到绝望.jpg
1 /*by SilverN*/ 2 #include<algorithm> 3 #include<iostream> 4 #include<cstring> 5 #include<cstdio> 6 #include<cmath> 7 #include<vector> 8 #define LL long long 9 using namespace std; 10 const int mxn=600010; 11 int read(){ 12 int x=0,f=1;char ch=getchar(); 13 while(ch<‘0‘ || ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();} 14 while(ch>=‘0‘ && ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} 15 return x*f; 16 } 17 struct edge{ 18 int v,nxt,w; 19 }e[mxn<<1],eg[mxn<<1]; 20 int hd1[mxn],hd2[mxn],mct=0,mct2=0; 21 void add_edge(int u,int v,int w){ 22 e[++mct].v=v;e[mct].nxt=hd1[u];e[mct].w=w;hd1[u]=mct;return; 23 } 24 void add(int u,int v){ 25 if(u==v)return; 26 eg[++mct2].v=v;eg[mct2].nxt=hd2[u];hd2[u]=mct2;return; 27 } 28 int fa[mxn][19],dep[mxn]; 29 LL dis[mxn]; 30 int in[mxn],out[mxn],dtime=0; 31 void DFS(int u,int ff){ 32 in[u]=++dtime; 33 dep[u]=dep[ff]+1; 34 for(int i=1;i<19;i++)fa[u][i]=fa[fa[u][i-1]][i-1]; 35 for(int i=hd1[u],v;i;i=e[i].nxt){ 36 if(e[i].v==ff)continue; 37 v=e[i].v; 38 fa[v][0]=u; 39 dis[v]=min(dis[u],(LL)e[i].w); 40 DFS(v,u); 41 } 42 out[u]=dtime; 43 return; 44 } 45 int cmp(int a,int b){return in[a]<in[b];}//按DFS序排序 46 int LCA(int x,int y){ 47 if(dep[x]<dep[y])swap(x,y); 48 for(int i=18;i>=0;i--) 49 if(dep[fa[x][i]]>=dep[y])x=fa[x][i]; 50 if(x==y)return y; 51 for(int i=18;i>=0;i--) 52 if(fa[x][i]!=fa[y][i])x=fa[x][i],y=fa[y][i]; 53 return fa[x][0]; 54 } 55 LL f[mxn]; 56 int DP(int u){ 57 f[u]=dis[u]; 58 LL res=0; 59 for(int i=hd2[u];i;i=eg[i].nxt){ 60 DP(eg[i].v); 61 res+=f[eg[i].v]; 62 } 63 if(res)f[u]=min(f[u],res); 64 hd2[u]=0; 65 return f[u]; 66 } 67 int a[mxn]; 68 int st[mxn],top=0; 69 void solve(){ 70 mct2=0;mct=0; 71 int i,j,K=read(); 72 for(i=1;i<=K;i++)a[i]=read(); 73 sort(a+1,a+K+1,cmp); 74 int cnt=0; 75 a[cnt=1]=a[1]; 76 for(i=2;i<=K;i++)if(in[a[i]]>out[a[cnt]])a[++cnt]=a[i]; 77 // for(i=2;i<=K;i++)if(LCA(a[cnt],a[i])!=a[cnt])a[++cnt]=a[i]; 78 st[top=1]=1; 79 for(i=1;i<=cnt;i++){ 80 int now=a[i]; 81 int tmp=LCA(a[i],st[top]); 82 while(1){ 83 if(dep[tmp]>=dep[st[top-1]]){ 84 add(tmp,st[top--]); 85 if(st[top]!=tmp)st[++top]=tmp; 86 break; 87 } 88 add(st[top-1],st[top]);--top; 89 } 90 if(st[top]!=now)st[++top]=now; 91 } 92 while(--top)add(st[top],st[top+1]); 93 DP(1); 94 printf("%lld\n",f[1]); 95 return; 96 } 97 int n,m; 98 int main(){ 99 int i,j,u,v,w; 100 n=read(); 101 for(i=1;i<n;i++){ 102 u=read();v=read();w=read(); 103 add_edge(u,v,w); 104 add_edge(v,u,w); 105 } 106 dis[1]=1e13; 107 DFS(1,0); 108 m=read(); 109 while(m--)solve(); 110 return 0; 111 }
标签：return cst silver mat 材质 ems algo swap 接下来
原文地址：http://www.cnblogs.com/SilverNebula/p/6628441.html
