Cod sursa(job #2450266)

Utilizator valentinchipuc123Valentin Chipuc valentinchipuc123 Data 22 august 2019 14:44:20
Problema BMatrix Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 5.35 kb
#include <iostream>
#include <fstream>
#include <string>
using namespace std;

int a[205][205],n,m,s[205],d[205],stiva[205],v[205][205],nr,mat[205][205],arie1=0,arie2=0,l1,l2,c1,c2,arie=0;
int mat1[205][205],mat2[205][205],mat3[205][205],mat4[205][205];
char ata[205][205];
int main()
{
ifstream f("bmatrix.in");
ofstream g("bmatrix.out");
f>>n>>m;
for(int i=1;i<=n;i++)
  for(int j=1;j<=m;j++)
 {
  f>>ata[i][j];
  if(ata[i][j]=='0') a[i][j]=0;
  else a[i][j]=1;
 }
for(int i=1;i<=n;i++)
 {
  for(int j=1;j<=m;j++)
    if(a[i][j]==0) v[i][j]=v[i-1][j]+1;
    else v[i][j]=0;
  for(int j=1;j<=m;j++)
   {
    if(j==1) {nr=1; stiva[1]=1;}
    else {
          nr++;
          while(v[i][stiva[nr-1]]>v[i][j])
         {
          nr--;
          d[stiva[nr]]=j;
         }
         stiva[nr]=j;
         }
   }
  for(int j=1;j<=nr;j++)
  d[stiva[j]]=m+1;
  for(int j=m;j>=1;j--)
  {
   if(j==m) {nr=1; stiva[1]=m;}
   else {
         nr++;
         while(v[i][stiva[nr-1]]>v[i][j])
         {
          nr--;
          s[stiva[nr]]=j;
         }
         stiva[nr]=j;
        }
  }
  for(int j=1;j<=nr;j++)
  s[stiva[j]]=0;
  for(int j=1;j<=m;j++)
    mat1[i][j]=v[i][j]*(d[j]-s[j]-1);
 }
for(int i=n;i>=1;i--)
 {
  for(int j=1;j<=m;j++)
    if(a[i][j]==0) v[i][j]=v[i+1][j]+1;
    else v[i][j]=0;
  for(int j=1;j<=m;j++)
   {
    if(j==1) {nr=1; stiva[1]=1;}
    else {
          nr++;
          while(v[i][stiva[nr-1]]>v[i][j])
         {
          nr--;
          d[stiva[nr]]=j;
         }
         stiva[nr]=j;
         }
   }
  for(int j=1;j<=nr;j++)
  d[stiva[j]]=m+1;
  for(int j=m;j>=1;j--)
  {
   if(j==m) {nr=1; stiva[1]=m;}
   else {
         nr++;
         while(v[i][stiva[nr-1]]>v[i][j])
         {
          nr--;
          s[stiva[nr]]=j;
         }
         stiva[nr]=j;
        }
  }
  for(int j=1;j<=nr;j++)
  s[stiva[j]]=0;
  for(int j=1;j<=m;j++)
    mat2[i][j]=v[i][j]*(d[j]-s[j]-1);
 }
for(int i=1;i<=m;i++)
 {
  for(int j=1;j<=n;j++)
    if(a[j][i]==0) v[j][i]=v[j][i-1]+1;
    else v[j][i]=0;
  for(int j=1;j<=n;j++)
   {
    if(j==1) {nr=1; stiva[1]=1;}
    else {
          nr++;
          while(v[stiva[nr-1]][i]>v[j][i])
         {
          nr--;
          d[stiva[nr]]=j;
         }
         stiva[nr]=j;
         }
   }
  for(int j=1;j<=nr;j++)
  d[stiva[j]]=n+1;
  for(int j=n;j>=1;j--)
  {
   if(j==n) {nr=1; stiva[1]=n;}
   else {
         nr++;
         while(v[stiva[nr-1]][i]>v[j][i])
         {
          nr--;
          s[stiva[nr]]=j;
         }
         stiva[nr]=j;
        }
  }
  for(int j=1;j<=nr;j++)
  s[stiva[j]]=0;
  for(int j=1;j<=n;j++)
    mat3[j][i]=v[j][i]*(d[j]-s[j]-1);
 }
for(int i=m;i>=1;i--)
 {
  for(int j=1;j<=n;j++)
    if(a[j][i]==0) v[j][i]=v[j][i+1]+1;
    else v[j][i]=0;
  for(int j=1;j<=n;j++)
   {
    if(j==1) {nr=1; stiva[1]=1;}
    else {
          nr++;
          while(v[stiva[nr-1]][i]>v[j][i])
         {
          nr--;
          d[stiva[nr]]=j;
         }
         stiva[nr]=j;
         }
   }
  for(int j=1;j<=nr;j++)
  d[stiva[j]]=n+1;
  for(int j=n;j>=1;j--)
  {
   if(j==n) {nr=1; stiva[1]=n;}
   else {
         nr++;
         while(v[stiva[nr-1]][i]>v[j][i])
         {
          nr--;
          s[stiva[nr]]=j;
         }
         stiva[nr]=j;
        }
  }
  for(int j=1;j<=nr;j++)
  s[stiva[j]]=0;
  for(int j=1;j<=n;j++)
    mat4[j][i]=v[j][i]*(d[j]-s[j]-1);
 }
for(int i=1;i<=n;i++)
 {
  for(int j=1;j<=m;j++)
   if(mat1[i][j]+mat2[i+1][j]>arie) arie=mat1[i][j]+mat2[i+1][j];
 }
for(int i=1;i<=m;i++)
 {
  for(int j=1;j<=n;j++)
   if(mat3[j][i]+mat4[j][i+1]>arie) arie=mat3[j][i]+mat4[j][i+1];
 }
for(int i=1;i<=n;i++)
 {
  for(int j=1;j<=m;j++)
    if(a[i][j]==0) v[i][j]=v[i-1][j]+1;
    else v[i][j]=0;
  for(int j=1;j<=m;j++)
   {
    if(j==1) {nr=1; stiva[1]=1;}
    else {
          nr++;
          while(v[i][stiva[nr-1]]>v[i][j])
         {
          nr--;
          d[stiva[nr]]=j;
         }
         stiva[nr]=j;
         }
   }
  for(int j=1;j<=nr;j++)
  d[stiva[j]]=m+1;
  for(int j=m;j>=1;j--)
  {
   if(j==m) {nr=1; stiva[1]=m;}
   else {
         nr++;
         while(v[i][stiva[nr-1]]>v[i][j])
         {
          nr--;
          s[stiva[nr]]=j;
         }
         stiva[nr]=j;
        }
  }
  for(int j=1;j<=m;j++)
   if(v[i][j]*(d[j]-s[j]-1)>arie1)
    {
     c1=s[j]+1;
     c2=d[j]-1;
     l1=i-v[i][j]+1;
     l2=i;
     arie1=v[i][j]*(d[j]-s[j]-1);
    }
 }
for(int i=l1;i<=l2;i++)
  for(int j=c1;j<=c2;j++)
 a[i][j]=1;
for(int i=1;i<=n;i++)
 {
  for(int j=1;j<=m;j++)
    if(a[i][j]==0) v[i][j]=v[i-1][j]+1;
    else v[i][j]=0;
  for(int j=1;j<=m;j++)
   {
    if(j==1) {nr=1; stiva[1]=1;}
    else {
          nr++;
          while(v[i][stiva[nr-1]]>v[i][j])
         {
          nr--;
          d[stiva[nr]]=j;
         }
         stiva[nr]=j;
         }
   }
  for(int j=1;j<=nr;j++)
  d[stiva[j]]=m+1;
  for(int j=m;j>=1;j--)
  {
   if(j==m) {nr=1; stiva[1]=m;}
   else {
         nr++;
         while(v[i][stiva[nr-1]]>v[i][j])
         {
          nr--;
          s[stiva[nr]]=j;
         }
         stiva[nr]=j;
        }
  }
  for(int j=1;j<=m;j++)
   if(v[i][j]*(d[j]-s[j]-1)>arie2)
    {
     arie2=v[i][j]*(d[j]-s[j]-1);
    }
 }
if(arie1+arie2>arie) arie=arie1+arie2;
g<<arie;
}