QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#699179 | #7998. 矩阵 | hos_lyric | 100 ✓ | 2121ms | 75040kb | C++14 | 8.4kb | 2024-11-02 03:03:59 | 2024-11-02 03:04:05 |
Judging History
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"