题意
求一条从1到N的路径,使路径上每条边A权值的最大值和B权值的最大值和最小。
刚开始写的暴力。
把边按A权值升序排序,每次操作加边,每次加边后用spfa维护图中的路径信息(可以不清空dis数组,这是个优化),这样可以暴力跑过去……还特别快……
#include <cstdio>
#include <string>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <queue>
#define N 50010
#define M 100010
using namespace std;
int n,m,cnt,Ans=1<<30;
int dis[N],V[N],G[N];
queue<int> Q;
struct edgef{
int a,b,x,y;
}e[M];
struct edge{
int t,nx,a,b;
}E[M<<2];
inline void reaD(int &x){
char Ch=getchar();x=0;
for(;Ch>'9'||Ch<'0';Ch=getchar());
for(;Ch>='0'&&Ch<='9';x=x*10+Ch-'0',Ch=getchar());
}
inline bool cmp(const edgef &A,const edgef &B){return A.a!=B.a?A.a<B.a:A.b<B.b;}
inline void InserT(int x,int y,int a,int b){
E[++cnt].t=y;E[cnt].a=a;E[cnt].b=b;E[cnt].nx=G[x];G[x]=cnt;
E[++cnt].t=x;E[cnt].a=a;E[cnt].b=b;E[cnt].nx=G[y];G[y]=cnt;
}
int main(){
reaD(n);reaD(m);
for(int i=1;i<=m;i++) reaD(e[i].x),reaD(e[i].y),reaD(e[i].a),reaD(e[i].b);
sort(e+1,e+1+m,cmp);
for(int i=2;i<=n;i++) dis[i]=1<<30;
for(int i=1,A;i<=m;i++){
A=e[i].a;
InserT(e[i].x,e[i].y,e[i].a,e[i].b);
Q.push(e[i].x);Q.push(e[i].y);
V[e[i].x]=V[e[i].y]=1;
while(!Q.empty()){
int x=Q.front();Q.pop();V[x]=0;
for(int i=G[x];i;i=E[i].nx)
if(dis[E[i].t]>max(dis[x],E[i].b)){
dis[E[i].t]=max(dis[x],E[i].b);
if(!V[E[i].t]) Q.push(E[i].t),V[E[i].t]=1;
}
}
Ans=min(Ans,A+dis[n]);
}
if(Ans==1<<30) puts("-1"); else printf("%d\n",Ans);
return 0;
}
当然UOJ没有那么友好……所以改成lct来维护就好了
#include <cstdio>
#include <algorithm>
#define N 200010
using namespace std;
int n,m,Ans=1<<30;
int fa[N],f[N],ch[N][2],v[N],maxn[N],rev[N],sta[N],tp;
struct edgef{
int a,b,x,y;
friend bool operator <(edgef A,edgef B){
return A.a<B.a;
}
}e[N];
inline char C(){
static char buf[100000],*p1=buf,*p2=buf;
if(p1==p2){
p2=(p1=buf)+fread(buf,1,100000,stdin);
if(p1==p2) return EOF;
}
return *p1++;
}
inline void reaD(int &x){
char Ch=C();x=0;
for(;Ch>'9'||Ch<'0';Ch=C());
for(;Ch>='0'&&Ch<='9';x=x*10+Ch-'0',Ch=C());
}
int fifa(int x){return fa[x]==x?x:fa[x]=fifa(fa[x]);}
inline int isl(int x){return ch[f[x]][1]==x;}
inline int isr(int x){return ch[f[x]][0]!=x&&ch[f[x]][1]!=x;}
inline void upd(int x){
if(!x)return;maxn[x]=x;
if(ch[x][0]&&v[maxn[ch[x][0]]]>v[maxn[x]]) maxn[x]=maxn[ch[x][0]];
if(ch[x][1]&&v[maxn[ch[x][1]]]>v[maxn[x]]) maxn[x]=maxn[ch[x][1]];
}
inline void pushdown(int x){
if(!(x&&rev[x])) return ;
swap(ch[x][0],ch[x][1]);
if(ch[x][0]) rev[ch[x][0]]^=1;
if(ch[x][1]) rev[ch[x][1]]^=1;
rev[x]=0;
}
inline void rot(int x){
int y=f[x],z=f[y],lor=isl(x);
if(!isr(y)) ch[z][ch[z][1]==y]=x;
f[x]=z;ch[y][lor]=ch[x][lor^1];
if(ch[y][lor]) f[ch[y][lor]]=y;
ch[x][lor^1]=y;f[y]=x;
upd(y);upd(x);
}
inline void splay(int x){
sta[++tp]=x;
for(int i=x;!isr(i);i=f[i]) sta[++tp]=f[i];
while(tp) pushdown(sta[tp--]);
for(;!isr(x);rot(x))if(!isr(f[x]))
if(isl(f[x])^isl(x)) rot(x); else rot(f[x]);
}
inline void access(int x){
int t=0;
for(;x;t=x,x=f[x])
splay(x),ch[x][1]=t,upd(x);
}
inline void reverse(int x){access(x);splay(x);rev[x]^=1;}
inline int newd(int w,int u){v[u]=w;return u;}
inline void link(int x,int y){reverse(x);f[x]=y;splay(x);}
inline void cut(int x,int y){reverse(x);access(y);splay(y);ch[y][0]=f[x]=0;}
inline int query(int x,int y){reverse(x);access(y);splay(y);return maxn[y];}
int main(){
freopen("forest.in","r",stdin);
freopen("forest.out","w",stdout);
reaD(n);reaD(m);
for(int i=1;i<=m;i++)reaD(e[i].x),reaD(e[i].y),reaD(e[i].a),reaD(e[i].b);
sort(e+1,e+1+m);
for(int i=1;i<=n;i++) fa[i]=i;
for(int i=1;i<=m;i++){
int x=fifa(e[i].x),y=fifa(e[i].y),z;
if(x!=y){
link(e[i].x,newd(e[i].b,n+i));
link(e[i].y,n+i);
fa[x]=y;
}
else{
int k=query(e[i].x,e[i].y);
if(v[k]>e[i].b){
cut(k,e[k-n].x);
cut(k,e[k-n].y);
link(e[i].x,newd(e[i].b,n+i));
link(e[i].y,n+i);
}
}
if(fifa(1)==fifa(n)) Ans=min(Ans,e[i].a+v[query(1,n)]);
}
if(Ans==(1<<30)) puts("-1"); else printf("%d\n",Ans);
return 0;
}