QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#765121 | #9479. And DNA | Fido_Puppy | AC ✓ | 1ms | 4020kb | C++23 | 4.5kb | 2024-11-20 12:19:00 | 2024-11-20 12:19:01 |
Judging History
answer
#include <bits/stdc++.h>
#define all(x) x.begin(), x.end()
#define pb push_back
#define eb emplace_back
#define rz resize
#define Siz(a) (int)(a.size())
#define MP make_pair
#define MT make_tuple
#define IT iterator
#define fi first
#define se second
#define For(i, a, b) for (int i = (int)(a); i <= (int)(b); ++i)
#define Rep(i, a, b) for (int i = (int)(a); i >= (int)(b); --i)
#define CLR(a, v) memset(a, v, sizeof(a))
#define CPY(a, b) memcpy(a, b, sizeof(a))
#define debug cerr << "ztxakking\n"
#define y0 ztxaknoi
#define y1 ztxakioi
using namespace std;
using ll = long long;
using ld = long double;
using ull = unsigned long long;
using uint = unsigned int;
using i128 = __int128_t;
using u128 = __uint128_t;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
using pli = pair<ll, int>;
using pil = pair<int, ll>;
using vi = vector<int>;
template<typename T>
using V = vector<T>;
mt19937 rnd(chrono::steady_clock::now().time_since_epoch().count());
template<const int MOD>
struct Modular {
int val;
Modular(int x = 0) : val(x) {}
friend Modular operator + (Modular x, Modular y) { int res = x.val + y.val; if (res >= MOD) res -= MOD; return res; }
friend Modular operator - (Modular x, Modular y) { int res = x.val - y.val; if (res < 0) res += MOD; return res; }
friend Modular operator * (Modular x, Modular y) { int res = 1ll * x.val * y.val % MOD; return res; }
Modular &operator += (Modular x) { *this = *this + x; return *this; }
Modular &operator -= (Modular x) { *this = *this - x; return *this; }
Modular &operator *= (Modular x) { *this = *this * x; return *this; }
Modular &operator ++ () { *this = *this + 1; return *this; }
Modular &operator ++ (int) { *this = *this + 1; return *this; }
Modular &operator -- () { *this = *this - 1; return *this; }
Modular &operator -- (int) { *this = *this - 1; return *this; }
Modular operator - () { int res = val; if (res) res = MOD - res; return res; }
static Modular pow(Modular x, ll p) {
Modular ans(1);
for (; p; p /= 2, x *= x) {
if (p & 1) ans *= x;
}
return ans;
}
static Modular inv(Modular x) { return pow(x, MOD - 2); }
friend Modular operator / (Modular x, Modular y) { return x * inv(y); }
Modular &operator /= (Modular x) { *this = *this / x; return *this; }
friend bool operator == (Modular x, Modular y) { return x.val == y.val; }
friend bool operator != (Modular x, Modular y) { return x.val != y.val; }
friend bool operator < (Modular x, Modular y) { return x.val < y.val; }
friend bool operator > (Modular x, Modular y) { return x.val > y.val; }
friend bool operator <= (Modular x, Modular y) { return x.val <= y.val; }
friend bool operator >= (Modular x, Modular y) { return x.val >= y.val; }
friend istream &operator >> (istream &in, Modular &item) { in >> item.val; return in; }
friend ostream &operator << (ostream &out, Modular item) { out << item.val; return out; }
};
using mint = Modular<998244353>;
int n, m;
mint w[3], f[31][4];
struct Matrix {
int n, m;
mint a[3][3];
Matrix() {
n = m = 0;
For(i, 0, 2) For(j, 0, 2) a[i][j] = 0;
}
friend Matrix operator * (Matrix A, Matrix B) {
Matrix C; C.n = A.n, C.m = B.m;
For(i, 0, C.n - 1) For(k, 0, A.m - 1) For(j, 0, C.m - 1) {
C.a[i][j] += A.a[i][k] * B.a[k][j];
}
return C;
}
} X, T;
Matrix qpow(Matrix X, Matrix T, int p) {
while (p) {
if (p & 1) X = T * X;
T = T * T, p >>= 1;
}
return X;
}
void Main() {
cin >> n >> m;
w[0] = w[2] = 1;
T.n = T.m = 3;
T.a[1][0] = T.a[2][0] = T.a[2][1] = T.a[0][2] = 1;
X.n = 3, X.m = 1;
X.a[0][0] = X.a[1][0] = X.a[2][0] = 0, X.a[0][0] = 1;
X = qpow(X, T, n - 2);
w[1] += X.a[0][0] + X.a[1][0];
X.a[0][0] = X.a[1][0] = X.a[2][0] = 0, X.a[1][0] = 1;
X = qpow(X, T, n - 2);
w[1] += X.a[2][0];
X.a[0][0] = X.a[1][0] = X.a[2][0] = 0, X.a[2][0] = 1;
X = qpow(X, T, n - 2);
w[1] += X.a[0][0] + X.a[2][0];
f[30][0] = 1;
Rep(i, 29, 0) For(j, 0, 3) {
int t = j << 1;
if (m >> i & 1) ++t;
For(k, 0, 2) {
if (t - k >= 0 && t - k <= 3) {
f[i][t - k] += f[i + 1][j] * w[k];
}
}
}
cout << f[0][0] << '\n';
}
int main() {
#ifndef ONLINE_JUDGE
freopen("C.in", "r", stdin) && freopen("C.out", "w", stdout);
#endif
ios::sync_with_stdio(0), cin.tie(0);
int t = 1;
// cin >> t;
while (t--) Main();
cerr << (double)(clock()) / CLOCKS_PER_SEC << '\n';
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 1ms
memory: 3992kb
input:
3 2
output:
4
result:
ok "4"
Test #2:
score: 0
Accepted
time: 0ms
memory: 3984kb
input:
3 0
output:
1
result:
ok "1"
Test #3:
score: 0
Accepted
time: 0ms
memory: 4020kb
input:
100 100
output:
343406454
result:
ok "343406454"
Test #4:
score: 0
Accepted
time: 0ms
memory: 3924kb
input:
686579306 119540831
output:
260602195
result:
ok "260602195"
Test #5:
score: 0
Accepted
time: 0ms
memory: 3968kb
input:
26855095 796233790
output:
492736984
result:
ok "492736984"
Test #6:
score: 0
Accepted
time: 0ms
memory: 3784kb
input:
295310488 262950628
output:
503008241
result:
ok "503008241"
Test #7:
score: 0
Accepted
time: 1ms
memory: 3912kb
input:
239670714 149827706
output:
994702969
result:
ok "994702969"
Test #8:
score: 0
Accepted
time: 0ms
memory: 3980kb
input:
790779949 110053353
output:
898252188
result:
ok "898252188"
Test #9:
score: 0
Accepted
time: 0ms
memory: 3992kb
input:
726600542 795285932
output:
183726777
result:
ok "183726777"
Test #10:
score: 0
Accepted
time: 0ms
memory: 3908kb
input:
957970519 585582861
output:
298814018
result:
ok "298814018"
Test #11:
score: 0
Accepted
time: 0ms
memory: 3980kb
input:
93349859 634036506
output:
110693443
result:
ok "110693443"
Test #12:
score: 0
Accepted
time: 0ms
memory: 3980kb
input:
453035113 34126396
output:
306244220
result:
ok "306244220"
Test #13:
score: 0
Accepted
time: 1ms
memory: 3980kb
input:
31994526 100604502
output:
930620701
result:
ok "930620701"
Test #14:
score: 0
Accepted
time: 0ms
memory: 3984kb
input:
234760741 249817734
output:
440195858
result:
ok "440195858"
Test #15:
score: 0
Accepted
time: 0ms
memory: 3980kb
input:
542621111 646412689
output:
207377992
result:
ok "207377992"
Test #16:
score: 0
Accepted
time: 0ms
memory: 3928kb
input:
28492783 602632297
output:
65234506
result:
ok "65234506"
Test #17:
score: 0
Accepted
time: 0ms
memory: 3928kb
input:
213500301 768820204
output:
540205113
result:
ok "540205113"
Test #18:
score: 0
Accepted
time: 0ms
memory: 3980kb
input:
697808101 753041955
output:
93561295
result:
ok "93561295"
Test #19:
score: 0
Accepted
time: 1ms
memory: 4016kb
input:
585126464 450455977
output:
91109717
result:
ok "91109717"
Test #20:
score: 0
Accepted
time: 1ms
memory: 3972kb
input:
236696315 482334538
output:
925575199
result:
ok "925575199"
Test #21:
score: 0
Accepted
time: 0ms
memory: 3984kb
input:
632719214 298704996
output:
358926097
result:
ok "358926097"
Test #22:
score: 0
Accepted
time: 0ms
memory: 3848kb
input:
869119333 933404114
output:
318501108
result:
ok "318501108"
Test #23:
score: 0
Accepted
time: 0ms
memory: 3744kb
input:
6977994 814763202
output:
585332987
result:
ok "585332987"
Test #24:
score: 0
Accepted
time: 0ms
memory: 3968kb
input:
3 1
output:
3
result:
ok "3"
Test #25:
score: 0
Accepted
time: 0ms
memory: 3984kb
input:
3 1000000000
output:
279679226
result:
ok "279679226"
Test #26:
score: 0
Accepted
time: 0ms
memory: 3976kb
input:
4 0
output:
1
result:
ok "1"
Test #27:
score: 0
Accepted
time: 0ms
memory: 3908kb
input:
4 1
output:
2
result:
ok "2"
Test #28:
score: 0
Accepted
time: 0ms
memory: 3980kb
input:
4 1000000000
output:
1755648
result:
ok "1755648"
Test #29:
score: 0
Accepted
time: 0ms
memory: 3992kb
input:
1000000000 0
output:
1
result:
ok "1"
Test #30:
score: 0
Accepted
time: 0ms
memory: 3908kb
input:
1000000000 1
output:
696798313
result:
ok "696798313"
Test #31:
score: 0
Accepted
time: 0ms
memory: 3984kb
input:
1000000000 1000000000
output:
703786301
result:
ok "703786301"
Extra Test:
score: 0
Extra Test Passed