Submission #370134

#TimeUsernameProblemLanguageResultExecution timeMemory
370134AriaHRaspad (COI17_raspad)C++11
61 / 100
6096 ms100732 KiB
/** 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 (stderr)

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);
      |   ~~~~~^~~~~~~~~
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...