poj2155树状数组（二维）

poj2155解题报告

题目介绍：给定一个01矩阵（初始时为0），对它进行两种操作：

                1.  C：改变一个矩阵里的数（0变1,1变0）；

                2.  Q：查询一个点的值；

输入：第一行一个T，表示T组数据，对于每组数据，输入一组n和m，分别表示n*你的矩阵和m次操作，接下来m行，每行输入

          一个字符（C或Q），对于C，接着输入x1,x2,y1,y2（0<x1<=x2  ,  0<y1<=y2<n）表示矩阵；对于Q，输入x1，y1；

输出 : 对于每一个Q，输出矩阵中（x1，y1）表示的点上的值（0或1）；

 

解题思路：容斥原理，可以令sum[i][j]=a[1][1]+a[1][2]+……+a[i][j],所以进行更改时只需要对sum[x1][y1]，sum[x2+1][y1]，sum[x1][y2+1]，sum[x2+1][y2+1]进行修改，这个过程可以用树状数组来维护。将会对红色部分操作1次，蓝色2次，绿色4次（因为是01矩阵，所以相当于不变）。

代码如下：

#include<stdio.h>
#include<math.h>
#include<string.h>
#define MAXN 81
using namespace std;
typedef long long ll;
struct lll{
	char s[200];
};
lll sum[MAXN],f[MAXN][MAXN],pre[MAXN];
int m,n,a[MAXN];

template <class T>void read(T &x)
{
	x=0;int f=0;char ch=getchar();
	while(ch<'0'||ch>'9'){f|=(ch=='9'); ch=getchar();}
	while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
	x=f?-x:x;
}

lll max(lll a,lll b)
{
	int la=strlen(a.s),lb=strlen(b.s);
	if(la>lb)return a;
	if(lb>la)return b;
	for(int i=la-1;i>=0;i--)
	{
		if(a.s[i]>b.s[i])return a;
		if(a.s[i]<b.s[i])return b;
	}
	return a;
}

lll mul(int a,lll b)
{
	lll res;
	memset(res.s,'0',sizeof(res));
	char x[200];
	int la=0,lb=strlen(b.s),p;
	while(a)
	{
		x[la++]=a%10+48;
		a/=10;
	}
	for(int i=0;i<la;i++)
	{
		for(int j=0;j<lb;j++)
		{
			int u=(i+1)*(j+1);
			p=((x[i]-48)*(b.s[i]-48));
			res.s[u-1]+=((x[i]-48)*(b.s[i]-48));
			p=res.s[u-1]-48;
			if(res.s[u-1]>'9')
			{
				res.s[u-1]=p%10+48;
				res.s[u]+=p/10;
			}
		}
	}
	for(int i=0;i<=la*lb;i++)
	{
		if(res.s[i]>'9')
		{
			p=res.s[i]-48;
			res.s[i]=p%10+48;
			res.s[i+1]+=p/10;
		}
	}
	return res;
}

lll add(lll a,lll b)
{
	lll res;
	int la=strlen(a.s),lb=strlen(a.s),i,L;
	if(la>lb)L=lb;
	int rst=0;
	for(i=0;i<L;i++)
	{
		int p=a.s[i]+b.s[i]-96;
		if(p>9)res.s[i]=p%10+48,rst=p/10;
		else   res.s[i]=p+48,rst=0;
	}
	while(la<lb)
	{
		res.s[++la]=b.s[i]+rst;
		if(res.s[la]>'9')
		{
			rst=(res.s[la]-48)/10;
			res.s[la]=(res.s[la]-48)%10+48;
		}
	}
	while(lb<la)
	{
		res.s[++lb]=a.s[i]+rst;
		if(res.s[lb]>'9')
		{
			rst=(res.s[lb]-48)/10;
			res.s[lb]=(res.s[lb]-48)%10+48;
		}
	}
	return res;
}

int main()
{
	read(n),read(m);
	for(int i=2;i<=m;i++)
		pre[i]=mul(2,pre[i-1]);
	lll ans;
	strcpy(ans.s,"0");
	for(int w=1;w<=n;w++)
	{
		for(int i=1;i<=m;i++)
		{
			scanf("%d",a[i]);
			f[i][i]=mul(a[i],pre[m]);
		}
		for(int i=1;i<m;i++)
		{
			for(int j=1;j<=m-i;j++)
			{
				int k=i+j;
				f[j][k]=max(add(f[j+1][k],mul(a[j],pre[m-i])),add(f[j][k-1],mul(a[k],pre[m-i])));
			}
		}
		ans=add(ans,f[1][m]);
	}
	printf("%s\n",ans.s);
	return 0;
}