/** Im the best because i work as hard as i possibly can **/
#pragma GCC optimize("Ofast")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx")
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
#define all(x) (x).begin(),(x).end()
#define F first
#define S second
#define Mp make_pair
#define fast_io ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define file_io freopen("in.txt" , "r+" , stdin) ; freopen("out.txt" , "w+" , stdout);
#define endl "\n"
#define SZ(x) (int)x.size()
const int N = 1e5 + 10;
const ll mod = 1e9 + 7;
const ll mod2 = 998244353;
const ll inf = 8e18;
const int LOG = 22;
const int M = 52;
char C[N];
int n, m, A[N][M], I[N][M], par[N * M], arr[N][M], par2[N * M];
pii RI[N * M];
ll tot, ps[N], cnt[N];
int get(int x) { return ((x == par[x])? x : par[x] = get(par[x])); }
inline int ok(int i)
{
return A[RI[i].F][RI[i].S];
}
int unify(int v, int u)
{
if(!ok(v) || !ok(u)) return 0;
if(u < v) swap(u, v);
v = get(v), u = get(u);
if(v == u) return 0;
par[u] = v;
return 1;
}
int get2(int x) { return ((x == par2[x])? x : par2[x] = get2(par2[x])); }
int unify2(int v, int u)
{
if(!ok(v) || !ok(u)) return 0;
v = get2(v), u = get2(u);
if(v == u) return 0;
par2[u] = v;
return 1;
}
void solve(int l, int r)
{
if(l > r) return;
if(l == r)
{
for(int j = 1; j <= m; j ++)
{
tot += A[l][j];
tot -= (A[l][j] + A[l][j - 1] == 2);
}
return;
}
int mid = (l + r) >> 1;
solve(l, mid); solve(mid + 1, r);
cnt[l - 1] = cnt[r + 1] = ps[l - 1] = ps[r + 1] = 0;
for(int i = l; i <= r; i ++)
{
ps[i] = cnt[i] = 0;
for(int j = 1; j <= m; j ++)
{
par[I[i][j]] = I[i][j];
arr[i][j] = 0;
}
}
vector < ll > merges;
merges.push_back(mid + 1);
for(int i = mid + 1; i <= r; i ++)
{
if(i > mid + 1) cnt[i] = cnt[i - 1];
for(int j = 1; j <= m; j ++)
{
cnt[i] += A[i][j];
cnt[i] -= unify(I[i][j], I[i][j - 1]);
}
if(i > mid + 1)
{
for(int j = 1; j <= m; j ++)
{
cnt[i] -= unify(I[i][j], I[i - 1][j]);
}
}
for(int j = 1; j <= m; j ++) arr[i][j] = get(I[mid + 1][j]);
if(i > mid + 1)
{
int flag = 0;
for(int j = 1; j <= m; j ++)
{
if(A[mid + 1][j] == 0) continue;
if(arr[i][j] != arr[i - 1][j]) flag = 1;
}
if(flag) merges.push_back(i);
}
ps[i] = cnt[i] + (i > mid + 1? ps[i - 1] : 0);
}
merges.push_back(r + 1);
ll cu = 0;
for(int i = mid; i >= l; i --)
{
for(int j = 1; j <= m; j ++)
{
cu += A[i][j];
cu -= unify(I[i][j], I[i][j - 1]);
}
if(i < mid)
{
for(int j = 1; j <= m; j ++)
{
cu -= unify(I[i][j], I[i + 1][j]);
}
}
for(int id = 0; id < SZ(merges) - 1; id ++)
{
for(int j = 1; j <= m; j ++)
{
par2[get(I[mid][j])] = get(I[mid][j]);
par2[arr[merges[id]][j]] = arr[merges[id]][j];
}
ll yal = 0;
for(int j = 1; j <= m; j ++)
{
if((A[mid][j] + A[mid + 1][j]) < 2) continue;
ll v = get(I[mid][j]), u = arr[merges[id]][j];
yal += unify2(v, u);
}
tot += 1ll * (cu - yal) * (merges[id + 1] - merges[id]) + (ps[merges[id + 1] - 1] - ps[merges[id] - 1]);
}
}
}
int main()
{
scanf("%d%d", &n, &m);
RI[0] = Mp(0, 0);
for(int i = 1; i <= n; i ++)
{
for(int j = 1; j <= m; j ++)
{
I[i][j] = (i - 1) * m + j;
RI[I[i][j]] = Mp(i, j);
}
}
for(int i = 1; i <= n; i ++)
{
scanf("%s", C);
for(int j = 0; j < m; j ++)
{
A[i][j + 1] = (C[j] == '1');
}
}
solve(1, n);
printf("%lld\n", tot);
return 0;
}
Compilation message
raspad.cpp: In function 'int main()':
raspad.cpp:156:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
156 | scanf("%d%d", &n, &m);
| ~~~~~^~~~~~~~~~~~~~~~
raspad.cpp:168:8: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
168 | scanf("%s", C);
| ~~~~~^~~~~~~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
492 KB |
Output is correct |
2 |
Correct |
5 ms |
492 KB |
Output is correct |
3 |
Correct |
2 ms |
492 KB |
Output is correct |
4 |
Correct |
2 ms |
512 KB |
Output is correct |
5 |
Correct |
3 ms |
492 KB |
Output is correct |
6 |
Correct |
4 ms |
492 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
492 KB |
Output is correct |
2 |
Correct |
5 ms |
492 KB |
Output is correct |
3 |
Correct |
2 ms |
492 KB |
Output is correct |
4 |
Correct |
2 ms |
512 KB |
Output is correct |
5 |
Correct |
3 ms |
492 KB |
Output is correct |
6 |
Correct |
4 ms |
492 KB |
Output is correct |
7 |
Correct |
19 ms |
1644 KB |
Output is correct |
8 |
Correct |
1 ms |
748 KB |
Output is correct |
9 |
Correct |
106 ms |
1772 KB |
Output is correct |
10 |
Correct |
17 ms |
1516 KB |
Output is correct |
11 |
Correct |
36 ms |
1772 KB |
Output is correct |
12 |
Correct |
33 ms |
1516 KB |
Output is correct |
13 |
Correct |
44 ms |
1644 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
673 ms |
42636 KB |
Output is correct |
2 |
Correct |
1245 ms |
86532 KB |
Output is correct |
3 |
Correct |
3618 ms |
86640 KB |
Output is correct |
4 |
Correct |
670 ms |
77916 KB |
Output is correct |
5 |
Correct |
326 ms |
26348 KB |
Output is correct |
6 |
Correct |
1442 ms |
86508 KB |
Output is correct |
7 |
Correct |
1749 ms |
86480 KB |
Output is correct |
8 |
Correct |
1080 ms |
69428 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
492 KB |
Output is correct |
2 |
Correct |
5 ms |
492 KB |
Output is correct |
3 |
Correct |
2 ms |
492 KB |
Output is correct |
4 |
Correct |
2 ms |
512 KB |
Output is correct |
5 |
Correct |
3 ms |
492 KB |
Output is correct |
6 |
Correct |
4 ms |
492 KB |
Output is correct |
7 |
Correct |
19 ms |
1644 KB |
Output is correct |
8 |
Correct |
1 ms |
748 KB |
Output is correct |
9 |
Correct |
106 ms |
1772 KB |
Output is correct |
10 |
Correct |
17 ms |
1516 KB |
Output is correct |
11 |
Correct |
36 ms |
1772 KB |
Output is correct |
12 |
Correct |
33 ms |
1516 KB |
Output is correct |
13 |
Correct |
44 ms |
1644 KB |
Output is correct |
14 |
Correct |
673 ms |
42636 KB |
Output is correct |
15 |
Correct |
1245 ms |
86532 KB |
Output is correct |
16 |
Correct |
3618 ms |
86640 KB |
Output is correct |
17 |
Correct |
670 ms |
77916 KB |
Output is correct |
18 |
Correct |
326 ms |
26348 KB |
Output is correct |
19 |
Correct |
1442 ms |
86508 KB |
Output is correct |
20 |
Correct |
1749 ms |
86480 KB |
Output is correct |
21 |
Correct |
1080 ms |
69428 KB |
Output is correct |
22 |
Correct |
3247 ms |
100732 KB |
Output is correct |
23 |
Execution timed out |
6096 ms |
87876 KB |
Time limit exceeded |
24 |
Halted |
0 ms |
0 KB |
- |