E. Qpwoeirut and Vertices(思维 + MST + lca + 线段树)

题目链接

我焯,这个是真的难…

参考题解

分析:

1.

首先,题目给的是一个图,但是通过样例我们发现并不需要那么多的边;
贪心的想,我们只需要最靠前的边使得所有节点连通就行了;
因此我们把边的编号当作边权,求一下MST(最小生成树);
那么我们的答案一定是在这个MST上找的;

2.

对于具体询问,我们转换一下问题:
原问题就转换成:给[l,r],求l~r之间所有点形成的连通块的最大边权;
注意这里是在MST上的,所以给一个[l,r],连通块都是确定的。(到这里就不会了…)

3.(注意这里我们都是在MST上搞的!)

这个有个知识盲点我不会…连续区间求lca一定是连通的;
先说结论吧:
如果我们能预处理出相邻两个点(i-1,i)路径上的边权的最大值mx[i];
那么每次查询[l,r],答案就是max(mx[i]);
因此我们可以用线段树维护mx[i]最值,每次询问[l,r]直接查一下就行;

4. 如何算出相邻两个点(i-1,i)路径上的边权的最大值mx[i]???

注意到任意两个点在树上的路径是唯一的,可以用lca来求!
因此我们可以在求lca的时候顺便求出两点路径上的边权mx。
这个就考到了对lca的理解了…不能光贴板子…
我们可以维护一个f_max[i][j],表示i这个点往上跳j步过程中的一个最大边权;
然后在bfs时 和 求lca时同时维护这个数组就行,这就比较简单了。

为什么维护相邻的点的lca这样做是对的?
首先就是任意两个点在树上的路径是确定的,唯一的;
其次连续区间内任意两点的lca一定包含在若干个[i,i+1]的lca内部,因为是一个连通图。

究极缝合怪,细节还是非常多的。

代码

#include<iostream>
#include<queue>
#include<cstring>
#include<vector>
#include<stdio.h>
#include<map>
#include<algorithm>
#include<deque>
#include<stack>
#include<set>
// #include <unordered_map>
#include<math.h>
#include<string.h>
#define IOS ios::sync_with_stdio(false),cin.tie(0);
using namespace std;
 
#define pb push_back
#define coutl cout<<"------------"<<endl;
#define fi first
#define se second
  
#define ire(x) scanf("%d",&x)
#define iire(a,b) scanf("%d %d",&a,&b)
#define lre(x) scanf("%lld",&x)
#define llre(a,b) scanf("%lld %lld",&a,&b)
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define endl "\n"
#define PI acos(-1.0)

typedef long long ll;
typedef unsigned long long ull;
      
typedef pair<int, int> PII;
typedef pair<double, int> PDI;
typedef pair<ll, ll> PLL;
typedef pair<double, double> PDD;
typedef pair<double, pair<int, double> > PDID;
typedef pair<char, char> PCC;
typedef pair<char, pair<int, int> > PCII;
typedef pair<int, pair<int, int> > PIII;
typedef pair<int, pair<int, pair<int, int> > > PIIII;
typedef pair<ll, pair<int, int> > PLII;
 
const int maxn = 2e5 + 7;
const int N = 1e6 + 7;
const int M = 1e6 + 7;
const int mod = 998244353;
const int inf = 0x3f3f3f3f;
// const ll INF = 0x3f3f3f3f3f3f3f3f;
const double pi = acos(-1);
const double eps = 1e-8;
  
ll gcd(ll a,ll b) {return b==0 ? a : gcd(b,a%b);}
ll lcm(ll a,ll b) {return a*b / gcd(a,b);}
ll qmi(ll a,ll b,ll p) {ll ans = 1; while(b) { if(b & 1) ans = ans * a % p; a = a * a % p; b >>= 1; } return ans;}
int lowbit(int x) {return x & (-x);}

int n,m,q;
int h[maxn],ne[maxn*2],e[maxn*2],w[maxn*2],idx;
int ff[maxn];	//并查集
int v[maxn];	//v[i]表示i到i-1路径上边权最大值

struct node	//线段树
{
	int l,r;
	int mx;
}tr[maxn*4];

//lca
int depth[maxn],f[maxn][22];    //2^15 > 4e4,所以15就可以跳完所有边
int f_max[maxn][22];	//每个点往上跳2^j过程中的最大边权

void add(int a,int b,int c)
{
    e[idx]=b; w[idx]=c; ne[idx]=h[a]; h[a]=idx++;
}

void bfs(int root)  //预处理depth[] 和 f[][]数组, 从根开始向下更新 
{
    for(int i=0;i<=n;i++) depth[i] = inf;
    depth[0] = 0;       //设置哨兵,处理跳出根节点了 
    depth[root] = 1;
	
    queue<int> q;
    q.push(root);

    while(q.size())
    {
        int no = q.front(); //父节点 
        q.pop();

        for(int i=h[no]; i!=-1; i=ne[i])    //子节点 
        {
            int j = e[i];

            //更新两个数组
            if(depth[j] > depth[no] + 1)    //还没更新过
            {
                depth[j] = depth[no] + 1;

                q.push(j);

                f[j][0] = no;   //子节点跳一步到父节点
                f_max[j][0] = w[i];
    
                for(int k=1; k<=19; k++)    //枚举往上跳2^k步 
                {
                	f[j][k] = f[f[j][k-1]][k-1];    //2^k = 2^(k-1) + 2^(k-1)
                	f_max[j][k] = max(f_max[j][k-1],f_max[f[j][k-1]][k-1]);
				}          
            }
        }   
    }
}

int lca(int a,int b)	//边求lca,边求最大值
{
	int ans = 0;
    //第一步,两个点先跳到同一层
    if(depth[a] < depth[b]) swap(a,b);  //让a跳到b那一层
    for(int k=19; k>=0; k--)    //从大到小枚举跳的步数
        if(depth[f[a][k]] >= depth[b])  //如果a跳完之后不超过b,则可以跳
        {
        	ans = max(ans,f_max[a][k]);	//注意先更新ans!
        	a = f[a][k];
		}

    if(a == b) return ans;    //跳到同一层后相等说明b就是公共祖先了

    //否则,a b同时往上跳
    for(int k=19; k>=0; k--)
        if(f[a][k] != f[b][k])  //跳完之后不相等一直跳 
        {
        	ans = max(ans,max(f_max[a][k],f_max[b][k]));
            a = f[a][k];
            b = f[b][k];    //跳到最近公共祖先的下一层为止 
        }
	
	//最后再往上跳一步,两个都要更新!
	ans = max(ans,max(f_max[a][0],f_max[b][0]));
	
    return ans; //在往上跳一步就是最近公共祖先了
}

void init()
{
	for(int i=0;i<=n;i++)
	{
		h[i] = -1;
		ff[i] = i;
		v[i] = 0;
		for(int j=0;j<22;j++) f[i][j] = f_max[i][j] = 0;
	}
	idx = 0;
}

void pushup(int u)
{
	tr[u].mx = max(tr[u<<1].mx, tr[u<<1|1].mx);
}

void build(int u,int l,int r)
{
	tr[u] = {l,r};
	
	if(l == r) 
	{
		tr[u].mx = v[l];
		return;
	}
	
	int mid = l + r >> 1;
	build(u<<1,l,mid);
	build(u<<1|1,mid+1,r);
	
	pushup(u);
}

int query(int u,int l,int r)
{
	if(l <= tr[u].l && r >= tr[u].r) return tr[u].mx;
	
	int mid = tr[u].l + tr[u].r >> 1;
	
	int ans = 0;
	if(l <= mid) ans = max(ans,query(u<<1, l, r));
	if(r > mid) ans = max(ans,query(u<<1|1, l, r));
	
	return ans;
}

int find(int x)
{
	return ff[x] == x ? x : ff[x]=find(ff[x]);
}

void solve()
{
	iire(n,m);
	ire(q);
	
	init();	//初始化
	
	//重构树,求MST
	for(int i=1;i<=m;i++)
	{
		int a,b;
		iire(a,b);
		
		int fa = find(a);
		int fb = find(b);
		
		if(fa != fb)
		{
			ff[fa] = fb;
			add(a,b,i);
			add(b,a,i);
		}
	}
	
	bfs(1);  //预处理depth[] 和 f[][]数组
	
	for(int i=2;i<=n;i++) v[i] = lca(i-1,i);	//边求lca,边求最大值
	
	build(1,1,n);	//建线段树,维护区间最大值
	
	while(q--)
	{
		int l,r;
		iire(l,r);
		
		if(l == r) cout<<"0 ";
		else cout<<query(1,l+1,r)<<' ';
	}
	cout<<endl;
}

int main()
{
	int t;
	ire(t);
	while(t--)
	{
		solve();
	}
	
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值