#include <bits/stdc++.h>
#define maxn (int)(1e5+51)
#define all(x) x.begin(), x.end()
#define sz(x) (int) x.size()
#define endl '\n'
#define ll long long
#define pb push_back
#define ull unsigned long long
#define ii pair<int,int>
#define iii tuple<int,int,int>
#define inf 2000000001
#define mod 1000000007 //998244353
#define _ ios_base::sync_with_stdio(false);cin.tie(NULL);cout.tie(NULL);
using namespace std;
template<typename X, typename Y> bool ckmin(X& x, const Y& y) { return (y < x) ? (x=y,1):0; }
template<typename X, typename Y> bool ckmax(X& x, const Y& y) { return (x < y) ? (x=y,1):0; }
int grid[1001][1001], pr[1001][1001], le[1001], ri[1001], n, m, k;
int best(int x){
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j)
pr[i][j] = grid[i][j]>=x?1:0;
// cout<<x<<endl;
// for(int i=1;i<=n;++i)
// for(int j=1;j<=m;++j)
// cout<<pr[i][j]<<" \n"[j==m];
for(int j=1;j<=m;++j)
for(int i=1;i<=n;++i)
pr[i][j]+=pr[i-1][j];
int resp = 0;
for(int i=1;i<=n;++i){
stack<ii>l, r;
for(int j=1;j<=m;++j){
while(!l.empty()&&l.top().first>pr[i][j])l.pop();
if(!l.empty()&&l.top().first==pr[i][j])le[j] = le[l.top().second] + j-l.top().second;
else if(l.empty())le[j] = j-1;
else le[j] = j - l.top().second - 1;
l.push({pr[i][j],j});
}
for(int j=m;j>=1;--j){
while(!r.empty()&&r.top().first>pr[i][j])r.pop();
if(!r.empty()&&r.top().first==pr[i][j])ri[j] = ri[r.top().second] + r.top().second-j;
else if(r.empty())ri[j] = m-j;
else ri[j] = r.top().second - j - 1;
r.push({pr[i][j],j});
}
for(int j=1;j<=m;++j){
// cout<<le[j]<<" \n"[j==m];
ckmax(resp,(le[j]+ri[j]+1)*pr[i][j]);
}
}
return resp;
}
int main(){_
// binary search the cheapest square
// grid[i][j] = 0 if c[i][j] < x, 1 if otherwise
// use monotonic stacks to find largest rectangle
// if size of largest rectangle >= k, increase l
// otherwise lower r
// total complexity: O(NMlog(1e9))
cin>>n>>m>>k;
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j)
cin>>grid[i][j];
int l = 1, r = 100, ans = 1, siz = 0;
while(l<=r){
int md = l + (r-l+1)/2;
int b = best(md);
if(b>=k)ans = md, siz = b, l = md + 1;
else r = md - 1;
}
cout<<ans<<" "<<siz<<endl;
// 0 0 0 0 0 0
// 0 1 0 0 0 0
// 0 0 1 1 0 0
// 0 1 1 1 0 0
// 0 1 1 0 0 0
// 0 0 0 0 0 0
}
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Incorrect |
1 ms |
308 KB |
Output isn't correct |
| 2 |
Halted |
0 ms |
0 KB |
- |
