可以把修改当成删除再插入一个新的,

线性基不容易删除,就用线段树分治

好像手写的bitset在这里很慢

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=(a),i##_end=(b);i<=i##_end;++i)
#define For(i,a,b) for(int i=(a),i##_end=(b);i<i##_end;++i)
#define per(i,a,b) for(int i=(b),i##_st=(a);i>=i##_st;--i)
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define dbg(x) cerr<<#x" = "<<x<<endl
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define Es(x,i) for(Edge*i=G[x];i;i=i->nxt)
typedef long long ll;
typedef pair<int,int> pii;
const int inf=~0u>>1,MOD=1e9+7;
inline int rd() {
int x,c,f=1;while(!isdigit(c=getchar()))f=c!='-';x=c-'0';
while(isdigit(c=getchar()))x=x*10+c-'0';return f?x:-x;
} const int N=1011,MX=16;
typedef unsigned long long ui;
struct bs{
#define Q 64
ui w[MX];
bs (int _t=1){if(!_t)memset(w,0,sizeof w);}
bool operator<(const bs&t)const{
per(i,0,MX-1)if(w[i]!=t.w[i])return w[i]<t.w[i];
}
void operator^=(const bs&b){
For(i,0,MX)this->w[i]^=b.w[i];
}
bs operator^(const bs&b)const{
static bs c;
memset(c.w,0,sizeof c.w);
For(i,0,MX)c.w[i]=w[i]^b.w[i];
return c;
}
inline bool co(int i){
return w[i>>6]>>(i&63)&1;
}
inline void print(){
static int p[Q],i,j;
ui x;
for(i=MX-1;~i&&!w[i];--i);
if(i==-1){
putchar('0');
}
else{
j=0,x=w[i];
while(x)p[j++]=(x&1),x>>=1;
per(k,0,j-1)putchar(p[k]+'0');
for(--i;i>=0;--i){
ui x=w[i];
For(k,0,Q)p[k]=x&1,x>>=1;
per(k,0,Q-1)putchar(p[k]+'0');
}
}
putchar('\n');
}
#undef Q
};
int LL;
struct xxk{
bs a[N];
bool b[N];
xxk(int _t=1){if(!_t)memset(a,0,sizeof a);}
void add(bs x){
per(i,0,LL)if(x.co(i)){
if(!b[i]){
b[i]=1,a[i]=x;
break;
}else x^=a[i];
}
}
void qmx(){
static bs tmp,p;
memset(tmp.w,0,sizeof tmp.w);
per(i,0,LL)if(b[i]){
p=tmp^a[i];
if(tmp<p)tmp=p;
}
tmp.print();
}
} A[N<<2];
int ql,qr,qx;
bs ooo,p[N];
vector<bs> q[N<<2];
inline void tran(int o,int l,int r){
if(ql<=l&&r<=qr){
q[o].pb(p[qx]);
return;
}
int m=l+r>>1;
if(ql<=m)tran(o<<1,l,m);
if(qr>m)tran(o<<1|1,m+1,r);
}
inline void solve(int o,int l,int r){
A[o]=A[o>>1];
for(vector<bs>::iterator it=q[o].begin();it!=q[o].end();++it){
A[o].add(*it);
}
if(l==r){
A[o].qmx();
return;
}
int m=l+r>>1;
solve(o<<1,l,m),solve(o<<1|1,m+1,r);
}
char s[N];
int n,m;
int b[N];
int main(){
#ifdef flukehn
freopen("test.txt","r",stdin);
#endif
ios::sync_with_stdio(0),cin.tie(0);
cin>>n;
cin>>n>>m;
rep(i,1,n)b[i]=1;
rep(i,1,m){
int u,v;
cin>>u>>v>>s;
if(u==v)continue;
memset(ooo.w,0,sizeof ooo.w);
int j=strlen(s)-1,p=0;
LL=max(LL,j);
while(j>=63){
ui &x=ooo.w[p];
rep(k,j-64+1,j)x=(x<<1)+(s[k]-'0');
j-=64,++p;
}
ui&x=ooo.w[p];
rep(k,0,j)x=(x<<1)+(s[k]-'0');
if(b[u]){
ql=b[u],qr=i-1,qx=u;
if(i!=1)tran(1,1,m);
b[u]=i,::p[u]^=ooo;
}
if(b[v]){
ql=b[v],qr=i-1,qx=v;
if(i!=1)tran(1,1,m);
b[v]=i,::p[v]^=ooo;
}
}
rep(i,1,n){
ql=b[i],qr=m,qx=i;
tran(1,1,m);
}
solve(1,1,m);
}

BZOJ 4644的更多相关文章

  1. [BZOJ]4644: 经典傻逼题

    某天我觉得一切题目都是那么不可做,于是百度了一下"傻逼题"-- 题目大意:对于图中的任意一个点集(可以为空或者全集),所有恰好有一个端点在这个点集中的边组成的集合被称为割.一个割的 ...

  2. 解题:BZOJ 4644 经典砂比题(雾

    题面 初见线段树分治 (对我来说可不是什么经典题=.=) 把时间轴建出来一棵线段树,然后在对应的区间上打标记,最后把整棵树DFS一遍,到叶节点输出答案即可 (把最终答案开成全局的了调了半天 #incl ...

  3. 【BZOJ】【4146】 【AMPPZ2014】Divisors

    暴力 由于值的范围很小($ \leq 2*10^6$),所以用一个cnt数组统计每个值有多少个数,然后从小到大,统计每个数的倍数即可. 根据调和数?的神奇性质= =这样是$O(nlogn)$的…… / ...

  4. BZOJ 2127: happiness [最小割]

    2127: happiness Time Limit: 51 Sec  Memory Limit: 259 MBSubmit: 1815  Solved: 878[Submit][Status][Di ...

  5. BZOJ 3275: Number

    3275: Number Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 874  Solved: 371[Submit][Status][Discus ...

  6. BZOJ 2879: [Noi2012]美食节

    2879: [Noi2012]美食节 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 1834  Solved: 969[Submit][Status] ...

  7. bzoj 4610 Ceiling Functi

    bzoj 4610 Ceiling Functi Description bzoj上的描述有问题 给出\(n\)个长度为\(k\)的数列,将每个数列构成一个二叉搜索树,问有多少颗形态不同的树. Inp ...

  8. BZOJ 题目整理

    bzoj 500题纪念 总结一发题目吧,挑几道题整理一下,(方便拖板子) 1039:每条线段与前一条线段之间的长度的比例和夹角不会因平移.旋转.放缩而改变,所以将每条轨迹改为比例和夹角的序列,复制一份 ...

  9. 【sdoi2013】森林 BZOJ 3123

    Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数.第三行包含N个非负整数 ...

随机推荐

  1. 11.11CTF答案

    签到题                             flag{0ca175b9c0f726a831d895e269332461 } 下载WIN HEX,用WIN HEX打开“海贼王图片”, ...

  2. halcon预处理函数

    颜色空间变换 trans_from_rgb(ImageRed, ImageGreen, ImageBlue :ImageResult1, ImageResult2, ImageResult3 : Co ...

  3. kubernetes云平台管理实战: 故障自愈实战(四)

    一.创建实验文件 [root@k8s-master ~]# cat myweb-rc.yml apiVersion: v1 kind: ReplicationController metadata: ...

  4. Dubbo学习笔记8:Dubbo的线程模型与线程池策略

    Dubbo默认的底层网络通讯使用的是Netty,服务提供方NettyServer使用两级线程池,其中 EventLoopGroup(boss) 主要用来接受客户端的链接请求,并把接受的请求分发给 Ev ...

  5. cmd 命令添加防火墙端口

    windows dos 命令添加防火墙端口. 示例 123 端口: netsh firewall add portopening protocol = UDP port = name = NTPSER ...

  6. 题解-HNOI2017 抛硬币

    Problem loj2023 题意概述:甲抛掷 \(a\) 次硬币,乙抛掷 \(b\) 次硬币,问有多少种情况甲正面向上的次数比乙多,答案对 \(10^k\) 取模 对于 \(10\%\) 的数据, ...

  7. ActiveMQ依赖JDK版本关系

    1.如何查看官方发布的activeMQ依赖的JDK版本1)以ActiveMQ 5.15.2 Release为例:在下载页面的Change Log处, 2)打开下载号的jar包,以activemq-al ...

  8. 六.ansible批量管理服务

    期中集群架构-第六章-ansible批量管理服务介绍====================================================================== 01. ...

  9. 【python】多进程共享变量Manager

    Manager的复杂结构赋值问题 Manager的字典类型: 如果value是简单类型,比如int,可以直接赋值给共享变量,并可以后续直接修改 如果value是复杂类型 ,比如list,dict,则必 ...

  10. Web Application Vulnerablities

    1. File inclusion berfoe start this caption  i make a conclusion for install third-part as follow I ...