QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#388876#7872. 崩坏天际线valeriu30 4292ms54088kbC++205.2kb2024-04-13 21:03:122024-04-13 21:03:13

Judging History

你现在查看的是最新测评结果

  • [2024-04-13 21:03:13]
  • 评测
  • 测评结果:30
  • 用时:4292ms
  • 内存:54088kb
  • [2024-04-13 21:03:12]
  • 提交

answer

#include <bits/stdc++.h>
#define all(x) (x).begin(),(x).end()
using namespace std;

using ll = long long;
using ld = long double;

//#define int ll
#define sz(x) ((int)(x).size())

using pii = pair<int,int>;
using tii = tuple<int,int,int>;

const int inf = 1e9 + 6, nmax = 5e4 + 5;

const int mod = 998244353;
struct Mint {
  int val;
  Mint(ll x = 0): val((x % mod + mod) % mod) {;}
  Mint operator +(const Mint& x) const { return Mint(val + x.val); }
  Mint operator -(const Mint& x) const { return Mint(val - x.val); }
  Mint operator *(const Mint& x) const { return Mint((ll)val * x.val); }
  Mint operator +=(const Mint& x) { return *this = Mint(val + x.val); }
  Mint operator -=(const Mint& x) { return *this = Mint(val - x.val); }
  Mint operator *=(const Mint& x) { return *this = Mint((ll)val * x.val); }
  Mint operator ^(const int& _b) const {
    Mint accum = 1, a = *this;
    if(val == 2 && _b == mod - 2) return (mod + 1) / 2;
    int b = _b;
    while(b) {
      accum = (b & 1? accum * a : accum);
      a *= a;
      b >>= 1;
    }
    return accum;
  }
  Mint operator /(const Mint& x) { return Mint((ll)val * (x ^ (mod - 2)).val); }
  Mint operator /=(const Mint& x) { return *this = Mint((ll)val * (x ^ (mod - 2)).val); }
};
Mint ip2[nmax], p2[nmax];


struct Segm { int l, r; Mint prob; };

vector<Segm> brutOver(vector<tii> oper) { // lumea nu este gata pentru a aprecia frumusetea acestui brut (gen, cand tineam mapuri pe bucati); bucatile se schimba conform timpului
                                          // tot le-as normaliza cu radix sau cv
   vector<Segm> open;
   for(auto [t, a, b] : oper) {
      if(t == 1)
         open.emplace_back(a, b, 1);
      else {
         const int X = a;         
         int D = sz(open);
         for(int i = 0; i < D; i++) {
            auto [l, r, P] = open[i];
            if(l < X && X < r) {
               P /= 2;
               open.emplace_back(X, r, P);
               open[i].r = X;
               open[i].prob /= 2;
            }
         }
      }
   }
   
   return open;
}

Mint extbatch(int cmax, vector<Segm> segm, vector<int> wcuts) {
   
   vector<int> cut(cmax + 1, 0);
   for(int i = 0; i < sz(wcuts); i++) {
      int x = wcuts[i];
      if(cut[x] == 0) cut[x] = i + 1;
   }
   
   vector<Mint> rez(sz(segm));
   vector<bool> fals(sz(segm), 1);
   int n = sz(segm);
   
   auto sweep = [&]() {
      vector<pair<int, Mint>> st;
      vector<Mint> spart;
      cut[0] = -inf;
      st.emplace_back(0, 0);
      spart.emplace_back(0);
      
      vector<vector<pii>> qs(cmax + 1, vector<pii>());
      int II = 0;
      for(auto [l, r, P] : segm)
         qs[r].emplace_back(l, II++);
      
      for(int i = 1; i <= cmax; i++) {
         bool revert = 0;
         if(cut[i] == 0)
            revert = 1;
         
         Mint sum = i - st.back().first;
         spart.emplace_back(spart.back() + sum * ip2[sz(st) - 1]);
         st.emplace_back(i, sum);
         
         for(auto [l, idx] : qs[i]) {
            auto it = distance(st.begin(), upper_bound(all(st), make_pair(l, Mint(0)), [&](pair<int, Mint> a, pair<int, Mint> b) { return a.first < b.first; })) - 1;
            
            if(it + 1 != spart.size() - 1) fals[idx] = 0;
            
            rez[idx] += (spart.back() - spart[it + 1]) * p2[it];
         }
         
         if(revert) {
            st.pop_back();
            spart.pop_back();
         }
         else {
            st.pop_back();
            spart.pop_back();
            while(cut[i] < cut[st.back().first]) {
               sum = sum / 2 + st.back().second / 2;
               st.pop_back();
               spart.pop_back();
            }
            spart.emplace_back(spart.back() + sum * ip2[sz(st)]);
            st.emplace_back(i, sum);
         }
      }
   };
   
   sweep();
   for(auto& [l, r, P] : segm) {
      r = cmax - r + 1;
      l = cmax - l + 1;
      swap(l, r);
   }
   reverse(cut.begin() + 1, cut.begin() + cmax + 1);
   sweep();
   
   Mint total;
   for(int i = 0; i < n; i++) {
      total += rez[i] * segm[i].prob;
      if(fals[i]) total += Mint(segm[i].r - segm[i].l) * segm[i].prob;
   }
   return total;
}

vector<int> rev(vector<int> D) {
   reverse(all(D));
   return D;
}

signed main() {
   ip2[0] = 1;
   p2[0] = 1;
   for(int i = 1; i < nmax; i++)
      ip2[i] = ip2[i - 1] * ((mod + 1) / 2),
      p2[i] = p2[i - 1] * 2;
      
   cin.tie(0) -> sync_with_stdio(0);
   int n, q;
   cin >> n >> q;
   
   vector<tii> oper(q);
   vector<int> cuts;
   for(auto &[t, l, r] : oper) {
      cin >> t >> l; 
      if(t == 1) cin >> r;
      else cuts.emplace_back(l);
   }
   reverse(all(cuts));
   
   auto R = brutOver(oper);
   
   Mint rez = 0;
   
   vector<tii> currst;
   const int B = 500;
   for(auto [t, l, r] : oper) {
      currst.emplace_back(t, l, r);
      if(t == 2) cuts.pop_back();
      if(sz(currst) >= B) {
         rez += extbatch(n, brutOver(currst), rev(cuts));
         currst.clear();
      }
   }
   
   rez += extbatch(n, brutOver(currst), rev(cuts));
   cout << rez.val << '\n';
}

/**
      Istenem! Nu poate fi real.
-- Surse verificate
*/


详细

Subtask #1:

score: 10
Accepted

Test #1:

score: 10
Accepted
time: 0ms
memory: 4460kb

input:

500 500
1 119 258
2 134
2 417
2 176
2 61
2 60
2 110
1 335 336
1 96 111
2 202
1 138 344
2 358
2 134
1 29 54
1 73 381
1 179 495
2 490
2 418
2 186
2 183
1 168 340
2 78
1 15 27
2 373
1 245 498
1 372 495
2 244
2 63
1 174 490
2 282
2 417
1 272 408
1 109 134
2 303
2 345
1 238 401
1 36 480
1 21 483
2 10
1 3...

output:

855279801

result:

ok single line: '855279801'

Test #2:

score: 0
Accepted
time: 3ms
memory: 4356kb

input:

495 466
1 35 393
2 236
1 4 335
2 455
1 36 470
1 23 61
2 195
2 109
2 451
1 282 491
2 238
2 117
2 468
1 2 60
1 439 487
2 238
1 209 294
2 321
2 309
1 113 183
2 409
2 87
2 130
2 124
2 176
2 448
2 379
1 181 446
2 146
2 450
1 171 423
2 355
2 332
1 123 387
1 151 269
1 17 417
2 122
1 324 494
1 265 439
2 225...

output:

294468977

result:

ok single line: '294468977'

Test #3:

score: 0
Accepted
time: 2ms
memory: 4248kb

input:

441 467
2 180
1 51 344
2 180
1 16 345
1 39 419
1 64 432
2 176
1 35 372
2 426
1 8 415
1 1 439
1 17 430
2 433
1 89 369
1 83 353
2 292
1 1 421
1 63 430
1 33 345
1 69 421
1 49 373
1 77 343
1 24 393
1 90 375
1 8 425
2 322
2 61
2 112
2 209
1 39 406
1 12 426
1 29 430
1 50 374
1 47 394
1 9 387
2 234
1 19 35...

output:

526117259

result:

ok single line: '526117259'

Test #4:

score: 0
Accepted
time: 5ms
memory: 4544kb

input:

500 500
2 442
1 12 414
1 40 435
2 138
1 79 448
1 16 464
2 163
1 94 492
2 97
2 335
1 7 452
1 25 474
1 78 442
2 286
1 93 430
1 78 438
2 469
2 354
2 270
2 292
2 108
2 301
1 100 480
2 258
1 17 487
2 2
2 409
2 385
2 338
1 83 454
1 41 490
1 95 475
1 43 442
1 66 445
2 406
2 168
1 10 406
2 330
2 20
1 90 491...

output:

810270061

result:

ok single line: '810270061'

Test #5:

score: 0
Accepted
time: 8ms
memory: 5280kb

input:

500 500
1 29 407
1 89 480
1 31 497
1 28 494
1 21 492
1 91 465
1 13 467
1 89 425
1 22 444
1 20 430
1 48 445
1 33 441
1 61 435
1 69 427
1 89 485
1 90 446
1 23 488
1 6 424
1 76 425
1 36 460
1 16 421
1 20 500
1 3 487
1 99 481
1 53 412
1 96 456
1 39 436
1 28 436
1 4 409
1 9 486
1 22 484
1 88 413
1 26 467...

output:

419428992

result:

ok single line: '419428992'

Test #6:

score: 0
Accepted
time: 12ms
memory: 5264kb

input:

500 500
1 85 442
1 20 473
1 10 441
1 31 426
1 95 478
1 60 454
1 54 491
1 97 464
1 14 443
1 88 474
1 28 462
1 97 410
1 99 496
1 96 493
1 62 479
1 12 466
1 64 471
1 43 490
1 50 411
1 85 448
1 48 433
1 30 456
1 39 462
1 46 409
1 63 494
1 39 409
1 36 436
1 27 463
1 37 498
1 69 464
1 8 441
1 99 436
1 84 ...

output:

519347055

result:

ok single line: '519347055'

Subtask #2:

score: 20
Accepted

Dependency #1:

100%
Accepted

Test #7:

score: 20
Accepted
time: 662ms
memory: 17016kb

input:

5000 5000
2 2254
2 4832
2 208
1 335 3080
2 481
1 527 3659
1 2645 3803
1 855 3544
2 3824
2 347
1 1567 4426
1 2184 4493
2 142
2 2451
1 995 4170
2 576
2 999
2 2726
1 278 3540
2 3218
1 922 3302
2 3253
2 4161
2 4505
1 4201 4534
1 1827 3540
2 3241
2 1909
2 2667
1 723 2453
2 3123
1 1017 4791
1 2953 3384
1 ...

output:

275175220

result:

ok single line: '275175220'

Test #8:

score: 0
Accepted
time: 552ms
memory: 17320kb

input:

4753 4704
1 589 2183
1 922 2210
2 2885
2 171
2 1597
2 3601
1 1906 4730
1 411 3615
2 1665
1 87 801
2 3525
2 2426
2 2723
1 323 4345
2 3950
2 460
2 4165
1 1156 2642
1 1490 3965
1 329 4081
1 1206 2077
2 4216
1 996 2254
2 2219
2 1035
2 4074
2 714
1 952 2726
2 3097
2 409
1 3320 4713
2 4061
1 1765 2040
1 2...

output:

840227126

result:

ok single line: '840227126'

Test #9:

score: 0
Accepted
time: 482ms
memory: 28988kb

input:

4141 4610
2 3761
2 2872
1 334 3247
1 273 3914
1 307 3651
1 607 4105
1 458 3269
1 270 3782
2 311
1 533 3332
2 2495
1 991 3573
1 376 3593
1 239 3682
1 259 3350
1 213 3380
2 1904
1 591 3512
1 845 3785
1 189 3335
1 817 3362
1 335 3288
2 3633
1 747 3586
2 4062
2 3812
1 487 3333
1 740 4002
1 847 3937
1 53...

output:

597472157

result:

ok single line: '597472157'

Test #10:

score: 0
Accepted
time: 1493ms
memory: 29380kb

input:

5000 5000
2 2864
1 473 4676
2 858
2 2672
2 4473
2 800
2 3259
2 470
2 3859
2 2228
1 491 4536
1 700 4378
2 498
1 769 4837
1 80 4861
1 109 4201
1 908 4094
1 9 4706
2 1017
2 737
2 4155
1 270 4290
2 4434
2 1867
1 148 4119
1 299 4194
2 4076
2 1863
2 1570
2 4855
1 1000 4834
1 637 4827
2 1961
2 4518
1 811 4...

output:

251906928

result:

ok single line: '251906928'

Test #11:

score: 0
Accepted
time: 4292ms
memory: 54088kb

input:

5000 5000
1 327 4388
1 768 4973
1 438 4243
1 288 4244
1 105 4460
1 862 4894
1 125 4611
1 934 4115
1 631 4349
1 635 4088
1 250 4629
1 873 4204
1 977 4296
1 391 4821
1 107 4589
1 86 4810
1 615 4072
1 221 4113
1 745 4771
1 806 4983
1 675 4334
1 709 4428
1 587 4180
1 494 4949
1 904 4253
1 901 4527
1 717...

output:

845230417

result:

ok single line: '845230417'

Test #12:

score: 0
Accepted
time: 3835ms
memory: 53668kb

input:

5000 5000
1 902 4097
1 263 4218
1 502 4305
1 798 4433
1 392 4689
1 479 4006
1 518 4269
1 764 4295
1 48 4834
1 966 4574
1 374 4970
1 950 4925
1 54 4860
1 987 4144
1 448 4504
1 329 4838
1 734 4807
1 403 4387
1 275 4396
1 731 4769
1 206 4348
1 282 4258
1 676 4956
1 274 4943
1 892 4146
1 337 4962
1 798 ...

output:

678724707

result:

ok single line: '678724707'

Subtask #3:

score: 0
Time Limit Exceeded

Test #13:

score: 0
Time Limit Exceeded

input:

50000 50000
1 24367 33007
1 14396 42256
1 6375 22327
1 7892 42501
1 10100 37998
1 6284 48524
1 7357 18164
1 16200 46424
1 18972 34131
1 16849 32591
1 1917 3018
1 19897 30272
1 45044 45753
1 18999 25448
1 5167 31033
1 6182 35335
1 7270 37270
1 12651 39965
1 28896 38022
1 13853 35426
1 35516 48244
1 1...

output:


result:


Subtask #4:

score: 0
Skipped

Dependency #1:

100%
Accepted

Dependency #2:

100%
Accepted

Dependency #3:

0%