牛客71E 组一组

https://ac.nowcoder.com/acm/contest/71/E

看见是位运算，就知道按位考虑，看见一段连续的全是1，一段1的个数>=1，一段全是0，一段1的个数<r-l+1，就想到对位做前缀和的差分约束，然而对于差分约束的认识还是差了点。。

本来我以为由于可以确定哪些地方是1哪些地方是0，就按照单向边从小到大，就不用建回来的边，题目保证有解，那么这样就一定会出来一组合法解。

但通过这题的实验，我发现少了任何一种关系的建边，都会WA掉

而且这题如果直接跑还会超时，从yyb那里学了一手确定是1的地方i，直接从0连到i一条sum[i]的边，sum[i]表示已经有多少,这样i会在比较早的时候达到一个比较大的水平，进队列的次数就会少很多。

前缀和 dis[j]-dis[j-1]>=0  dis[j-1]-dis[j]>=-1

一段全是1  -> dis[r]-dis[l-1]>=r-l+1 dis[l-1]-dis[r]>=-(r-l+1)

一段至少有1个1  -> dis[r]-dis[l-1]>=1

一段全是0  ->dis[r]-dis[l-1]>=0  dis[l-1]-dis[r]>=0

一段不全是1 -> dis[l-1]-dis[r]>=l-r

这些连边一条也不能少，不然就会答案错误。。。

然后再加上用来优化的边才能过

#include<bits/stdc++.h>
using namespace std;

const int maxl=1e5+10;

int n,m,cnt;
int ans[maxl],sum[maxl],l[maxl],r[maxl],op[maxl],x[maxl];
int in1[maxl],out1[maxl],ehead[maxl],dis[maxl];
struct ed
{
	int to,nxt,l;
}e[maxl*5];
queue<int> q;
bool in[maxl];

inline void add(int u,int v,int l)
{
	e[++cnt].to=v;e[cnt].l=l;
	e[cnt].nxt=ehead[u];ehead[u]=cnt;
}

inline void prework()
{
	scanf("%d%d",&n,&m);
	//rd(n);rd(m);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d%d",&op[i],&l[i],&r[i],&x[i]);
		//rd(op[i]);rd(l[i]);rd(r[i]);rd(x[i]);
	}
}

inline void spfa(int id)
{
	for(int i=0;i<=n;i++)
		dis[i]=-2e9;
	int u,v;
	dis[0]=0;q.push(0);in[0]=true;
	while(!q.empty())
	{
		u=q.front();q.pop();
		for(int i=ehead[u];i;i=e[i].nxt)
		{
			v=e[i].to;
			if(dis[u]+e[i].l>dis[v])
			{
				dis[v]=dis[u]+e[i].l;
				if(!in[v])
				{
					in[v]=true;
					q.push(v);
				}
			}
		}
		in[u]=false;
	}
	for(int i=1;i<=n;i++)
	if(dis[i]>=dis[i-1]+1)
		ans[i]|=1<<id;
}

inline void mainwork()
{
	for(int i=0;i<20;i++)
	{
		for(int j=0;j<=n;j++)
			ehead[j]=0,in1[j]=0,out1[j]=0;
		cnt=0;
		for(int j=1;j<=m;j++)
		if(op[j]==1)
		{
			if(x[j]&(1<<i))
				add(l[j]-1,r[j],1);
			else
				add(l[j]-1,r[j],0),add(r[j],l[j]-1,0);
		}
		else
		{
			if(x[j]&(1<<i))
			{
				add(l[j]-1,r[j],r[j]-l[j]+1);
				add(r[j],l[j]-1,-(r[j]-l[j]+1));
				in1[l[j]]++,out1[r[j]]++;
			}
			else
				add(r[j],l[j]-1,l[j]-r[j]);	
		}
		int sum1=0;
		for(int j=1;j<=n;j++)
		{
			sum1+=in1[j];
			sum[j]=sum[j-1]+(sum1>0); 
			add(0,j,sum[j]);
			/*if(sum1>0)
				add(j-1,j,1);
			else
				add(j-1,j,0);*/
			add(j-1,j,0);
			add(j,j-1,-1);
			sum1-=out1[j];
		}
		spfa(i);
		for(int j=1;j<=n;j++)
		if(dis[j]>dis[j-1])
			ans[j]|=1<<i;
	}
}

inline void print()
{
	for(int i=1;i<=n;i++)
		printf("%d%c",ans[i],(i==n)?'\n':' ');
}

int main()
{
	prework();
	mainwork();
	print();
	return 0;
}