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个非负整数 ...
随机推荐
- host-only局域网络
在网桥配置的情况下设置: 网桥配置链接 配置前: 配置后: vi /etc/sysconfig/network-scripts/ifcfg-eth0 这里是使用static静态配置,具体网络配置及声明 ...
- ACM-ICPC 2018 徐州赛区网络预赛 J Maze Designer(最大生成树+LCA)
https://nanti.jisuanke.com/t/31462 题意 一个N*M的矩形,每个格点到其邻近点的边有其权值,需要构建出一个迷宫,使得构建迷宫的边权之和最小,之后Q次查询,每次给出两点 ...
- 解压unzip的用法
1.把文件解压到当前目录下 [root@instance-q6z0ksfb xmerge_test]# unzip db1.zip 2.把文件解压到指定的目录下,需要用到-d参数. unzip -d ...
- Element老司机开车了
orm Select选择器 坑: 1.选择器视图层一直渲染最后一个元素(value-key作为唯一标识符是关键)2.视图不更新,没有在data函数中声明变量,或者其他地方重置了已经声明过得变量 o ...
- Map和Collection
Map:key---Value(一对儿数据) HashMap:无序存放,key不允许重复 HashTable:无序存放,key不允许重复 key是set集合,value是collection集合 Co ...
- L1-Day4
L1-Day4 1.这消息使她非常悲伤. [我的翻译]The message makes she very sad. [标准答案]The news made her very sad. [对比分析]( ...
- django drf 基础学习2
DRF基本程序调用一 models初步编写 1 编写model.py from django.db import models 导入 class dbinfo(models.Model) ...
- 人体姿势识别,Convolutional pose machines文献阅读笔记。
开源实现 https://github.com/shihenw/convolutional-pose-machines-release(caffe版本) https://github.com/psyc ...
- Quartz+TopShelf实现定时任务
转自 https://www.cnblogs.com/frozenzhang/archive/2016/04/29/5443778.html 1.创建控制台程序 2.添加引用 添加TopShelf的引 ...
- unity读取txt,并且生成vector3的数组
读取txt的脚本: public void ReadFileTX(){ var str = File.ReadAllText(Application.streamingAssetsPath + &qu ...