#include <bits/stdc++.h>
using namespace std;
#define fast ios::sync_with_stdio(0);cin.tie(0);
#define s second
#define f first
typedef long long ll;
const ll MOD = 1e9 + 7;
const ll LOGN = 18;
const ll INF = 1e15;
const ll MAXN = 2005;
vector<int> occ[MAXN];
int N, M;
int emp[MAXN][MAXN], sg[4*MAXN], lazy[4*MAXN];
ll total = 0, ans = 0;
void push(int k, int a, int b) {
if (lazy[k] != -1) {
sg[k] = lazy[k];
if (a != b) {
lazy[2*k] = lazy[k];
lazy[2*k+1] = lazy[k];
}
lazy[k] = -1;
}
}
void update(int k, int a, int b, int q_l, int q_r, int new_val) {
push(k, a, b);
if (q_r < a || q_l > b)
return ;
if (q_l <= a && b <= q_r) {
lazy[k] = new_val;
push(k, a, b);
return ;
}
update(2*k, a, (a+b)/2, q_l, q_r, new_val);
update(2*k+1, (a+b)/2+1, b, q_l, q_r, new_val);
sg[k] = max(sg[2*k], sg[2*k+1]);
}
ll calc(ll x) {
return (x * x * x + 3 * x * x + 2 * x) / 6;
}
int get(int k, int a, int b, int pos) {
push(k, a, b);
if (b < pos || a > pos)
return 0;
if (a == b)
return sg[k];
return get(2*k, a, (a+b)/2, pos) + get(2*k+1, (a+b)/2+1, b, pos);
}
int query(int k, int a, int b, int val) {
push(k, a, b);
if (a == b) {
if (sg[k] >= val)
return a;
return N+1;
}
push(2*k, a, (a+b)/2);
if (sg[2*k] >= val)
return query(2*k, a, (a+b)/2, val);
return query(2*k+1, (a+b)/2 + 1, b, val);
}
void split(int pos) {
int Q = get(1, 0, N, pos);
int l = query(1, 0, N, Q);
int r = query(1, 0, N, Q + 1) - 1;
ans += calc(pos - l - 1) + calc(r - pos) - calc(r - l);
update(1, 0, N, pos, r, pos);
}
string s;
int main() {
fast
cin >> N >> M;
for (int i = 1; i <= N; i++) {
cin >> s;
for (int j = M; j >= 1; j--) {
if (s[j-1] == '.')
emp[i][j] = emp[i][j+1] + 1;
}
}
ll temp = calc(N);
for (int j = 1; j <= M; j++) {
memset(sg, 0, sizeof(sg));
memset(lazy, -1, sizeof(lazy));
ans = temp;
for (int i = 1; i <= N; i++)
occ[emp[i][j]].push_back(i);
for (ll q = 0; q <= 2000; q++) {
total += ans * q;
while (occ[q].size()) {
split(occ[q].back());
occ[q].pop_back();
}
}
}
cout << total << "\n";
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
600 KB |
Output is correct |
2 |
Correct |
1 ms |
604 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
604 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
28 ms |
2136 KB |
Output is correct |
2 |
Correct |
29 ms |
2140 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
28 ms |
2148 KB |
Output is correct |
2 |
Correct |
29 ms |
2140 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
21 ms |
2648 KB |
Output is correct |
2 |
Correct |
29 ms |
2136 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
351 ms |
7004 KB |
Output is correct |
2 |
Correct |
912 ms |
12500 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
871 ms |
11024 KB |
Output is correct |
2 |
Execution timed out |
1068 ms |
15992 KB |
Time limit exceeded |
3 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
571 ms |
7536 KB |
Output is correct |
2 |
Execution timed out |
1004 ms |
13224 KB |
Time limit exceeded |
3 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
97 ms |
8900 KB |
Output is correct |
2 |
Execution timed out |
1080 ms |
16548 KB |
Time limit exceeded |
3 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
1036 ms |
20048 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |