BZOJ 4644
可以把修改当成删除再插入一个新的,
线性基不容易删除,就用线段树分治
好像手写的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的更多相关文章
- [BZOJ]4644: 经典傻逼题
某天我觉得一切题目都是那么不可做,于是百度了一下"傻逼题"-- 题目大意:对于图中的任意一个点集(可以为空或者全集),所有恰好有一个端点在这个点集中的边组成的集合被称为割.一个割的 ...
- 解题:BZOJ 4644 经典砂比题(雾
题面 初见线段树分治 (对我来说可不是什么经典题=.=) 把时间轴建出来一棵线段树,然后在对应的区间上打标记,最后把整棵树DFS一遍,到叶节点输出答案即可 (把最终答案开成全局的了调了半天 #incl ...
- 【BZOJ】【4146】 【AMPPZ2014】Divisors
暴力 由于值的范围很小($ \leq 2*10^6$),所以用一个cnt数组统计每个值有多少个数,然后从小到大,统计每个数的倍数即可. 根据调和数?的神奇性质= =这样是$O(nlogn)$的…… / ...
- BZOJ 2127: happiness [最小割]
2127: happiness Time Limit: 51 Sec Memory Limit: 259 MBSubmit: 1815 Solved: 878[Submit][Status][Di ...
- BZOJ 3275: Number
3275: Number Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 874 Solved: 371[Submit][Status][Discus ...
- BZOJ 2879: [Noi2012]美食节
2879: [Noi2012]美食节 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 1834 Solved: 969[Submit][Status] ...
- bzoj 4610 Ceiling Functi
bzoj 4610 Ceiling Functi Description bzoj上的描述有问题 给出\(n\)个长度为\(k\)的数列,将每个数列构成一个二叉搜索树,问有多少颗形态不同的树. Inp ...
- BZOJ 题目整理
bzoj 500题纪念 总结一发题目吧,挑几道题整理一下,(方便拖板子) 1039:每条线段与前一条线段之间的长度的比例和夹角不会因平移.旋转.放缩而改变,所以将每条轨迹改为比例和夹角的序列,复制一份 ...
- 【sdoi2013】森林 BZOJ 3123
Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数.第三行包含N个非负整数 ...
随机推荐
- java.lang.OutOfMemoryError: Java heap space内存不足问题
今晚,在定义一个new int[19001][13001]的数组时候内存不够:特转了一下方法: Exception in thread "main" java.lang.OutOf ...
- Redis学习笔记(三)使用Lua脚本实现分布式锁
Redis在2.6推出了脚本功能,允许开发者使用Lua语言编写脚本传到Redis中执行. 使用Lua脚本的好处如下: 1.减少网络开销:本来5次网络请求的操作,可以用一个请求完成,原先5次请求的逻辑放 ...
- LFYZ-OJ ID: 1009 阶乘和
思路 循环n次,每次计算i的阶乘i!,并加入sum中. n的范围从1~100,这里一定要使用高精度运算,涉及到"高精度乘低精度","高精度加高精度". 避免每次 ...
- 二维、三维 Laplace 算子的极坐标表示
(1) 设 $(r,\theta)$ 是 $\bbR^2$ 的极坐标, 即 $$\bex x=r\cos\theta,\quad y=r\sin \theta. \eex$$ 证明 Laplace 算 ...
- mysql中使用enum,如何获取所有可能的值
SELECT column_type FROM information_schema. COLUMNS WHERE TABLE_SCHEMA = "数据库名" AND DATA_T ...
- mysql root password
"""centos:mysql忘记root密码解决 1.修改MySQL的登录设置: # vim /etc/my.cnf 在[mysqld]的段中加上一句:skip-gra ...
- 项目改bug期间总结
kotlin中使用labda 表到时注意参数,另外setOnClickListener() 和 setOnClickListener{} 差别很大的,使用要细心, 因为这个问题拿了大半天闷, 最后还是 ...
- mybatis LIKE动态参数 sql语句
@Select({ "select id, vedio_name, vedio_path,vedio_duration, vedio_classify_id, crt_user_id, cr ...
- Android app中的so库和CPU架构
一.android目前有几种cpu架构? 早期的Android系统几乎只支持ARMv5的CPU架构,目前支持七种CPU架构:ARMv5,ARMv7 (从2010年起),x86 (从2011年起),MI ...
- ES6 + Webpack + React + Babel 如何在低版本浏览器上愉快的玩耍
https://blog.csdn.net/a324539017/article/details/52824189