QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#699179#7998. 矩阵hos_lyric100 ✓2121ms75040kbC++148.4kb2024-11-02 03:03:592024-11-02 03:04:05

Judging History

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

  • [2024-11-02 03:04:05]
  • 评测
  • 测评结果:100
  • 用时:2121ms
  • 内存:75040kb
  • [2024-11-02 03:03:59]
  • 提交

answer

#include <cassert>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iostream>
#include <limits>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>

using namespace std;

using Int = long long;

template <class T1, class T2> ostream &operator<<(ostream &os, const pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> ostream &operator<<(ostream &os, const vector<T> &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; }
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; }
template <class T> bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; }
template <class T> bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; }
#define COLOR(s) ("\x1b[" s "m")

////////////////////////////////////////////////////////////////////////////////
template <unsigned M_> struct ModInt {
  static constexpr unsigned M = M_;
  unsigned x;
  constexpr ModInt() : x(0U) {}
  constexpr ModInt(unsigned x_) : x(x_ % M) {}
  constexpr ModInt(unsigned long long x_) : x(x_ % M) {}
  constexpr ModInt(int x_) : x(((x_ %= static_cast<int>(M)) < 0) ? (x_ + static_cast<int>(M)) : x_) {}
  constexpr ModInt(long long x_) : x(((x_ %= static_cast<long long>(M)) < 0) ? (x_ + static_cast<long long>(M)) : x_) {}
  ModInt &operator+=(const ModInt &a) { x = ((x += a.x) >= M) ? (x - M) : x; return *this; }
  ModInt &operator-=(const ModInt &a) { x = ((x -= a.x) >= M) ? (x + M) : x; return *this; }
  ModInt &operator*=(const ModInt &a) { x = (static_cast<unsigned long long>(x) * a.x) % M; return *this; }
  ModInt &operator/=(const ModInt &a) { return (*this *= a.inv()); }
  ModInt pow(long long e) const {
    if (e < 0) return inv().pow(-e);
    ModInt a = *this, b = 1U; for (; e; e >>= 1) { if (e & 1) b *= a; a *= a; } return b;
  }
  ModInt inv() const {
    unsigned a = M, b = x; int y = 0, z = 1;
    for (; b; ) { const unsigned q = a / b; const unsigned c = a - q * b; a = b; b = c; const int w = y - static_cast<int>(q) * z; y = z; z = w; }
    assert(a == 1U); return ModInt(y);
  }
  ModInt operator+() const { return *this; }
  ModInt operator-() const { ModInt a; a.x = x ? (M - x) : 0U; return a; }
  ModInt operator+(const ModInt &a) const { return (ModInt(*this) += a); }
  ModInt operator-(const ModInt &a) const { return (ModInt(*this) -= a); }
  ModInt operator*(const ModInt &a) const { return (ModInt(*this) *= a); }
  ModInt operator/(const ModInt &a) const { return (ModInt(*this) /= a); }
  template <class T> friend ModInt operator+(T a, const ModInt &b) { return (ModInt(a) += b); }
  template <class T> friend ModInt operator-(T a, const ModInt &b) { return (ModInt(a) -= b); }
  template <class T> friend ModInt operator*(T a, const ModInt &b) { return (ModInt(a) *= b); }
  template <class T> friend ModInt operator/(T a, const ModInt &b) { return (ModInt(a) /= b); }
  explicit operator bool() const { return x; }
  bool operator==(const ModInt &a) const { return (x == a.x); }
  bool operator!=(const ModInt &a) const { return (x != a.x); }
  friend std::ostream &operator<<(std::ostream &os, const ModInt &a) { return os << a.x; }
};
////////////////////////////////////////////////////////////////////////////////

using Mint998 = ModInt<998244353>;
using Mint = ModInt<1000000007>;


int N, Q;
Mint A[3010][3010];
Mint B[3010][3010];

struct Rect {
  int x0, x1, y0, y1;
  int ox, oy;
  int rot;
  Mint lz;
  friend ostream &operator<<(ostream &os, const Rect &r) {
    return os << "Rect([" << r.x0 << "," << r.x1 << ")*[" << r.y0 << "," << r.y1 << "); " << r.ox << "," << r.oy << "; " << r.rot << "; " << r.lz << ")";
  }
  bool any() const {
    return (x0 < x1 && y0 < y1);
  }
};
vector<Rect> rs;

Rect tmp[5];
bool cut(const Rect &r, int x0, int y0, int x1, int y1) {
  if (x1 <= r.x0 || r.x1 <= x0) return false;
  if (y1 <= r.y0 || r.y1 <= y0) return false;
  chmax(x0, r.x0); chmin(x1, r.x1);
  chmax(y0, r.y0); chmin(y1, r.y1);
  fill(tmp, tmp + 5, r);
  tmp[0].x1 = x0;
  tmp[1].x0 = x1;
  tmp[2].x0 = x0; tmp[2].x1 = x1; tmp[2].y1 = y0;
  tmp[3].x0 = x0; tmp[3].x1 = x1; tmp[3].y0 = y1;
  tmp[4].x0 = x0; tmp[4].x1 = x1; tmp[4].y0 = y0; tmp[4].y1 = y1;
  if (r.rot == 0) {
    tmp[1].ox += (x1 - r.x0);
    tmp[2].ox += (x0 - r.x0);
    tmp[3].ox += (x0 - r.x0); tmp[3].oy += (y1 - r.y0);
    tmp[4].ox += (x0 - r.x0); tmp[4].oy += (y0 - r.y0);
  } else if (r.rot == 1) {
    tmp[0].oy += (r.x1 - x0);
    tmp[2].oy += (r.x1 - x1);
    tmp[3].oy += (r.x1 - x1); tmp[3].ox += (y1 - r.y0);
    tmp[4].oy += (r.x1 - x1); tmp[4].ox += (y0 - r.y0);
  } else if (r.rot == 2) {
    tmp[0].ox += (r.x1 - x0);
    tmp[2].ox += (r.x1 - x1); tmp[2].oy += (r.y1 - y0);
    tmp[3].ox += (r.x1 - x1); 
    tmp[4].ox += (r.x1 - x1); tmp[4].oy += (r.y1 - y1);
  } else if (r.rot == 3) {
    tmp[1].oy += (x1 - r.x0);
    tmp[2].oy += (x0 - r.x0); tmp[2].ox += (r.y1 - y0);
    tmp[3].oy += (x0 - r.x0);
    tmp[4].oy += (x0 - r.x0); tmp[4].ox += (r.y1 - y1);
  } else {
    assert(false);
  }
  return true;
}

void build() {
  rs = {Rect{0, N, 0, N, 0, 0, 0, Mint(0)}};
}
void dliub() {
  for (const Rect &r : rs) {
    const int m = r.x1 - r.x0;
    const int n = r.y1 - r.y0;
    for (int dx = 0; dx < m; ++dx) for (int dy = 0; dy < n; ++dy) {
      int cx, cy;
      if (r.rot == 0) {
        cx = dx;
        cy = dy;
      } else if (r.rot == 1) {
        cx = dy;
        cy = m - 1 - dx;
      } else if (r.rot == 2) {
        cx = m - 1 - dx;
        cy = n - 1 - dy;
      } else if (r.rot == 3) {
        cx = n - 1 - dy;
        cy = dx;
      } else {
        assert(false);
      }
      B[r.x0 + dx][r.y0 + dy] = A[r.ox + cx][r.oy + cy] + r.lz;
    }
  }
// cerr<<COLOR("96");for(int x=0;x<N;++x)pv(B[x],B[x]+N);cerr<<COLOR();
  for (int x = 0; x < N; ++x) for (int y = 0; y < N; ++y) A[x][y] = B[x][y];
}

int main() {
  for (; ~scanf("%d%d", &N, &Q); ) {
    for (int x = 0; x < N; ++x) {
      Mint998 v = 1;
      for (int y = 0; y < N; ++y) {
        v *= (x + 2);
        A[x][y] = v.x;
      }
    }
    build();
    
    const int BLOCK = sqrt(N);
    for (int q = 0; q < Q; ++q) {
      int O, X0, X1, Y0, Y1;
      scanf("%d%d%d%d%d", &O, &X0, &Y0, &X1, &Y1);
      --X0;
      --Y0;
      if (O == 1) {
        const int D = X1 - X0;
        assert(D == Y1 - Y0);
// cerr<<COLOR("93")<<"["<<X0<<", "<<X1<<") * ["<<Y0<<", "<<Y1<<")"<<COLOR()<<endl;
        vector<Rect> rrs;
        for (const auto &r : rs) {
          if (cut(r, X0, Y0, X1, Y1)) {
            const int x0 = tmp[4].x0;
            const int x1 = tmp[4].x1;
            const int y0 = tmp[4].y0;
            const int y1 = tmp[4].y1;
            tmp[4].x0 = X0 + (D - (y1 - Y0));
            tmp[4].x1 = X0 + (D - (y0 - Y0));
            tmp[4].y0 = Y0 + (x0 - X0);
            tmp[4].y1 = Y0 + (x1 - X0);
            ++tmp[4].rot &= 3;
            for (int i = 0; i < 5; ++i) if (tmp[i].any()) {
              rrs.push_back(tmp[i]);
            }
          } else {
            rrs.push_back(r);
          }
        }
        rs.swap(rrs);
      } else if (O == 2) {
        Mint D;
        scanf("%u", &D.x);
// cerr<<COLOR("94")<<"["<<X0<<", "<<X1<<") * ["<<Y0<<", "<<Y1<<") "<<D<<COLOR()<<endl;
        vector<Rect> rrs;
        for (const auto &r : rs) {
          if (cut(r, X0, Y0, X1, Y1)) {
            tmp[4].lz += D;
            for (int i = 0; i < 5; ++i) if (tmp[i].any()) {
              rrs.push_back(tmp[i]);
            }
          } else {
            rrs.push_back(r);
          }
        }
        rs.swap(rrs);
      } else {
        assert(false);
      }
// cerr<<"rs = "<<rs<<endl;
      if ((q + 1) % BLOCK == 0) {
        dliub();
        build();
      }
    }
    dliub();
    
    Mint key = 0;
    for (int x = N; --x >= 0; ) for (int y = N; --y >= 0; ) {
      (key += A[x][y].x) *= 12345;
    }
    printf("%u\n", key.x);
  }
  return 0;
}

详细


Pretests


Final Tests

Test #1:

score: 10
Accepted
time: 7ms
memory: 8432kb

input:

100 3000
1 4 3 100 99
1 3 8 83 88
1 4 7 96 99
1 2 7 90 95
2 6 5 98 98 762240996
1 9 8 98 97
1 3 3 100 100
2 18 7 97 84 230446266
1 2 5 94 97
1 2 18 82 98
1 9 2 100 93
1 3 9 91 97
1 1 2 99 100
1 7 6 98 97
1 2 4 98 100
1 3 2 99 98
1 4 2 100 98
2 9 74 28 83 703323901
1 6 4 100 98
2 2 6 99 93 882702932
...

output:

535773693

result:

ok 1 number(s): "535773693"

Test #2:

score: 10
Accepted
time: 1630ms
memory: 74828kb

input:

3000 3000
2 40 215 2496 2298 67590963
2 100 211 2684 2855 670656770
2 417 399 2612 2930 826127245
2 116 779 2407 2596 977282152
2 512 366 2768 2949 737671139
2 191 594 2671 2978 404448620
2 3 11 2720 2633 516813112
2 95 294 2553 2820 112269230
2 284 58 2934 2754 225336228
2 306 53 2815 2984 49218367...

output:

343427775

result:

ok 1 number(s): "343427775"

Test #3:

score: 10
Accepted
time: 1435ms
memory: 74744kb

input:

3000 2000
1 853 610 2674 2431
1 24 415 2490 2881
1 19 150 2842 2973
1 597 79 2977 2459
1 307 408 2397 2498
1 292 316 2907 2931
1 7 42 2953 2988
1 36 30 2953 2947
1 55 438 2609 2992
1 97 204 2664 2771
1 607 254 2803 2450
1 59 159 2808 2908
1 58 12 2991 2945
1 287 400 2506 2619
1 105 143 2949 2987
1 2...

output:

743809971

result:

ok 1 number(s): "743809971"

Test #4:

score: 10
Accepted
time: 1432ms
memory: 74844kb

input:

3000 2000
1 646 55 2815 2224
1 47 78 2964 2995
1 48 7 2882 2841
1 9 8 2999 2998
1 128 36 2984 2892
1 14 5 2990 2981
1 52 81 2934 2963
1 231 317 2612 2698
1 33 7 2960 2934
1 111 165 2732 2786
1 9 2 2994 2987
1 131 224 2866 2959
1 19 49 2917 2947
1 46 44 2997 2995
1 1 14 2969 2982
1 69 49 2975 2955
1 ...

output:

517390609

result:

ok 1 number(s): "517390609"

Test #5:

score: 10
Accepted
time: 2121ms
memory: 74744kb

input:

3000 3000
1 6 5 2997 2996
1 14 362 2624 2972
1 657 108 2405 1856
1 20 27 2983 2990
1 245 236 2903 2894
1 184 100 2957 2873
1 61 47 2954 2940
1 39 89 2935 2985
1 43 49 2956 2962
1 132 138 2763 2769
1 4 4 2998 2998
1 16 59 2752 2795
1 30 24 2997 2991
1 143 11 2959 2827
1 57 44 2887 2874
1 40 17 2991 2...

output:

632994577

result:

ok 1 number(s): "632994577"

Test #6:

score: 10
Accepted
time: 2077ms
memory: 74792kb

input:

3000 3000
1 544 5 2907 2368
1 60 63 2978 2981
1 228 310 2881 2963
1 57 71 2958 2972
1 101 565 2467 2931
1 979 1004 2882 2907
1 37 3 2913 2879
1 418 427 2922 2931
1 30 62 2937 2969
1 83 104 2979 3000
1 4 18 2939 2953
1 67 203 2834 2970
1 357 882 1000 1525
1 34 10 2991 2967
1 8 4 2997 2993
1 2 4 2991 ...

output:

275056325

result:

ok 1 number(s): "275056325"

Test #7:

score: 10
Accepted
time: 1440ms
memory: 75040kb

input:

3000 2000
1 82 42 2883 2843
1 2 48 2939 2985
2 188 176 2376 2976 715286463
1 15 24 2984 2993
1 106 50 2811 2755
2 45 132 2745 2432 339749026
1 5 100 2779 2874
1 343 134 2478 2269
1 106 227 2820 2941
2 10 206 2947 2798 365598855
2 560 106 2582 2903 92795190
1 74 100 2960 2986
1 5 33 2944 2972
2 873 3...

output:

191562107

result:

ok 1 number(s): "191562107"

Test #8:

score: 10
Accepted
time: 1451ms
memory: 74816kb

input:

3000 2000
1 20 5 2990 2975
1 14 10 2999 2995
1 192 89 2950 2847
2 547 646 2654 2940 439095457
1 18 16 3000 2998
1 35 105 2754 2824
1 17 52 2942 2977
2 356 121 2384 2963 959830720
2 76 78 2412 1842 361763869
2 199 46 2624 2931 743439747
2 180 320 2981 2964 772483137
2 173 122 2732 2901 759079379
1 47...

output:

787136396

result:

ok 1 number(s): "787136396"

Test #9:

score: 10
Accepted
time: 2063ms
memory: 74756kb

input:

3000 3000
1 220 131 2785 2696
2 1443 1396 2159 2850 578335729
1 536 1184 1545 2193
1 26 23 2999 2996
1 43 56 2881 2894
1 29 4 2872 2847
2 69 155 2681 2931 88398671
1 66 44 2945 2923
1 53 512 2515 2974
2 402 58 2752 2782 641047796
2 458 36 2859 2756 694998888
2 145 192 2995 2973 483103477
1 117 315 2...

output:

973923858

result:

ok 1 number(s): "973923858"

Test #10:

score: 10
Accepted
time: 2086ms
memory: 74956kb

input:

3000 3000
2 7 47 2887 2544 535526356
2 160 639 2988 2995 343766215
1 122 40 2985 2903
2 62 240 2858 2177 841981272
1 169 18 2682 2531
1 61 286 2729 2954
1 180 27 2994 2841
1 238 330 2881 2973
2 57 188 2880 2914 203515190
1 28 61 2903 2936
2 102 269 2876 2397 681849386
1 3 53 2896 2946
1 101 56 2985 ...

output:

328382756

result:

ok 1 number(s): "328382756"