hdu5465

这道题明显是博弈+二维线段树，但是不知道怎么写出了时间排名倒几的程序。

已ac的代码：

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
#define N 510
#pragma comment(linker, "/STACK:1024000000,1024000000")

int ornum[N*N*10];
int allx1,ally1,allx2,ally2;
int changeto;
int cnum[N][N];

void create(int root,int x1,int y1,int x2,int y2){
    if(x1==x2&&y1==y2){
        ornum[root]=cnum[x1][y1];
        //printf("%d %d %d %d %d\n",x1,y1,x2,y2,ornum[root]);
    }
    else{
        int xmiddle=(x1+x2)>>1;
        int ymiddle=(y1+y2)>>1;

        create((root<<2)+1,x1,y1,xmiddle,ymiddle);
        if(y1!=y2){
            create((root<<2)+2,x1,ymiddle+1,xmiddle,y2);
        }
        else{
            ornum[(root<<2)+2]=0;
        }

        if(x1!=x2){
            create((root<<2)+3,xmiddle+1,y1,x2,ymiddle);
        }
        else{
            ornum[(root<<2)+3]=0;
        }

        if(x1!=x2&&y1!=y2){
            create((root<<2)+4,xmiddle+1,ymiddle+1,x2,y2);
        }
        else{
            ornum[(root<<2)+4]=0;
        }

        ornum[root]=ornum[(root<<2)+1]^ornum[(root<<2)+2]^ornum[(root<<2)+3]^ornum[(root<<2)+4];//ornum[(root<<2)+2]写成了ornum[(root)<<2+2],于是就一直有问题

    }

    return;
}

int query(int root,int x1,int y1,int x2,int y2){
    if(allx1<=x1&&allx2>=x2&&ally1<=y1&&ally2>=y2){//把x，y和allx，ally之间的关系想反了，所以这里写反了，大于等于和小于等于左右两边写反了。
        return ornum[root];
    }
    else{
        int xmiddle=(x1+x2)>>1;
        int ymiddle=(y1+y2)>>1;
        int ans=0;

        if(allx1<=xmiddle){
            if(ally1<=ymiddle){
                ans^=query((root<<2)+1,x1,y1,xmiddle,ymiddle);
            }
            if(ally2>ymiddle&&y1!=y2){
                ans^=query((root<<2)+2,x1,ymiddle+1,xmiddle,y2);//这里写成了xmiddle+1,y1，所以一直栈溢出？
            }
        }
        if(allx2>xmiddle&&x1!=x2){
            if(ally1<=ymiddle){
                ans^=query((root<<2)+3,xmiddle+1,y1,x2,ymiddle);
            }
            if(ally2>ymiddle&&y1!=y2){
                ans^=query((root<<2)+4,xmiddle+1,ymiddle+1,x2,y2);
            }
        }

        return ans;
    }
}

void change(int root,int x1,int y1,int x2,int y2){
    if(x1==x2&&y1==y2){
        ornum[root]=changeto;
    }
    else{
        int xmiddle=(x1+x2)>>1;
        int ymiddle=(y1+y2)>>1;

        if(allx1<=xmiddle&&ally1<=ymiddle){
            change((root<<2)+1,x1,y1,xmiddle,ymiddle);
        }
        else if(allx1<=xmiddle&&ally1>ymiddle){
            change((root<<2)+2,x1,ymiddle+1,xmiddle,y2);
        }
        else if(allx1>xmiddle&&ally1<=ymiddle){
            change((root<<2)+3,xmiddle+1,y1,x2,ymiddle);
        }
        else{
            change((root<<2)+4,xmiddle+1,ymiddle+1,x2,y2);
        }

        ornum[root]=ornum[(root<<2)+1]^ornum[(root<<2)+2]^ornum[(root<<2)+3]^ornum[(root<<2)+4];
    }

    return;
}

int main(){
    int t;
    int n,m,q;
    int op;
    //int x1,x2,y1,y2;

    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d",&n,&m,&q);

        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&cnum[i][j]);//二维线段树不像一维的，可以在建树是输入，因为二维的如果这样做，会错乱。
            }
        }

        create(0,1,1,n,m);
        //printf("wo shi da hao ren");
        for(int i=0;i<q;i++){
            scanf("%d",&op);
            if(op==1){
                scanf("%d%d%d%d",&allx1,&ally1,&allx2,&ally2);

                //printf("%d %d %d %d\n",allx1,ally1,allx2,ally2);
                /*printf("Yes\n");*/
                //printf("%d %d\n",query(0,1,1,n,m),ornum[0]);
                if(query(0,1,1,n,m)==0){//猜测：query使得栈溢出
                    printf("No\n");
                }
                else{
                    printf("Yes\n");
                }
            }
            else{
                scanf("%d%d%d",&allx1,&ally1,&changeto);
                change(0,1,1,n,m);
            }
        }
    }

    return 0;
}

参考链接：

http://blog.csdn.net/dingding_tao/article/details/48596601

看别人的题解使用二维树状数组做的，于是就又写了一遍，结果发现我理解的二维数组只不过是很多一维数组的集合罢了，将近超时。

将近超时的ac的代码：

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
#define N 510

int bit[N][N];
int cnum[N][N];
int n,m;

int lowbit(int x){
    return x&((~x)+1);
}

void change(int x,int y,int from,int changeto){
    while(y<=m){
        bit[x][y]^=from;
        bit[x][y]^=changeto;
        y=y+lowbit(y);
    }

    return;
}

int query(int x1,int y1,int x2,int y2){
    int ans=0;

    for(int i=x1;i<=x2;i++){
        for(int j=y1-1;j>0;j-=lowbit(j)){
            ans^=bit[i][j];
        }
        for(int j=y2;j>0;j-=lowbit(j)){
            ans^=bit[i][j];
        }
    }

    return ans;
}

int main(){
    int t;
    int q;
    int tempnum;
    int op;
    int x1,y1,x2,y2;
    int ans;
    int changeto;

    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d",&n,&m,&q);

        memset(bit,0,sizeof(bit));
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&tempnum);
                change(i,j,0,tempnum);
                cnum[i][j]=tempnum;
            }
        }

        for(int i=0;i<q;i++){
            scanf("%d",&op);

            if(op==1){
                scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
                ans=0;

                ans=ans^query(x1,y1,x2,y2);
                if(ans==0){
                    printf("No\n");
                }
                else{
                    printf("Yes\n");
                }
            }
            else{
                scanf("%d%d%d",&x1,&y1,&changeto);
                change(x1,y1,cnum[x1][y1],changeto);
                cnum[x1][y1]=changeto;
            }
        }
    }

    return 0;
}

真正的二维数组，果然速度快好多。

ac的代码：

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
#define N 510

int bit[N][N];
int cnum[N][N];
int n,m;

int lowbit(int x){
    return ((~x)+1)&x;
}

int query(int x,int y){
    int ans=0;

    for(int i=x;i;i=i-lowbit(i)){
        for(int j=y;j;j=j-lowbit(j)){
            ans^=bit[i][j];
        }
    }

    return ans;
}

void change(int x,int y,int changeto){
    for(int i=x;i<=n;i=i+lowbit(i)){
        for(int j=y;j<=m;j=j+lowbit(j)){
            bit[i][j]^=changeto;
        }
    }

    return;
}

int main(){
    int t;
    int q;
    int op;
    int x1,y1,x2,y2;
    int ans;
    int changeto;

    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d",&n,&m,&q);

        memset(bit,0,sizeof(bit));
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&cnum[i][j]);
                change(i,j,cnum[i][j]);
            }
        }

        for(int i=0;i<q;i++){
            scanf("%d",&op);

            if(op==1){
                scanf("%d%d%d%d",&x1,&y1,&x2,&y2);

                ans=query(x2,y2)^query(x1-1,y2)^query(x2,y1-1)^query(x1-1,y1-1);
                if(ans==0){
                    printf("No\n");
                }
                else{
                    printf("Yes\n");
                }
            }
            else{
                scanf("%d%d%d",&x1,&y1,&changeto);
                change(x1,y1,changeto^cnum[x1][y1]);
                cnum[x1][y1]=changeto;
            }
        }
    }

    return 0;
}