[JZOJ4962]J

题目大意

一个$n$行$m$列的棋盘，这$n\times m$个位置全部放好了棋。棋有正反两面（用$0/1$表示）。
有总共$n+m$种魔法，前$i(1\le i\le n)$种魔法的效果是将第$i$行的所有棋翻面，第$i+n(1\le i\le m)$种魔法的效果是将第$i$列的所有棋翻面。当棋盘上形成某一个给定的形状时，便会触发一个魔法阵。
现在有$R$个人，这些人会依次触发这些魔法中的一种（两个人的魔法不一定相同），如果在某个人触发魔法后触发了魔法阵，那么这个人就会被嘿嘿嘿。
现在假定每个人最开始都只会触发第一种魔法，所有棋子都是$0$。有$Q$个操作，可能是下面三种中的一种：
$\bullet$把第$x$个人触发的魔法改为第$y$个魔法。
$\bullet$询问第$l$个人到第$r$个人中有多少个会被嘿嘿嘿。注意即使只是询问着一段，在这一段前面的人也会触发魔法。
$\bullet$把第$l$个人到第$r$个人触发的魔法全部修改为第$x$种。

对于每一个询问你都要输出答案。
$1\le n\le2,1\le m\le3,1\le R\le10^6,1\le Q\le1.2\times10^5$。

---

题目分析

注意到棋盘的大小其实很小，棋盘的形态我们可以使用二进制状态来表示。
那么现在相当于询问有多少个位置前缀异或和为某个数。修改有单点和区间两种。显然我们可以使用分块，每个块内维护块内前缀异或和为某个数的前缀位置个数。知道了是分块以后各种操作都很简单，在这里不累赘了。
时间复杂度$\mathrm O(Q\sqrt R+2^{nm}Q)$。

---

代码实现

#include <iostream>
#include <cstdio>
#include <cctype>
#include <cmath>

using namespace std;

int read()
{
    int x=0,f=1;
    char ch=getchar();
    while (!isdigit(ch)) f=ch=='-'?-1:f,ch=getchar();
    while (isdigit(ch)) x=x*10+ch-'0',ch=getchar();
    return x*f;
}

const int MAXR=1000050;
const int B=2005;
const int S=1<<6;
const int MG=6;

int size[B],tag[B],xsum[B];
int magtp[MAXR];
int magic[MG+5];
int ways[B][S];
int n,m,R,Q,bs,cnt,aim,all;

void pre()
{
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            magic[i]+=1<<(i-1)*m+j-1;
    for (int i=1;i<=m;i++)
        for (int j=1;j<=n;j++)
            magic[i+n]+=1<<(j-1)*m+i-1;
    for (int i=1;i<=R;i++) magtp[i]=1;
    bs=trunc(sqrt(R))+1;
    for (int i=1;i*bs<=R;i++) size[++cnt]=bs,ways[cnt][0]=size[cnt]-(ways[cnt][magic[1]]=size[cnt]+1>>1),xsum[cnt]=magic[1]*(size[cnt]&1);
    if (R%bs) size[++cnt]=R%bs,ways[cnt][0]=size[cnt]-(ways[cnt][magic[1]]=size[cnt]+1>>1),xsum[cnt]=magic[1]*(size[cnt]&1);
}

void rebuild(int bid)
{
    for (int s=0;s<all;s++) ways[bid][s]=0;
    int st=(bid-1)*bs+1,en=min(bid*bs,R);
    xsum[bid]=0;
    for (int i=st;i<=en;i++) magtp[i]=tag[bid]?tag[bid]:magtp[i],++ways[bid][xsum[bid]^=magic[magtp[i]]];
    tag[bid]=0;
}

int query(int l,int r)
{
    int st=0,en=0,xs=0,ret=0;
    for (int bid=1;bid<=cnt;bid++)
    {
        st=en+1,en=st+size[bid]-1;
        if (r<st) break;
        if (st<l&&l<=en)
        {
            if (tag[bid]) rebuild(bid);
            for (int i=st;i<l;i++) xs^=magic[magtp[i]];
            if (r<=en)
            {
                for (int i=l;i<=r;i++) ret+=(xs^=magic[magtp[i]])==aim;
                break;
            }
            else for (int i=l;i<=en;i++) ret+=(xs^=magic[magtp[i]])==aim;
            continue;
        }
        else if (st<=r&&r<en)
        {
            if (tag[bid]) rebuild(bid);
            for (int i=st;i<=r;i++) ret+=(xs^=magic[magtp[i]])==aim;
            break;
        }
        else if (l<=en) ret+=tag[bid]?((aim^xs)==magic[tag[bid]]?(size[bid]+1>>1):(!(aim^xs)?size[bid]-(size[bid]+1>>1):0)):ways[bid][aim^xs];
        xs^=tag[bid]?(magic[tag[bid]]*(size[bid]&1)):xsum[bid];
    }
    return ret;
}

void modify(int l,int r,int x)
{
    int st=0,en=0,xs=0,ret=0;
    for (int bid=1;bid<=cnt;bid++)
    {
        st=en+1,en=st+size[bid]-1;
        if (en<l) continue;
        if (r<st) break;
        if (st<l&&l<=en)
        {
            if (tag[bid]) rebuild(bid);
            for (int i=l;i<=en&&i<=r;i++) magtp[i]=x;
            rebuild(bid);
        }
        else if (st<=r&&r<en)
        {
            if (tag[bid]) rebuild(bid);
            for (int i=st;i<=r;i++) magtp[i]=x;
            rebuild(bid);
            break;
        }
        else tag[bid]=x;
    }
}

int main()
{
    freopen("ji.in","r",stdin),freopen("ji.out","w",stdout);
    n=read(),m=read(),all=1<<n*m;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            aim+=read()*(1<<(i-1)*m+j-1);
    for (R=read(),Q=read(),pre();Q--;)
    {
        int op=read();
        switch (op)
        {
            case 0:
            {
                int x=read(),y=read();
                rebuild((x-1)/bs+1),magtp[x]=y,rebuild((x-1)/bs+1);
                break;
            }
            case 1:
            {
                int l=read(),r=read();
                printf("%d\n",query(l,r));
                break;
            }
            case 2:
            {
                int l=read(),r=read(),x=read();
                modify(l,r,x);
                break;
            }
        }
    }
    fclose(stdin),fclose(stdout);
    return 0;
}