[BZOJ]2594 水管局长数据加强版(Wc2006)

  失踪人口回归。

  LCT一直是小C的弱项,特别是这种维护链的信息的,写挂了就会调代码调到心态爆炸。

  不过还好这一次的模板练习没有出现太多的意外。

 

Description

  SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的路径,接着通过信息化的控制中心通知路径上的水管进入准备送水状态,等到路径上每一条水管都准备好了,供水公司就可以开始送水了。嘟嘟一次只能处理一项送水任务,等到当前的送水任务完成了,才能处理下一项。
  在处理每项送水任务之前,路径上的水管都要进行一系列的准备操作,如清洗、消毒等等。嘟嘟在控制中心一声令下,这些水管的准备操作同时开始,但由于各条管道的长度、内径不同,进行准备操作需要的时间可能不同。供水公司总是希望嘟嘟能找到这样一条送水路径,路径上的所有管道全都准备就绪所需要的时间尽量短。嘟嘟希望你能帮助他完成这样的一个选择路径的系统,以满足供水公司的要求。另外,由于MY市的水管年代久远,一些水管会不时出现故障导致不能使用,你的程序必须考虑到这一点。
  不妨将MY市的水管网络看作一幅简单无向图(即没有自环或重边):水管是图中的边,水管的连接处为图中的结点。

Input

  输入文件第一行为3个整数:N, M, Q分别表示管道连接处(结点)的数目、目前水管(无向边)的数目,以及你的程序需要处理的任务数目(包括寻找一条满足要求的路径和接受某条水管坏掉的事实)。
以下M行,每行3个整数x, y和t,描述一条对应的水管。x和y表示水管两端结点的编号,t表示准备送水所需要的时间。我们不妨为结点从1至N编号,这样所有的x和y都在范围[1, N]内。
  以下Q行,每行描述一项任务。其中第一个整数为k:若k=1则后跟两个整数A和B,表示你需要为供水公司寻找一条满足要求的从A到B的水管路径;若k=2,则后跟两个整数x和y,表示直接连接x和y的水管宣布报废(保证合法,即在此之前直接连接x和y尚未报废的水管一定存在)。

Output

  按顺序对应输入文件中每一项k=1的任务,你需要输出一个数字和一个回车/换行符。该数字表示:你寻找到的水管路径中所有管道全都完成准备工作所需要的时间(当然要求最短)。

Sample Input

  4 4 3
  1 2 2
  2 3 3
  3 4 2
  1 4 2
  1 1 4
  2 1 4
  1 1 4

Sample Output

  2
  3

HINT

  N ≤ 100000
  M ≤ 1000000
  Q ≤ 100000
  任何时候我们考虑的水管网络都是连通的,即从任一结点A必有至少一条水管路径通往任一结点B。

 

Solution

  很容易看出题目要求我们维护无向图的最小生成树,查询两点树上路径最大权值,并支持删边操作。

  删边不好处理,而且题目没有要求我们在线,我们很容易想到倒过来做,把删边改为加边。

  于是用hash判断删去了哪些边(注意hash可能会撞),对删去所有边之后的图做最小生成树。

  这样就是LCT维护最小生成树的裸题了。

  对于每条新加入的边,查找两端点之间路径最大权值,如果小于那个权值,就删去那条边,并加入这条边。

 

#include <cstdio>
#include <cstring>
#include <algorithm>
#define l(a) son[a][0]
#define r(a) son[a][1]
#define nrt(a) (l(fa[a])==a || r(fa[a])==a)
#define MM 1000005
#define MN 100005
#define mod 10000007
using namespace std;
struct meg1{int x,y,z;}a[MM];
struct meg2{int g,x,y,z,pos;}b[MN];
struct edge{int nex,to,wt,pos;}e[MN<<1];
struct dmp{dmp *nex; int x,y,val;}*mp[mod];
int son[MN][2],fa[MN],tfa[MN],tson[MN],tg[MN],mx[MN];
int hr[MN],hua[MN],pin;
int n,m,p;

inline int read()
{
    int n=0,f=1; char c=getchar();
    while (c<'0' || c>'9') {if(c=='-')f=-1; c=getchar();}
    while (c>='0' && c<='9') {n=n*10+c-'0'; c=getchar();}
    return n*f;
}

void revs(int x) {tg[x]^=1; swap(l(x),r(x)); swap(tfa[x],tson[x]);}
void down(int x) {if (tg[x]) {revs(l(x)); revs(r(x)); tg[x]=0;}}
int L(int x) {down(x); return l(x)?L(l(x)):x;}
void aldown(int x) {if (nrt(x)) aldown(fa[x]); down(x);}
void upd(int &x,int y) {if (a[y].z>a[x].z) x=y;}
void update(int x) {mx[x]=0; upd(mx[x],tfa[x]); upd(mx[x],tson[x]); upd(mx[x],mx[l(x)]); upd(mx[x],mx[r(x)]);}
void rotate(int x)
{
    register int y,z,l,r;
    y=fa[x]; z=fa[y]; l=r(y)==x; r=l^1;
    if (nrt(y)) son[z][r(z)==y]=x;
    fa[y]=x; fa[x]=z; fa[son[x][r]]=y;
    son[y][l]=son[x][r]; son[x][r]=y;
    update(y);
}
void splay(int x)
{
    aldown(x);
    register int y,z;
    for (;nrt(x);rotate(x))
        {y=fa[x]; if (nrt(y)) z=fa[y],rotate(r(z)==y^r(y)==x?x:y);}
    update(x);
}
void access(int x)
{
    for (int i=0;x;x=fa[i=x])
    {
        splay(x); tson[x]=tfa[L(i)];
        r(x)=i; update(x);
    }
}

int getmx(int x,int y)
{
    access(x); splay(x); revs(x);
    access(y); splay(y); return mx[y];
}
void link(int x,int y,int z)
{
    access(x); splay(x);
    access(y); splay(y); revs(y);
    fa[y]=x; r(x)=y; tfa[y]=tson[x]=z;
    update(y); update(x);
}
void cut(int x,int y)
{
    access(x); splay(x); revs(x);
    access(y); splay(y);
    l(y)=fa[l(y)]=0; tson[x]=tfa[y]=0;
    update(x); update(y);
}

inline void ins(int x,int y,int z,int w)
{
    e[++pin]=(edge){hr[x],y,z,w}; hr[x]=pin;
    e[++pin]=(edge){hr[y],x,z,w}; hr[y]=pin;
}
void dfs(int x,int fat)
{
    fa[x]=fat;
    for (register int i=hr[x];i;i=e[i].nex)
        if (e[i].to!=fat) mx[e[i].to]=tfa[e[i].to]=e[i].pos,dfs(e[i].to,x);
}
inline bool gather(int x,int y)
{
    if (x==y) return false;
    hua[x]=y; return true;
}
int getf(int x) {return hua[x]?hua[x]=getf(hua[x]):x;}
bool cmp(const meg1& a,const meg1& b) {return a.z<b.z;}

int seamp(int h,int x,int y)
{
    for (dmp* q=mp[h];q;q=q->nex) if (q->x==x && q->y==y) return q->val;
    return 0;
}
void insmp(int h,int x,int y,int z)
{
    dmp* q=new(dmp);
    q->nex=mp[h]; q->x=x; q->y=y; q->val=z; mp[h]=q;
}
int geth(int x,int y) {return (1LL*x*300007+y)%mod;}

int main()
{
    register int i,lt;
    n=read(); m=read(); p=read();
    for (i=1;i<=m;++i)
    {
        a[i].x=read(); a[i].y=read(); a[i].z=read();
        if (a[i].x>a[i].y) swap(a[i].x,a[i].y);
    }
    for (i=1;i<=p;++i)
    {
        b[i].g=read(); b[i].x=read(); b[i].y=read();
        if (b[i].x>b[i].y) swap(b[i].x,b[i].y);
        if (b[i].g==2) insmp(geth(b[i].x,b[i].y),b[i].x,b[i].y,i);
    }
    sort(a+1,a+m+1,cmp);
    for (i=1;i<=m;++i)
    {
        if (lt=seamp(geth(a[i].x,a[i].y),a[i].x,a[i].y)) {b[lt].pos=i; continue;}
        if (gather(getf(a[i].x),getf(a[i].y))) ins(a[i].x,a[i].y,a[i].z,i);
    }
    dfs(1,0);
    for (i=p;i;--i)
    {
        if (b[i].g==1) b[i].z=a[getmx(b[i].x,b[i].y)].z;
        else if (b[i].g==2) {if (a[b[i].pos].z<a[lt=getmx(b[i].x,b[i].y)].z) cut(a[lt].x,a[lt].y),link(b[i].x,b[i].y,b[i].pos);}
    }
    for (i=1;i<=p;++i) if (b[i].g==1) printf("%d\n",b[i].z);
}

 

Last Word

  由于link和cut部分可能有点难以理解,LCT一直是令小C头疼的算法之一,但是小C还是成功把模板默出来了。

  然后迷之RE,并且发现树上部分节点的关系违背了伦理道德。调完hash之后,结果发现是在找最左边的节点的节点时忘记down操作了。

  果然LCT在维护信息时,考虑哪里需要update,哪里需要down,顺序是什么是最烦的。

  不得不说,D-map真的好用。

posted @ 2017-06-30 15:56  ACMLCZH  阅读(278)  评论(0编辑  收藏  举报