運輸計劃

一.審題-關鍵點:

1.只能建造一個蟲洞,也就是隻能將一條邊的權值改成0

2.求的實質上是最長路最短

二.思路……

1.所有的都是一棵樹上的兩點的距離->自然想到lca求兩點距離

2.然後求完了距離要求最值。因爲是最長路最短,很自然想到二分答案

3.二分固定最長路最短,也就是說所有的路徑都要小於這個值x。(下面是check函數)

4.找出所有長度>x的路徑,存下來,開始去邊。

5.我們去掉的邊一定滿足兩個條件:

  1)所有>x的路徑都經過這條邊(不然肯定超啊。。。)

  2)超的最多的那條邊減掉這條邊的長度<=x

6.由條件1就能想到把所有路徑經過的道路標一標,經過次數==超過x的路徑數的存下來。這裏就能想到樹狀差分求一下就行啦。

7.最後肯定是這些存下來的邊中的權值最大的最優啊(這不是廢話嘛。。。)找到權值最大的,判斷一下條件2就做完咯。

代碼在這裏(敲黑板)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>

using namespace std;
const int maxn=300005;
vector<int> more;
int n,m,cnt; 
int v[maxn*2],w[maxn*2],nxt[maxn*2],head[maxn],len[maxn];
int d[maxn],dep[maxn],p[maxn][33],dis[maxn],anc[maxn];
int c[maxn];
struct node
{
	int x,y;
}a[maxn];
int cmp(int x,int y)
{
	return dis[x]>dis[y];
}
void add_edge(int x,int y,int z)
{
	cnt++;v[cnt]=y;w[cnt]=z;
	nxt[cnt]=head[x];head[x]=cnt;
}
int lca(int x,int y)
{
	if(d[x]<d[y]) swap(x,y);
	int i;for(i=0;(1<<i)<=d[x];i++);i--;
	for(int j=i;j>=0;j--) if(d[x]-(1<<j)>=d[y]) x=p[x][j];
	if(x==y) return x;
	for(int j=i;j>=0;j--)
	{
		if(p[x][j]!=p[y][j]) {x=p[x][j];y=p[y][j];}
	}
	return p[x][0];
}
void find(int x,int fa)
{
	d[x]=d[fa]+1;p[x][0]=fa;
	for(int i=head[x];i!=-1;i=nxt[i])
	{
		if(v[i]!=fa)
		{
			dep[v[i]]=dep[x]+w[i];
			len[v[i]]=w[i];
			find(v[i],x);
		}
	}
}
void prework()
{
	for(int j=1;j<=32;j++)
	{
		for(int i=1;i<=n;i++) p[i][j]=p[p[i][j-1]][j-1];
	}
} 
void dfs(int x)
{
	for(int i=head[x];i!=-1;i=nxt[i])
	{
		if(v[i]!=p[x][0]) 
		{
			dfs(v[i]);
			c[x]+=c[v[i]];	
		}
	}
}
bool ok(int x)
{
	more.clear();
	for(int i=1;i<=m;i++) {if(dis[i]>x) more.push_back(i);}
	if(more.empty()) return 1;
	memset(c,0,sizeof(c));
	sort(more.begin(),more.end(),cmp);
	for(int i=0;i<more.size();i++)
	{
		int tmp=lca(a[more[i]].x,a[more[i]].y);
		c[a[more[i]].x]++;c[a[more[i]].y]++;c[tmp]-=2;
	}
	dfs(1);
	int maxv=0;
	for(int i=1;i<=n;i++)
	{
		if(c[i]>=more.size())
		{
			maxv=max(maxv,len[i]);
		}
	}
	return dis[more[0]]-maxv<=x;
}
int main()
{
	memset(head,-1,sizeof(head));
	scanf("%d%d",&n,&m);
	for(int i=1;i<n;i++)
	{
		int x,y,z;scanf("%d%d%d",&x,&y,&z);
		add_edge(x,y,z);add_edge(y,x,z); 
	}
	find(1,0);prework();
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&a[i].x,&a[i].y);
		anc[i]=lca(a[i].x,a[i].y);
		dis[i]=dep[a[i].x]+dep[a[i].y]-2*dep[anc[i]];
	}
	int l=0,r=1e9,ans=0;
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(ok(mid))
		{
			ans=mid;r=mid-1;
		}
		else l=mid+1;
	}
	printf("%d",ans);
	return 0;
} 

 

 

 

 

 

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章