#ifndef Local
#pragma GCC optimize("O3,unroll-loops")
const int lim=2e3+5;
#else
const int lim=1e3+3;
#endif
#include "bits/stdc++.h"
using namespace std;
//#define int int64_t
#define pb push_back
int mod=1e9+7;
//const int mod=(1ll<<61)-1;
using pii=pair<int,int>;
char s[lim][lim];
int n,m,gn=0;
bool vis[lim*lim];
int parent[lim*lim];
int val[lim*lim/2],ll[lim*lim/2],rr[lim*lim/2];
inline void dfs(int node,int par){
vis[node]=1;
int x=node/lim,y=node%lim;
parent[node]=parent[par];
val[gn]+=s[x][y]-'0';
ll[gn]=min(ll[gn],y);
rr[gn]=max(rr[gn],y);
for(int dd:{1,-1,lim,-lim}){
if(!vis[node+dd]){
dfs(node+dd,node);
}
}
}
#define P(i,j) ((i)*lim+(j))
int g[lim][lim];
int dp[lim][lim],curl;
int ansss=INT_MIN;
inline void dnc(int l,int r,int optl,int optr){
int mid=(l+r)>>1;
int opt=min(optr,mid-1);
dp[curl][mid]=g[mid][opt]+dp[curl-1][opt];
for(int i=opt;optl<=i;i--){
if(dp[curl][mid]<g[mid][i]+dp[curl-1][i]){
dp[curl][mid]=g[mid][i]+dp[curl-1][i];
opt=i;
}
}
ansss=max(ansss,dp[curl][mid]);
if(l<=mid-1)dnc(l,mid-1,optl,opt);
if(mid+1<=r)dnc(mid+1,r,opt,optr);
}
struct inpout{
static const int sz=1<<17;
char inp[sz];
int curc=0,len=0;
inline char gc(){
if(curc==len){
len=(int)fread(inp,1,sz,stdin);
curc=0;
if(!len)return -1;
}
return inp[curc++];
}
inline void readstr(char*c,int ll){
for(int i=0;i<ll;i++){
c[i]=gc();
}
gc();
}
inline int readint(){
char c=gc();
bool isneg=(c=='-');
int ans=0;
if(!isneg)ans=c-'0';
c=gc();
while('0'<=c&&c<='9'){
ans=10*ans+c-'0';
c=gc();
}
if(isneg)return -ans;
return ans;
}
}io;
inline void solve(){
n=io.readint(),m=io.readint();
for(int i=0;i<m+2;i++){
s[0][i]=s[n+1][i]='.';
vis[P(0,i)]=vis[P(n+1,i)]=1;
}
for(int i=1;i<=n;i++){
vis[P(i,0)]=vis[P(i,m+1)]=1;
for(int j=1;j<=m;j++){
if((s[i][j]=io.gc())=='.'){
vis[P(i,j)]=1;
}
}
io.gc();
}
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
if(!vis[P(i,j)]){
val[gn]=0;
ll[gn]=INT_MAX;
rr[gn]=-1;
parent[P(i,j)]=gn;
dfs(P(i,j),-1);
gn++;
}
}
}
for(int i=0;i<gn;i++){
g[ll[i]][0]+=val[i];
g[ll[i]][ll[i]]-=val[i];
g[rr[i]+1][0]-=val[i];
g[rr[i]+1][ll[i]]+=val[i];
}
for(int i=0;i<=m;i++){
for(int j=1;j<=m;j++){
g[i][j]+=g[i][j-1];
}
}
for(int i=0;i<=m;i++){
for(int j=1;j<=m;j++){
g[j][i]+=g[j-1][i];
}
}
dp[0][0]=0;
for(int i=1;i<=m;i++){
dp[0][i]=g[i][0];
ansss=max(ansss,dp[0][i]);
}
cout<<ansss<<"\n";
for(int i=1;i<m;i++){
ansss=INT_MIN;
curl=i;
dnc(i,m,1,m);
cout<<ansss<<"\n";
}
}
signed main(){
#ifdef Local
freopen(".in","r",stdin);
freopen(".out","w",stdout);
#endif
solve();
}
Compilation message
nafta.cpp: In function 'void solve()':
nafta.cpp:29:28: warning: array subscript -1 is below array bounds of 'int [4020025]' [-Warray-bounds]
29 | parent[node]=parent[par];
| ~~~~~~~~~~^
nafta.cpp:23:5: note: while referencing 'parent'
23 | int parent[lim*lim];
| ^~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
16988 KB |
Output is correct |
2 |
Correct |
3 ms |
16988 KB |
Output is correct |
3 |
Correct |
2 ms |
16988 KB |
Output is correct |
4 |
Correct |
3 ms |
16732 KB |
Output is correct |
5 |
Correct |
2 ms |
16732 KB |
Output is correct |
6 |
Correct |
2 ms |
16732 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
16988 KB |
Output is correct |
2 |
Correct |
3 ms |
16988 KB |
Output is correct |
3 |
Correct |
2 ms |
16988 KB |
Output is correct |
4 |
Correct |
3 ms |
16732 KB |
Output is correct |
5 |
Correct |
2 ms |
16732 KB |
Output is correct |
6 |
Correct |
2 ms |
16732 KB |
Output is correct |
7 |
Correct |
5 ms |
23132 KB |
Output is correct |
8 |
Correct |
5 ms |
23132 KB |
Output is correct |
9 |
Correct |
6 ms |
21596 KB |
Output is correct |
10 |
Correct |
5 ms |
23132 KB |
Output is correct |
11 |
Correct |
4 ms |
21084 KB |
Output is correct |
12 |
Correct |
4 ms |
21156 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
16988 KB |
Output is correct |
2 |
Correct |
3 ms |
16988 KB |
Output is correct |
3 |
Correct |
2 ms |
16988 KB |
Output is correct |
4 |
Correct |
3 ms |
16732 KB |
Output is correct |
5 |
Correct |
2 ms |
16732 KB |
Output is correct |
6 |
Correct |
2 ms |
16732 KB |
Output is correct |
7 |
Correct |
5 ms |
23132 KB |
Output is correct |
8 |
Correct |
5 ms |
23132 KB |
Output is correct |
9 |
Correct |
6 ms |
21596 KB |
Output is correct |
10 |
Correct |
5 ms |
23132 KB |
Output is correct |
11 |
Correct |
4 ms |
21084 KB |
Output is correct |
12 |
Correct |
4 ms |
21156 KB |
Output is correct |
13 |
Correct |
101 ms |
67200 KB |
Output is correct |
14 |
Correct |
120 ms |
63060 KB |
Output is correct |
15 |
Correct |
151 ms |
80696 KB |
Output is correct |
16 |
Correct |
77 ms |
63048 KB |
Output is correct |
17 |
Correct |
58 ms |
63060 KB |
Output is correct |
18 |
Correct |
62 ms |
63072 KB |
Output is correct |