行列接頭辞+と剪定の魔法
行列プレフィックスと
示すように、A [2] [2]及び接頭辞は2 + 3 + 2 + 7、 [1]〜[4] 及びプレフィックスが2 + 3 + 4 + 8である
ようにと、それぞれの座標接頭;マトリックス及び数の左上隅の座標である
;したがって部分行列の、任意の数とすることができ、クエリ内のO(1)時間
のコードは
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
b[i][j]+=b[i][j-1];
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
b[i][j]+=b[i-1][j];//计算前缀和
int js(int i,int j,int l,int r)//i,j为查询矩阵左上角顶点,l,r为右下角顶点
{
sum=b[l][r]-b[i-1][r]-b[l][j-1]+b[i-1][j-1];
}//查询任意矩阵数字和
だから、元はかなり速く実行されている非常に良いアルゴリズムではないこと。。。
エントリー
3 2
4 0
-10 8
-2 -2
出力
4
#include<cstdio>
#include<iostream>
using namespace std;
bool js(int i,int j,int l,int r)
{
return b[l][r]-b[i-1][r]-b[l][j-1]+b[i-1][j-1]>0;
}
bool check(int x)
{
for(int i=1;i<=n;i++)
{
if(x%i)continue;
int ll=i,rr=x/i;
for(int k=1;k+ll-1<=n;k++)
for(int h=1;h+rr-1<=n;h++)
if(js(k,h,k+ll-1,h+rr-1))
return 1;
}
for(int i=1;i<=n;i++)
pr[i*x]=1;
return 0;
}
int main()
{
scanf("%d%d",&n,&m)
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",a[i][j]);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
b[i][j]+=b[i][j-1];
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
b[i][j]+=b[i-1][j];
for(int i=1;i<=n*m;i++)
{
if(!pi[i])
if(check(i))
ans=i;
}
printf("%d",ans);
}