NOIP模拟38:a

#include<bits/stdc++.h>
using namespace std;
namespace STD
{
	#define rr register
	#define scanf ybbb=scanf
	typedef long long ll;
	const int M=5e4+4;
	const int N=32;
	int n,m,l,r,ybbb;
	char s[M];
	int mat[N][M];
	int ton[30*M];
	int read()
	{
		rr int x_read=0,y_read=1;
		rr char c_read=getchar();
		while(c_read<'0'||c_read>'9')
		{
			if(c_read=='-') y_read=-1;
			c_read=getchar();
		}
		while(c_read<='9'&&c_read>='0')
		{
			x_read=(x_read<<3)+(x_read<<1)+(c_read^48);
			c_read=getchar();
		}
		return x_read*y_read;
	}
};
using namespace STD;
int main()
{
	n=read(),m=read();
	for(rr int i=1;i<=n;i++)
	{	
		scanf("%s",s+1);
		for(rr int j=1;j<=m;j++)
			if(s[j]=='1')
				mat[i][j]=1;
	}
	l=read(),r=read();
	for(rr int i=1;i<=n;i++)
		for(rr int j=1;j<=m;j++)
			mat[i][j]+=mat[i][j-1];
	for(rr int i=1;i<=n;i++)
		for(rr int j=1;j<=m;j++)
			mat[i][j]+=mat[i-1][j];
	ll ans=0;
	for(rr int i=1;i<=n;i++)
		for(rr int j=i;j<=n;j++)
		{
			for(rr int k=1;k<=m;k++)
			{
				if(!l)
					ans+=ton[mat[j][k]-mat[i-1][k]];
				ton[mat[j][k]-mat[i-1][k]]++;
			}
			for(rr int k=1;k<=mat[j][m]-mat[i-1][m];k++)
				ton[k]+=ton[k-1];
			for(rr int k=1;k<=m;k++)
			{
				int temp=mat[j][k]-mat[i-1][k];
				if(l<=temp&&temp<=r) ans++;
				if(l<=temp)
				{
					if(l)
					{
						if(temp>r)
							ans+=ton[temp-l]-ton[temp-r-1];
						else 
							ans+=ton[temp-l];
					}
					else
					{
						if(temp>r)
							ans+=ton[temp-l-1]-ton[temp-r-1];
						else 
							ans+=ton[temp-l-1];
					}
				}
			}
			for(rr int k=0;k<=mat[j][m]-mat[i-1][m];k++)
				ton[k]=0;
		}
	printf("%lld\n",ans);
}