其他分享
首页 > 其他分享> > [DTOJ3996]:Lesson5!(DP+拓扑+线段树)

[DTOJ3996]:Lesson5!(DP+拓扑+线段树)

作者:互联网

题目描述

  “最短的捷径就是绕远路,绕远路就是我最短的捷径”
  转眼就$Stage\ X$了,$Stage\ X$的比赛路线可以看做一个$n$个点$m$条边的有向无环图,每条边长度都是$1$。杰洛$\cdot$齐贝林会选择走最长的那一条路径。
  迪亚哥$\cdot$布兰度决定摧毁一个城市以及所有关于该城市的边,由于变成恐龙后脑子有点问题,他想要让摧毁后的$Stage$最长路径最短,他想知道要摧毁哪个城市,及摧毁后最长路径的长度,如果有多个城市答案相同,则输出编号最小的那一个。


输入格式

  本题包含多组数据,输入第一行一个整数$T$代表数据组数
  每组数据第一行两个整数$n,m$表示点数,边数。
  每组数据第$2\sim m+1$行每行两个整数$x_i,y_i$表示有一条连接$x_i,y_i$的边。


输出格式

  对于每组数据,输出一行两个整数,表示删除的城市编号及删除该城市后最长路径的长度。


样例

样例输入:

1
6 5
1 3
1 4
3 6
3 4
4 5

样例输出:

1 2


数据范围与提示

对于所有数据,满足$T\leqslant 10,1\leqslant n\leqslant 100,000,0\leqslant m\leqslant 500,000$。


题解

$Dijkstra$不能跑最长路!!!

解释一下。

众所周知$Dijkstra$不能跑带负边权的最短路,而跑最长路也就相当于是跑带负边权的最短路,所以它死了……

那么回来考虑这道题。

对于$DAG$,首先想到拓扑。

不妨先跑正反拓扑计算出$dis_s$和$dis_t$分别表示正反拓扑的最长路,思想类似$DP$。

那么边$i$对答案的贡献就是$dis_{s_{i_u}}+dis{t_{i_v}}+1$,删除一个点就相当与删掉了与它相连的边。

快速修改用线段树就好啦。

时间复杂度:$\Theta(m\log n)$。

期望得分:$100$分。

实际得分:$100$分。


代码时刻

#include<bits/stdc++.h>
#define L(x) x<<1
#define R(x) x<<1|1
using namespace std;
struct rec{int nxt,to;}e[1000001];
int head[2][100001],cnt;
int n,m;
int a[100001];
int dis[2][100001],sum[100001],tr[400001],in[100001],ou[100001];
pair<int,int> ans;
queue<int> q;
void pre_work()
{
	memset(head,0,sizeof(head));
	memset(dis,0,sizeof(dis));
	memset(sum,0,sizeof(sum));
	memset(tr,0,sizeof(tr));
	memset(in,0,sizeof(in));
	memset(ou,0,sizeof(ou));
	memset(a,0,sizeof(a));
	cnt=0;ans=make_pair(1,n);
}
void add(bool id,int x,int y)
{
	e[++cnt].nxt=head[id][x];
	e[cnt].to=y;
	head[id][x]=cnt;
}
void pushup(int x){tr[x]=max(tr[L(x)],tr[R(x)]);}
void add(int x,int l,int r,int k)
{
	if(l==r)
	{
		sum[l]++;
		if(sum[l]>0)tr[x]=l;
		else{sum[l]=0;tr[x]=-1;}
		return;
	}
	int mid=(l+r)>>1;
	if(k<=mid)add(L(x),l,mid,k);
	else add(R(x),mid+1,r,k);
	pushup(x);
}
void del(int x,int l,int r,int k)
{
	if(l==r)
	{
		sum[l]--;
		if(sum[l]>0)tr[x]=l;
		else{sum[l]=0;tr[x]=-1;}
		return;
	}
	int mid=(l+r)>>1;
	if(k<=mid)del(L(x),l,mid,k);
	else del(R(x),mid+1,r,k);
	pushup(x);
}
int main()
{
	freopen("johnny.in","r",stdin);
	freopen("johnny.out","w",stdout);
	int T;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d",&n,&m);
		pre_work();
		for(int i=1;i<=m;i++)
		{
			int x,y;
			scanf("%d%d",&x,&y);
			add(0,x,y);in[y]++;
			add(1,y,x);ou[x]++;
		}
		for(int i=1;i<=n;i++)if(!in[i]){a[++a[0]]=i;q.push(i);}
		while(q.size())
		{
			int x=q.front();q.pop();
			for(int i=head[0][x];i;i=e[i].nxt)
			{
				if(dis[0][e[i].to]<dis[0][x]+1)dis[0][e[i].to]=dis[0][x]+1;
				in[e[i].to]--;
				if(!in[e[i].to]){a[++a[0]]=e[i].to;q.push(e[i].to);}
			}
		}
		for(int i=1;i<=n;i++)if(!ou[i])q.push(i);
		while(q.size())
		{
			int x=q.front();q.pop();
			for(int i=head[1][x];i;i=e[i].nxt)
			{
				if(dis[1][e[i].to]<dis[1][x]+1)dis[1][e[i].to]=dis[1][x]+1;
				ou[e[i].to]--;
				if(!ou[e[i].to])q.push(e[i].to);
			}
		}
		for(int i=1;i<=n;i++)add(1,0,n,dis[1][i]);
		for(int x=1;x<=n;x++)
		{
			for(int i=head[1][a[x]];i;i=e[i].nxt)del(1,0,n,dis[0][e[i].to]+dis[1][a[x]]+1);
			del(1,0,n,dis[1][a[x]]);
			if(tr[1]<ans.second||(ans.second==tr[1]&&a[x]<ans.first))ans=make_pair(a[x],tr[1]);
			for(int i=head[0][a[x]];i;i=e[i].nxt)add(1,0,n,dis[0][a[x]]+dis[1][e[i].to]+1);
			add(1,0,n,dis[0][a[x]]);
		}
		printf("%d %d\n",ans.first,ans.second);
	}
	return 0;
}

rp++

标签:Lesson5,int,sum,memset,tr,DP,sizeof,DTOJ3996,dis
来源: https://www.cnblogs.com/wzc521/p/11804336.html