P4899-[IOI2018]werewolf 狼人【Kruskal重构树,主席树】

正题

题目链接:https://www.luogu.com.cn/problem/P4899


题目大意

nnn个点的一张无向图,每次询问(s,t,l,r)(s,t,l,r)(s,t,l,r)表示询问能否找到一条s∼ts\sim tst的路径使得该路径可以分割成点的序号在[l,n][l,n][l,n][1,r][1,r][1,r]的两段。


解题思路

首先对原图构造两棵KruskalKruskalKruskal生成树,分别是最小序号最大的和最大序号最小的路径,然后我们每次询问是可以分别在树上跳到一个最高的位置。

然后问题就转换为询问两个树的子树交集,考虑使用主席树,对于第dfn1xdfn1_xdfn1x棵主席树上的第dfn2xdfn2_xdfn2x个位置加权值(dfn1xdfn1_xdfn1xdfn2xdfn2_xdfn2x分别表示xxx在第1/21/21/2棵树上的dfsdfsdfs序)。然后主席树询问即可。

时间复杂度O(nlog⁡n)O(n\log n)O(nlogn)


codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=6e5,T=20;
struct edge{int x,y,w;
}e[N];
struct node{int to,next;
};
struct Kul{node a[N<<1];int tot,cnt,num,fa[N],ls[N],w[N],dfn[N],rfn[N],ed[N],f[N][T+1];void addl(int x,int y){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;//printf("%d %d\n",x,y);return;}int find(int x){return (fa[x]==x)?x:(fa[x]=find(fa[x]));}void addp(int x,int y,int val){if(find(x)==find(y))return;w[++num]=val; addl(num,find(x));addl(num,find(y));fa[find(y)]=num;fa[find(x)]=num;return;}void dfs(int x){dfn[++cnt]=x;rfn[x]=cnt;for(int i=ls[x];i;i=a[i].next){int y=a[i].to;f[y][0]=x;dfs(y);}ed[x]=cnt;}void init(){for(int i=num;i>=1;i--)if(!f[i][0])dfs(i);w[0]=num+1;for(int j=1;j<=T;j++)for(int i=1;i<=num;i++)f[i][j]=f[f[i][j-1]][j-1];return;}int jump(int x,int val){for(int i=T;i>=0;i--)if(w[f[x][i]]<=val)x=f[x][i];return x;}
}k1,k2;
int n,m,q,cnt,rt[N],x[N],y[N],w[N<<5],ls[N<<5],rs[N<<5];
int Insert(int y,int l,int r,int pos){int x=++cnt;w[x]=w[y]+1;if(l==r)return x;int mid=(l+r)>>1;if(pos<=mid)ls[x]=Insert(ls[y],l,mid,pos),rs[x]=rs[y];else rs[x]=Insert(rs[y],mid+1,r,pos),ls[x]=ls[y];return x;
}
int Ask(int x,int y,int L,int R,int l,int r){if(!(w[y]-w[x]))return 0;if(l==L&&R==r)return w[y]-w[x];int mid=(L+R)>>1;if(r<=mid)return Ask(ls[x],ls[y],L,mid,l,r);if(l>mid)return Ask(rs[x],rs[y],mid+1,R,l,r);return Ask(ls[x],ls[y],L,mid,l,mid)+Ask(rs[x],rs[y],mid+1,R,mid+1,r);
}
bool cmp(edge x,edge y)
{return x.w<y.w;}
int main()
{scanf("%d%d%d",&n,&m,&q);k1.num=k2.num=n;for(int i=1;i<=n+m;i++)k1.fa[i]=k2.fa[i]=i;for(int i=1;i<=m;i++)scanf("%d%d",&e[i].x,&e[i].y),e[i].x++,e[i].y++,e[i].w=max(e[i].x,e[i].y);sort(e+1,e+1+m,cmp);for(int i=1;i<=m;i++)k1.addp(e[i].x,e[i].y,e[i].w),e[i].w=min(e[i].x,e[i].y);sort(e+1,e+1+m,cmp);for(int i=m;i>=1;i--)k2.addp(e[i].x,e[i].y,n-e[i].w+1);k1.init();k2.init();for(int i=1;i<=k1.num;i++)if(k1.dfn[i]<=n)rt[i]=Insert(rt[i-1],1,k2.num,k2.rfn[k1.dfn[i]]);else rt[i]=rt[i-1];while(q--){int l,r,s,t;scanf("%d%d%d%d",&s,&t,&l,&r);s++;t++;l++;r++;r=k1.jump(t,r);l=k2.jump(s,n-l+1);//printf("%d %d\n",l,r);if(Ask(rt[k1.rfn[r]-1],rt[k1.ed[r]],1,k2.num,k2.rfn[l],k2.ed[l])) printf("1\n");else printf("0\n");}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/320049.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

min_25筛

用途 设f(x)f(x)f(x)是一个积性函数&#xff0c;min_25筛可以在O(n34log⁡n)O(\frac{n^{\frac{3}{4}}}{\log n})O(lognn43​​)内求f(x)f(x)f(x)的前缀和&#xff1a;∑i1Nf(i)\sum_{i1}^{N}f(i)∑i1N​f(i) 使用要求&#xff1a;f(p),f(pk)f(p),f(p^k)f(p),f(pk)的值可以快速…

Service Fabric 用 Powershell 部署应用到本地

前置说明安装 Service Fabric SDK&#xff0c;会在本机 C:\Program Files\Microsoft SDKs\Service Fabric\Tools\PSModule\ServiceFabricSDK 生成部署脚本&#xff0c;如下图&#xff1a;用VS创建一个 Service Fabric 应用后&#xff0c;有一个部署脚本&#xff0c;位置在 [应用…

【期望】路径长度(金牌导航 期望-1)

路径长度 金牌导航 期望-1 题目大意 给出一个图&#xff0c;问你从1走到n的期望路径长度 输入样例 4 4 1 2 1 1 3 2 2 3 3 3 4 4输出样例 7.00数据范围 1⩽n⩽1051\leqslant n \leqslant 10^51⩽n⩽105 1⩽m⩽2n1\leqslant m\leqslant 2\times n1⩽m⩽2n 1⩽u,v⩽n1\…

牛客网 【每日一题】6月8日 [SCOI2005]最大子矩阵

链接&#xff1a; 文章目录题目描述题目描述 这里有一个n*m的矩阵&#xff0c;请你选出其中k个子矩阵&#xff0c;使得这个k个子矩阵分值之和最大。 注意&#xff1a;选出的k个子矩阵 不能相互重叠。 输入描述: 第一行为n,m,k&#xff08;1 ≤ n ≤ 100,1 ≤ m ≤ 2,1 ≤ k ≤ …

CF1054D-Changing Array【贪心】

正题 题目链接:https://www.luogu.com.cn/problem/CF1054D 题目大意 一个长度为nnn的序列&#xff0c;每个数小于2k2^k2k&#xff0c;可以选择一些数xorxorxor上2k−12^k-12k−1。要求使得满足alxoral1xor...xorar0a_l\ xor\ a_{l1}\ xor...xor\ a_r0al​ xor al1​ xor...xo…

ElasticSearch入门 附.Net Core例子

1.什么是ElasticSearch?Elasticsearch是基于Lucene的搜索引擎。它提供了一个分布式&#xff0c;支持多租户的全文搜索引擎&#xff0c;它具有HTTP Web界面和无模式JSON文档。 Elasticsearch是用Java开发的&#xff0c;根据Apache许可条款作为开源发布。----来自维基百科的解释…

[XSY]Tree Ext(矩阵树定理,拉格朗日插值,最小生成树,二分)

Tree Ext 这道题相当于把3道题合了起来。 要求修复的边中恰好有 k 条白边&#xff1a; 五颜六色的幻想乡(附拉格朗日插值法求多项式系数 ) bzoj2654 tree(WQS二分 新科技get) 是最小生成树计数而非生成树计数&#xff1a; BZOJ1016」[JSOI2008] 最小生成树计数 具体可以看看…

【期望】乘坐电梯(金牌导航 期望-2)

乘坐电梯 金牌导航 期望-2 题目大意 有n个人&#xff0c;对于没一个单位时间有p的概率最前面的1个人进电梯&#xff0c;有&#xff08;1-p&#xff09;的概率不进&#xff0c;问你t个单位时间后&#xff0c;电梯中的期望人数 样例输入 1 1 0.50 1 样例输出 1 0.5 样例输入…

动态规划练习【一】 背包问题

详细讲解 背包问题大汇总 文章目录背包问题大汇总01背包问题&#xff1a;思路&#xff1a;空间优化复杂度代码总结&#xff1a;完全背包问题&#xff1a;思路&#xff1a;代码&#xff1a;优化多重背包问题&#xff1a;思路&#xff1a;代码&#xff1a;单调队列优化混合三种背…

P3329-[ZJOI2011]最小割【最小割树】

正题 题目链接:https://www.luogu.com.cn/problem/P3329 题目大意 nnn个点mmm条边的无向图&#xff0c;每次询问一个xxx表示最小割不超过xxx的点对数量。 解题思路 我们对于两个点sss到ttt完成网络流后的残量网络上&#xff0c;与sss联通的点属于点集SSS&#xff0c;与ttt联通…

《.NET 性能优化》送书活动结果公布

截止到9月7日18&#xff1a;00&#xff08;规则本是12&#xff1a;00&#xff0c;忙的忘记了这事&#xff0c;18点截的图&#xff09;&#xff0c;本次送书活动《.NET 性能优化》共收到100多位同学参与回复&#xff0c;本次很多同学在看到活动的书 &#xff0c;自行就到异步社区…

【期望】期望收益(金牌导航 期望-3)

期望收益 金牌导航 期望-3 题目大意 给你一个01串&#xff0c;有些位置是未知的&#xff0c;连续的x个1贡献为想x2x^2x2&#xff0c;现在问你该串的期望贡献 输入样例 4 ????输出样例 4.1250数据范围 1⩽n⩽31051\leqslant n \leqslant 3\times 10^51⩽n⩽3105 解题…

杜教筛技巧随记

常见完全积性函数&#xff1a; ϵ(n)[n1]ϵ(n)[n1]ϵ(n)[n1]&#xff08;元函数&#xff0c;满足f∗ϵff*ϵff∗ϵf&#xff09; I(n)1I(n)1I(n)1 id(n)nid(n)nid(n)n 常见卷积&#xff1a; μ∗Iϵ\mu* Iϵμ∗Iϵ ϕ∗Iid\phi*Iidϕ∗Iid f(n)∑i1niϕ(i)f(n)\sum_{i1}^{n}i…

hdu-2844 Coins (混合背包+二进制优化)

HDU链接 文章目录题目描述&#xff1a;题意&#xff1a;题解&#xff08;代码&#xff09;题目描述&#xff1a; 输入描述: 输出描述: For each test case output the answer on a single line. 输入 3 10 1 2 4 2 1 1 2 5 1 4 2 1 0 0输出 8 4题意&#xff1a; 有n种硬币&…

CF1451F-Nullify The Matrix【结论题,博弈论】

正题 题目链接:https://www.luogu.com.cn/problem/CF1451F 题目大意 n∗mn*mn∗m的网格&#xff0c;每个网格上有数字&#xff0c;先后手轮流操作 每次操作选择一个只有往右和往下的路径&#xff0c;让第一个格子减去一个正整数并且随意修改后面的格子。要求完成后所有格子非…

.NET Core WebApi中实现多态数据绑定

什么是多态数据绑定&#xff1f;我们都知道在ASP.NET Core WebApi中数据绑定机制&#xff08;Data Binding&#xff09;负责绑定请求参数&#xff0c; 通常情况下大部分的数据绑定都能在默认的数据绑定器&#xff08;Binder&#xff09;中正常的进行&#xff0c;但是也会出现少…

【期望】期望分数(金牌导航 期望-4)

期望分数 金牌导航 期望-4 题目大意 告诉你一个01串中每个位置是1的概率&#xff0c;对于连续的x个1&#xff0c;贡献为x3x^3x3&#xff0c;问你期望贡献是多少 输入样例 3 0.5 0.5 0.5输入样例 6.0数据范围 1⩽N⩽1051\leqslant N \leqslant 10^51⩽N⩽105 解题思路…

[51NOD1847]奇怪的数学题(杜教筛+min_25筛+第二类斯特林数)

f(x)f(x)f(x)表示xxx的次大约数&#xff0c;有f(x)xx的最小质因数f(x)\frac{x}{x的最小质因数}f(x)x的最小质因数x​&#xff0c;那么 ∑i1n∑j1nsgcd(i,j)k∑i1n∑j1nf(gcd(i,j))k∑d1nf(d)k∑i1n∑j1n[gcd(i,j)d]∑d1nf(d)k∑i1⌊nd⌋∑j1⌊nd⌋[gcd(i,j)1]∑d1nf(d)k⋅(2∑i…

hdu 1059 Dividing

Hdu链接 文章目录题目描述题意&#xff1a;题解&#xff1a;代码&#xff1a;题目描述 输入描述: 输出描述: 示例1 输入 1 0 1 2 0 0 1 0 0 0 1 1 0 0 0 0 0 0输出 Collection #1: Cant be divided.Collection #2: Can be divided.题意&#xff1a; 有价值分别是1~6的6种…

UOJ#284-快乐游戏鸡【长链剖分,线段树】

正题 题目链接:https://uoj.ac/problem/284 题目大意 nnn个点的一棵树&#xff0c;每个点有一个wiw_iwi​表示至少死亡wiw_iwi​次才能通过这个点&#xff0c;否则就会死亡。只能往子节点走&#xff0c;mmm此询问从sis_isi​走到tit_iti​至少要死多少次。 解题思路 也就算我…