TYVJ P1577 泥泞的道路
题目链接:http://www.tyvj.cn/p/1577#
描述
输入格式
随后m行每行三个正整数x、y、w,用来描述一条道路,它连接x和y景点并且泥泞程度为w。
随后Q行,每行四个参数p1、p2、q1、q2,含义如下:
设数列Si表示你求得的第i个询问的结果(s0=1),则对于第i个询问:
X=(Si-1+p1)*p2 mod n + 1;
Y=(si-1+q1)*q2 mod n + 1;
输出格式
测试样例1
输入
4 4
1 2 2
2 3 1
1 3 3
3 4 5
2
3 1 0 2
1 1 2 1
输出
2
5
说明:第一个询问是{x=1,y=3},第二个询问是{x=4,y=1}。
备注
对于30%的数据,n<=1000,m<=3000,Q<=1000;
对于100%的数据,n<=100000,m<=300000,Q<=100000;
对于100%的数据,0<=p1、p2、q1、q2<n,w<300000。
在线问。两点之间的边的最大值最小。
先MST构造一棵树。
然后LCA询问,可是这道题爆栈,加栈也没用,还有这种LCA也第一次写,WA不少。
DFS变成BFS应该可以A了
1 #pragma comment(linker, "/STACK:1024000000,1024000000")
2 #include<iostream>
3 #include<algorithm>
4 #include<cmath>
5 #include<string.h>
6 #include<vector>
7 #include<stdio.h>
8
9
10 using namespace std;
11 #define N 623456
12 struct node
13 {
14 int x,y,z;
15 }a[N];
16
17 int f[N];
18
19 int find(int x)
20 {
21 if (f[x]!=x) f[x]=find(f[x]);
22 return f[x];
23 }
24
25 struct nod
26 {
27 int v,c,next;
28 }e[N<<];
29 int head[N],tot;
30
31 int cmp(node a,node b)
32 {
33 return a.z<b.z;
34 }
35
36 void add(int u,int v,int z)
37 {
38 e[tot].v=v;
39 e[tot].c=z;
40 e[tot].next=head[u];
41 head[u]=tot++;
42 }
43 int dp[N][],dep[N];
44 int dpmax[N][];
45
46 void dfs(int u,int fa,int t,int z)
47 {
48 dp[u][]=fa;
49 dpmax[u][]=z;
50 dep[u]=t;
51 for (int i=head[u];i!=-;i=e[i].next)
52 {
53 int v=e[i].v;
54 if (v==fa) continue;
55 //dpmax[v][0]=e[i].c;
56 dfs(v,u,t+,e[i].c);
57 }
58 }
59
60 int lca(int u,int v)
61 {
62 if (dep[u]<dep[v]) swap(u,v);
63 int dif=dep[u]-dep[v];
64
65 int res=;
66 for (int i=;i<;i++)
67 if ((dif>>i)&) res=max(res,dpmax[u][i]),u=dp[u][i];
68
69 if (u==v) return res;
70 for (int i=;i>=;i--)
71 if (dp[u][i]!=dp[v][i]) {
72 res=max(res,dpmax[u][i]),
73 res=max(res,dpmax[v][i]);
74 u=dp[u][i],v=dp[v][i];
75 }
76 if (u!=v) res=max(res,dpmax[u][]),res=max(res,dpmax[v][]);
77 return res;
78 }
79 inline int read()
80 {
81 int x=,f=;char ch=getchar();
82 while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
83 while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
84 return x*f;
85 }
86 int main()
87 {
88 int n,m;
89 // freopen("input.txt","r",stdin);
90 memset(head,-,sizeof(head));
91 tot=;
92 // scanf("%d%d",&n,&m);
93 n=read();
94 m=read();
95 for (int i=;i<=m;i++) //scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);
96 a[i].x=read(),a[i].y=read(),a[i].z=read();
97
98 sort(a+,a+m+,cmp);
99 for (int i=;i<=n;i++) f[i]=i;
for (int i=;i<=m;i++)
{
int x=find(a[i].x);
int y=find(a[i].y);
if (x!=y)
{
add(a[i].x,a[i].y,a[i].z);
add(a[i].y,a[i].x,a[i].z);
f[x]=y;
}
}
dfs(,,,);
//for (int i=1;i<=n;i++) cout<<dpmax[i][0]<<endl;
for (int i=;i<;i++)
for (int j=;j<=n;j++){
dpmax[j][i]=max(dpmax[j][i-],dpmax[dp[j][i-]][i-]);
dp[j][i]=dp[dp[j][i-]][i-];
}
int Q;
// scanf("%d",&Q);
Q=read();
int ans=;
while (Q--)
{
int p1,p2,q1,q2;
// scanf("%d%d%d%d",&p1,&p2,&q1,&q2);
p1=read(),p2=read(),q1=read(),q2=read();
int x=1ll*(ans+p1)*p2%n+;
int y=1ll*(ans+q1)*q2%n+;
// cout<<x<<y<<endl;
ans=lca(x,y);
printf("%d\n",ans);
}
return ;
}
TYVJ P1577 泥泞的道路的更多相关文章
- Codevs 1183 泥泞的道路
1183 泥泞的道路 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 传送门 题目描述 Description CS有n个小区,并且任意小区之间都有两条单向道路 ...
- 【Codevs1183】泥泞的道路
Position: http://codevs.cn/problem/1183/ List Codevs1183 泥泞的道路 List Description Input Output Sample ...
- codevs1183 泥泞的道路(01分数规划)
1183 泥泞的道路 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题目描述 Description CS有n个小区,并且任意小区之间都有两 ...
- codevs1183 泥泞的道路
题目描述 Description CS有n个小区,并且任意小区之间都有两条单向道路(a到b,b到a)相连.因为最近下了很多暴雨,很多道路都被淹了,不同的道路泥泞程度不同.小A经过对近期天气和地形的科学 ...
- 泥泞的道路(codevs 1183)
题目描述 Description CS有n个小区,并且任意小区之间都有两条单向道路(a到b,b到a)相连.因为最近下了很多暴雨,很多道路都被淹了,不同的道路泥泞程度不同.小A经过对近期天气和地形的科学 ...
- codevs 1183 泥泞的道路 01分数规划
题目链接 题目描述 Description CS有n个小区,并且任意小区之间都有两条单向道路(a到b,b到a)相连.因为最近下了很多暴雨,很多道路都被淹了,不同的道路泥泞程度不同.小A经过对近期天气和 ...
- CODEVS——T1183 泥泞的道路
时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题解 查看运行结果 题目描述 Description CS有n个小区,并且任意小区之间都有两条单向道路(a到 ...
- codevs1183泥泞的道路
题意:给定一张有向稠密图和通过每条边的时间和路程,问从1到n的路程/时间 最大为多少 正解:SPFA+二分答案 开始做的时候,想直接跑图论,后来发现好像不对(不然数据范围怎么这么小) 但是显然要用到图 ...
- codevs 1183 泥泞的道路 (二分+SPFA+差分约束)
/* 二分答案(注意精度) 对于每一个答案 有(s1+s2+s3...)/(t1+t2+t3...)>=ans 时符合条件 这时ans有变大的空间 对于上述不等式如果枚举每一条路显得太暴力 化简 ...
随机推荐
- syslog(),closelog()与openlog()--日志操作函数 (1)
文章出处:http://blog.csdn.net/xx77009833/archive/2010/07/30/5776383.aspx 为了满足某些目的,进行日志记录是很有必要的. 在典型的 LIN ...
- prototype中的ajax异步加载
jquery前台处理: var param = {a:a}; $.post("*.do",param,function(data) { var columHtml = " ...
- struts1标签库
Struts提供了五个标签库,即:HTML.Bean.Logic.Template和Nested. HTML标签 : 用来创建能够和Struts 框架和其他相应的HTML 标签交互的HTML 输入表单 ...
- Ubuntu创建应用快捷方式
Ubuntu创建应用快捷方式 新建一个.desktop文件 vi eclipse.desktop 然后又进行编辑 [Desktop Entry] Encoding=UTF-8 Name=eclipse ...
- mysql 备份解密脚本
#!/bin/bash #by sk 备份解码脚本 echo "-------------------------------------------------" functio ...
- inotify+rsync sersync+rsync实时同步服务
中小型网站搭建-数据实时的复制-inotify/sersync inotify是一种强大的,细粒度的.异步的文件系统事件监控机制(软件),linux内核从2.6.13起,加入inotify支持,通过i ...
- RN安卓原生模块
https://facebook.github.io/react-native/docs/native-modules-android.html RN实际就是依附在原生平台上,把各种各样的RN组件展示 ...
- 【css】【动画】【转发】旋转动画
<!DOCTYPE HTML> <html> <head> <meta charset="utf-8"> <s ...
- Linux中断底半部机制
参考: Linux下半部处理之软中断 linux中断底半部机制 <深入理解Linux内核>软中断/tasklet/工作队列 软中断和tasklet介绍 详解操作系统中断 Linux内核:中 ...
- PAT Basic 1044
1044 火星数字 火星人是以 13 进制计数的: 地球人的 0 被火星人称为 tret. 地球人数字 1 到 12 的火星文分别为:jan, feb, mar, apr, may, jun, jly ...