BZOJ卡不过灰常蛋疼(毕竟人蠢自带巨大常数

这和动态维护最小生成树很像,但加边变成了删边,似乎没法做了。

然后根据之前的套路离线做,删边变成加边,就可以做了orz

二分查找的:(慢

// It is made by XZZ
#include<cstdio>
#include<algorithm>
#include<map>
#include<stack>
#define il inline
#define rg register
#define vd void
#define sta static
using namespace std;
il int gi(){
rg int x=0,f=1;rg char ch=getchar();
while(ch<'0'||ch>'9')f=ch=='-'?-1:f,ch=getchar();
while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
return x*f;
}
const int maxn=100001;
stack<int>stk;
namespace LCT{
typedef const int& fast;
int ch[maxn*2][2],fa[maxn*2],mx[maxn*2],w[maxn*2],index,A[maxn*2],B[maxn*2];bool rev[maxn*2];
il vd upd(fast x){
mx[x]=x;
if(w[mx[ch[x][0]]]>w[mx[x]])mx[x]=mx[ch[x][0]];
if(w[mx[ch[x][1]]]>w[mx[x]])mx[x]=mx[ch[x][1]];
}
il vd Rev(fast x){if(x)rev[x]^=1,swap(ch[x][0],ch[x][1]);}
il vd down(fast x){if(rev[x])Rev(ch[x][0]),Rev(ch[x][1]),rev[x]=0;}
il bool isrt(fast x){return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;}
il vd rotate(fast x){
sta int y,z,o;y=fa[x],z=fa[y],o=x==ch[y][1];
if(!isrt(y))ch[z][y==ch[z][1]]=x;fa[x]=z;
ch[y][o]=ch[x][!o],fa[ch[x][!o]]=y;
ch[x][!o]=y,fa[y]=x;
upd(y);
}
il vd splay(fast x){
sta int y,z,stk[maxn],top;stk[top=1]=x;
for(rg int i=x;!isrt(i);i=fa[i])stk[++top]=fa[i];
while(top)down(stk[top--]);
for(y=fa[x],z=fa[y];!isrt(x);rotate(x),y=fa[x],z=fa[y])
if(!isrt(y))rotate(((x==ch[y][0])^(y==ch[z][0]))?y:x);
upd(x);
}
il vd access(int x){for(rg int y=0;x;x=fa[y=x])splay(x),ch[x][1]=y,upd(x);}
il vd makert(fast x){access(x),splay(x),Rev(x);}
il vd link(fast x,fast y){makert(x),fa[x]=y;}
il vd cut(fast x,fast y){makert(x),access(y),splay(y);fa[x]=ch[y][0]=0;}
il int Query(fast x,fast y){makert(x),access(y),splay(y);return mx[y];}
il int find(int x){access(x),splay(x);while(ch[x][0])x=ch[x][0];return x;}
il vd Link(fast x,fast y,fast _w){
if(find(x)^find(y)){w[++index]=_w,link(A[index]=x,index),link(B[index]=y,index);return;}
sta int z;z=Query(x,y);
if(_w>=w[z])return;
cut(A[z],z),cut(B[z],z);
w[z]=_w,link(A[z]=x,z),link(B[z]=y,z);
}
};
namespace GGMap{
struct gzy{int x,y,*k;}s[maxn];
il bool cmp(const gzy&a,const gzy&b){return a.x==b.x?a.y<b.y:a.x<b.x;}
int tot=0;
il vd insert(int a,int b,int*c){++tot,s[tot].x=a,s[tot].y=b,s[tot].k=c;}
il vd prepare(){sort(s+1,s+tot+1,cmp);}
il int* find(int a,int b){
rg int l=1,r=tot,mid;
while(l<r){
mid=(l+r)>>1;
if(s[mid].x<a||(s[mid].x==a&&s[mid].y<b))l=mid+1;
else r=mid;
}
if(s[l].x==a&&s[l].y==b)return s[l].k;
else return NULL;
}
}
int opt[maxn],A[maxn],B[maxn],C[maxn];
int X[10000001],Y[10000001],Z[10000001];
int main(){
freopen("tube_strong.in","r",stdin);
freopen("tube_strong.out","w",stdout);
int n=gi(),m=gi(),q=gi();
LCT::index=n;
for(rg int i=1;i<=m;++i){X[i]=gi(),Y[i]=gi(),Z[i]=gi();if(X[i]>Y[i])swap(X[i],Y[i]);}
for(rg int i=1;i<=q;++i){
opt[i]=gi(),A[i]=gi(),B[i]=gi();if(A[i]>B[i])swap(A[i],B[i]);
if(opt[i]==2)GGMap::insert(A[i],B[i],&C[i]);
}
GGMap::prepare();
for(rg int i=1;i<=m;++i){
int *it=GGMap::find(X[i],Y[i]);
if(it==NULL)LCT::Link(X[i],Y[i],Z[i]);
else *it=Z[i];
}
while(q){
if(opt[q]==1)stk.push(LCT::w[LCT::Query(A[q],B[q])]);
else LCT::Link(A[q],B[q],C[q]);
--q;
}
while(!stk.empty())printf("%d\n",stk.top()),stk.pop();
return 0;
}

map+一丁点哈希优化:(较上面快,但还是慢

// It is made by XZZ
#include<cstdio>
#include<algorithm>
#include<map>
#include<stack>
#define il inline
#define rg register
#define vd void
#define sta static
using namespace std;
il int gi(){
rg int x=0,f=1;rg char ch=getchar();
while(ch<'0'||ch>'9')f=ch=='-'?-1:f,ch=getchar();
while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
return x*f;
}
const int maxn=100001;
map<pair<int,int>,int>yyb[32768];stack<int>stk;
namespace LCT{
typedef const int& fast;
int ch[maxn*2][2],fa[maxn*2],mx[maxn*2],w[maxn*2],index,A[maxn*2],B[maxn*2];bool rev[maxn*2];
il vd upd(fast x){
mx[x]=x;
if(w[mx[ch[x][0]]]>w[mx[x]])mx[x]=mx[ch[x][0]];
if(w[mx[ch[x][1]]]>w[mx[x]])mx[x]=mx[ch[x][1]];
}
il vd Rev(fast x){if(x)rev[x]^=1,swap(ch[x][0],ch[x][1]);}
il vd down(fast x){if(rev[x])Rev(ch[x][0]),Rev(ch[x][1]),rev[x]=0;}
il bool isrt(fast x){return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;}
il vd rotate(fast x){
sta int y,z,o;y=fa[x],z=fa[y],o=x==ch[y][1];
if(!isrt(y))ch[z][y==ch[z][1]]=x;fa[x]=z;
ch[y][o]=ch[x][!o],fa[ch[x][!o]]=y;
ch[x][!o]=y,fa[y]=x;
upd(y);
}
il vd splay(fast x){
sta int y,z,stk[maxn],top;stk[top=1]=x;
for(rg int i=x;!isrt(i);i=fa[i])stk[++top]=fa[i];
while(top)down(stk[top--]);
for(y=fa[x],z=fa[y];!isrt(x);rotate(x),y=fa[x],z=fa[y])
if(!isrt(y))rotate(((x==ch[y][0])^(y==ch[z][0]))?y:x);
upd(x);
}
il vd access(int x){for(rg int y=0;x;x=fa[y=x])splay(x),ch[x][1]=y,upd(x);}
il vd makert(fast x){access(x),splay(x),Rev(x);}
il vd link(fast x,fast y){makert(x),fa[x]=y;}
il vd cut(fast x,fast y){makert(x),access(y),splay(y);fa[x]=ch[y][0]=0;}
il int Query(fast x,fast y){makert(x),access(y),splay(y);return mx[y];}
il int find(int x){access(x),splay(x);while(ch[x][0])x=ch[x][0];return x;}
il vd Link(fast x,fast y,fast _w){
if(find(x)^find(y)){w[++index]=_w,link(A[index]=x,index),link(B[index]=y,index);return;}
sta int z;z=Query(x,y);
if(_w>=w[z])return;
cut(A[z],z),cut(B[z],z);
w[z]=_w,link(A[z]=x,z),link(B[z]=y,z);
}
};
int opt[maxn],A[maxn],B[maxn],C[maxn];
int X[10000001],Y[10000001],Z[10000001];
int main(){
freopen("tube_strong.in","r",stdin);
freopen("tube_strong.out","w",stdout);
int n=gi(),m=gi(),q=gi(),k;
LCT::index=n;
for(rg int i=1;i<=m;++i){X[i]=gi(),Y[i]=gi(),Z[i]=gi();if(X[i]>Y[i])swap(X[i],Y[i]);}
for(rg int i=1;i<=q;++i){
opt[i]=gi(),A[i]=gi(),B[i]=gi();if(A[i]>B[i])swap(A[i],B[i]);
if(opt[i]==2)yyb[(A[i]+B[i]*23333ll)&32767][make_pair(A[i],B[i])]=-1;
}
for(rg int i=1;i<=m;++i){
sta map<pair<int,int>,int>::iterator it;
k=(X[i]+Y[i]*23333ll)&32767;
it=yyb[k].find(make_pair(X[i],Y[i]));
if(it==yyb[k].end())LCT::Link(X[i],Y[i],Z[i]);
else it->second=Z[i];
}
while(q){
if(opt[q]==1)stk.push(LCT::w[LCT::Query(A[q],B[q])]);
else
LCT::Link(A[q],B[q],yyb[(A[q]+B[q]*23333ll)&32767][make_pair(A[q],B[q])]);
--q;
}
while(!stk.empty())printf("%d\n",stk.top()),stk.pop();
return 0;
}

//orz各位巨犇

cogs1885 [WC2006]水管局长数据加强版的更多相关文章

  1. BZOJ2594: [Wc2006]水管局长数据加强版

    题解: 裸LCT+离线+二分+MST... 代码:(几乎摘抄自hzwer) #include<cstdio> #include<cstdlib> #include<cma ...

  2. bzoj 2594: [Wc2006]水管局长数据加强版 动态树

    2594: [Wc2006]水管局长数据加强版 Time Limit: 25 Sec  Memory Limit: 128 MBSubmit: 934  Solved: 291[Submit][Sta ...

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

    离线然后就是维护加边的动态MST, Link cut tree秒掉..不过我写+调了好久...时间复杂度O(NlogN + MlogM) ------------------------------- ...

  4. [bzoj2594][Wc2006]水管局长数据加强版 (lct)

    论蒟蒻的自我修养T_T.. 和noi2014魔法森林基本一样...然而数据范围大得sxbk...UPD:这题如果用lct判联通的话可能会被卡到O(mlogm)..所以最好还是用并查集吧 一开始数组开太 ...

  5. BZOJ 2594: [Wc2006]水管局长数据加强版 [LCT kruskal]

    2594: [Wc2006]水管局长数据加强版 Time Limit: 25 Sec  Memory Limit: 128 MBSubmit: 2917  Solved: 918[Submit][St ...

  6. BZOJ_2594_[Wc2006]水管局长数据加强版_LCT

    BZOJ_2594_[Wc2006]水管局长数据加强版_LCT Description SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供 ...

  7. [WC2006]水管局长数据加强版

    Description SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一 ...

  8. 【刷题】BZOJ 2594 [Wc2006]水管局长数据加强版

    Description SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一 ...

  9. 2594. [WC2006]水管局长数据加强版【LCT+最小生成树】

    Description SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一 ...

随机推荐

  1. Error loading XML document: dwz.frag.xml 处理方式

    问题:直接用IE打开index.html弹出一个对话框:Error loading XML document: dwz.frag.xml 方案一(已经验证): 转自:http://blog.csdn. ...

  2. 解析UIControl

    解析UIControl 从下图可以看出,UIControl继承自UIView,添加了响应事件功能. UIButton之所以能响应各种各样的事件是因为继承自UIControl 使用UIControl可以 ...

  3. UNIX高级环境编程(5)Files And Directories - 文件相关时间,目录文件相关操作

     1 File Times 每个文件会维护三个时间字段,每个字段代表的时间都不同.如下表所示: 字段说明: st_mtim(the modification time)记录了文件内容最后一次被修改的时 ...

  4. linux centos6.5 网络配置

    1.方法一.修改网络配置文件 ①cd /etc/sysconfig/network-scripts ②cp ifcfg-eth0  ./ifcfg-eth0.bak   //修改前先备份 第一个以太网 ...

  5. 铁乐学python_Day44_IO多路复用

    目录 IO模型介绍 阻塞IO(blocking IO) 非阻塞IO(non-blocking IO) 多路复用IO(IO multiplexing) 异步IO(Asynchronous I/O) IO ...

  6. .NET Reflector注册机激活方法

    .NET Reflector注册机是一款专门针对.NET Reflector(.NET反编译工具软件)而推出的一款破解辅助工具软件.因为官方破解版软件需要118美元才能用,不然只有14天的试用期,为此 ...

  7. ABAP很厉害是怎么一种体验?

    知乎上偶然看到这个问题,觉得很有意思,我也来回答一发. 我本科和研究生学的是计算机专业,做项目用C/C++,研究生三年项目的代码量大概在三到四万行左右.2007年大学毕业加入SAP成都研究院一直工作到 ...

  8. sqlserver 一个字段分成两个字段 截取,保存

      Ma,LeiTo(Ma,lT) 分割后为 姓       名 Ma,LeiTo  Ma,lT SELECT LEFT(作者, CHARINDEX('(', 作者) - 1) AS 姓 , repl ...

  9. php 实现hash表

    hash表又称散列表,通过把关键字key经过hash函数映射到hash表中某个位置获取记录. 存放记录的数组又称为hash表,映射函数称为hash函数 下面是php中实现hash表的方法 <?p ...

  10. 启动64位 IIS 32位应用程序的支持

    64位的系统可以让IIS在32位的环境下运行asp.net程序,设置方法如下: ASP.NET程序在编译的时候默认是Any CPU,即编译的程序可以在X86.X64系统平台上运行.若希望我们的ASP. ...