[luogu P1967][NOIp2013] 货车运输

题目描述

A 国有 n 座城市，编号从 1 到 n，城市之间有 m 条双向道路。每一条道路对车辆都有重量限制，简称限重。现在有 q 辆货车在运输货物，司机们想知道每辆车在不超过车辆限重的情况下，最多能运多重的货物。

输入输出格式

输入格式：

输入文件名为 truck.in。

输入文件第一行有两个用一个空格隔开的整数 n，m，表示 A 国有 n 座城市和 m 条道

路。接下来 m 行每行 3 个整数 x、 y、 z，每两个整数之间用一个空格隔开，表示从 x 号城市到 y 号城市有一条限重为 z 的道路。注意： x 不等于 y，两座城市之间可能有多条道路。

接下来一行有一个整数 q，表示有 q 辆货车需要运货。

接下来 q 行，每行两个整数 x、y，之间用一个空格隔开，表示一辆货车需要从 x 城市运输货物到 y 城市，注意： x 不等于 y 。

输出格式：

输出文件名为 truck.out。

输出共有 q 行，每行一个整数，表示对于每一辆货车，它的最大载重是多少。如果货

车不能到达目的地，输出-1。

输入输出样例

输入样例#1：

输出样例#1：

说明

对于 30%的数据，0 < n < 1,000，0 < m < 10,000，0 < q< 1,000；

对于 60%的数据，0 < n < 1,000，0 < m < 50,000，0 < q< 1,000；

对于 100%的数据，0 < n < 10,000，0 < m < 50,000，0 < q< 30,000，0 ≤ z ≤ 100,000。

因为这题在bzoj里又是那啥啥，所以就又交luogu了

代码这么长又比别人慢真是太尴尬了

这么慢思路肯定很easy对吧

->先求用kruskal最大生成树，可忽略重边问题

->以重心为根造树

->对于每一组询问(ss,tt)，求其lca

->我这边模拟了它从ss和tt走到lca，程序秒速变慢 QwQ(懒癌晚期)

 #include<cstdio>

 #include<cstring>

 #include<iostream>

 #include<vector>

 #include<algorithm>

 #define inf 0x3f3f3f3f

 using namespace std;

 #define ss first

 #define tt second

 inline int read(){

     char ch;

     int re=;

     bool flag=;

     while((ch=getchar())!='-'&&(ch<''||ch>''));

     ch=='-'?flag=:re=ch-'';

     while((ch=getchar())>=''&&ch<='')  re=re*+ch-'';

     return flag?-re:re;

 }

 struct odge{

     int from,to,w;

     odge(int from=,int to=,int w=):

         from(from),to(to),w(w){}

     bool operator < (const odge &n1) const {

         return w>n1.w;

     }

 };

 struct edge{

     int to,w,next;

     edge(int to=,int w=,int next=):

         to(to),w(w),next(next){}

 };

 const int maxn=,maxm=;

 const int maxlogn=;

 int n,m;

 vector<odge> odges;

 vector<edge> edges;

 int head[maxn];

 int par[maxn];

 int cnt=;

 int F[maxn],son[maxn];

 bool vis[maxn];

 int sum,root;

 int depth[maxn],parent[maxlogn][maxn];

 int q;

 typedef pair<int,int> PII;

 vector<PII> ques;

 int length[maxn];

 inline void add_edge(odge e){

     edges.push_back(edge(e.to,e.w,head[e.from]));

     head[e.from]=++cnt;

     edges.push_back(edge(e.from,e.w,head[e.to]));

     head[e.to]=++cnt;

 }

 int find(int x){

     return par[x]==x?x:par[x]=find(par[x]);

 }

 void kruskal(){

     memset(head,,sizeof head);

     edges.push_back(edge(,,));

     int cnt1=,cnt2=;

     while(cnt2<m){

         odge e=odges[cnt2++];

         int dx=find(e.from),dy=find(e.to);

         if(dx==dy)  continue;

         add_edge(e);

         par[dy]=dx;

         if(++cnt1==n)  break;

     }

 }

 void getroot(int x,int fa){

     son[x]=;  F[x]=;

     for(int ee=head[x];ee;ee=edges[ee].next)

         if(edges[ee].to!=fa){

             getroot(edges[ee].to,x);

             son[x]+=son[edges[ee].to];

             F[x]=max(F[x],son[edges[ee].to]);

         }

     F[x]=max(F[x],sum-son[x]);

     if(F[x]<F[root])  root=x;

 }

 void dfs(int x,int fa,int deep){

     parent[][x]=fa;

     depth[x]=deep;

     for(int ee=head[x];ee;ee=edges[ee].next)

         if(edges[ee].to!=fa){

             length[edges[ee].to]=edges[ee].w;

             dfs(edges[ee].to,x,deep+);

         }

 }

 void init(){

     n=read();  m=read();

     for(int i=;i<m;i++){

         int from=read(),to=read(),w=read();

         odges.push_back(odge(from,to,w));

     }

     for(int i=;i<=n;i++)

         par[i]=i;

     sort(odges.begin(),odges.end());

     kruskal();

     sum=F[root=]=n;

     getroot(,);

     dfs(root,-,);

     for(int k=;k+<maxlogn;k++)

         for(int v=;v<=n;v++)

             if(parent[k][v]<)  parent[k+][v]=-;

             else  parent[k+][v]=parent[k][parent[k][v]];

     q=read();

     for(int i=;i<q;i++){

         int from=read(),to=read();

         ques.push_back(make_pair(from,to));

     }

 }

 int getlca(int u,int v){

     if(depth[u]>depth[v])  swap(u,v);

     for(int k=;k<maxlogn;k++)

         if(depth[v]-depth[u]>>k&)  v=parent[k][v];

     if(u==v)  return u;

     for(int k=maxlogn-;k>=;k--)

         if(parent[k][u]!=parent[k][v]){

             u=parent[k][u];

             v=parent[k][v];

         }

     return parent[][u];

 }

 void solve(){

     for(int i=;i<q;i++){

         int u=ques[i].ss,v=ques[i].tt;

         if(find(u)!=find(v)){

             puts("-1");

             continue;

         }

         int ans=inf;

         int lca=getlca(ques[i].ss,ques[i].tt);

         while(u!=lca){

             ans=min(ans,length[u]);

             u=parent[][u];

         }

         while(v!=lca){

             ans=min(ans,length[v]);

             v=parent[][v];

         }

         printf("%d\n",ans);

     }

 }

 int main(){

     //freopen("temp.in","r",stdin);

     init();

     solve();

     return ;

 }

天亮的时候你还没回来

亲爱的我们的孩子醒来了