QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#767518 | #9479. And DNA | Add_Catalyst | AC ✓ | 0ms | 3680kb | C++14 | 3.9kb | 2024-11-20 21:10:44 | 2024-11-20 21:10:54 |
Judging History
answer
#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define ll long long
#define RCL(a,b,c,d) memset(a,b,sizeof(c)*(d))
#define tomin(a,...) ((a)=min({(a),__VA_ARGS__}))
#define tomax(a,...) ((a)=max({(a),__VA_ARGS__}))
#define FOR(i,a,b) for(int i(a); i<=(int)(b); ++i)
#define DOR(i,a,b) for(int i(a); i>=(int)(b); --i)
#define EDGE(g,i,x,y) for(int i(g.h[x]),y(g[i].v); ~i; y=g[i=g[i].nxt].v)
using namespace std;
constexpr int lN(30),lV(lN+2);
namespace IOEstream {
#define isdigit(ch) ('0'<=(ch)&&(ch)<='9')
#define DE(...) E(#__VA_ARGS__,__VA_ARGS__)
struct Istream {
char getc() {
return getchar();
}
template<class T>void operator ()(T &x) {
static bool sign(0);
static char ch(0);
sign=0,x=0;
while(ch=getc(),!isdigit(ch))if(ch=='-')sign=1;
do x=(x<<3)+(x<<1)+(ch^48);
while(ch=getc(),isdigit(ch));
if(sign)x=-x;
}
template<class T,class...Types>void operator ()(T &x,Types &...args) {
return (*this)(x),(*this)(args...);
}
} I;
struct Ostream {
void putc(char c) {
putchar(c);
}
template<class T>void operator ()(T x,const char lst='\n') {
static int top(0);
static char st[50];
if(x<0)putc('-'),x=-x;
do st[++top]=(x%10)^48,x/=10;
while(x);
while(top)putc(st[top--]);
putc(lst);
}
template<class T,class...Types>void operator ()(const T x,const char lst='\n',const Types ...args) {
return (*this)(x,lst),(*this)(args...);
}
} O;
struct Estream {
template<class T>void operator ()(const char *fmt,T x) {
cerr<<fmt<<':'<<x<<".\n";
}
template<class T,class...Types>void operator ()(const char *fmt,const T x,const Types ...args) {
while(*fmt^',')cerr<<*fmt++;
return cerr<<':'<<x<<", ",(*this)(++fmt,args...);
}
} E;
} using namespace IOEstream;
namespace Modular {
#define Mod 998244353
template<class T1,class T2>constexpr auto add(const T1 a, const T2 b) {
return a+b>=Mod?a+b-Mod:a+b;
}
template<class T1,class T2>constexpr auto mul(const T1 a,const T2 b) {
return (ll)a*b%Mod;
}
template<class T,class...Types>constexpr auto add(const T a,const Types...args) {
return add(a,add(args...));
}
template<class T,class...Types>constexpr auto mul(const T a,const Types...args) {
return mul(a,mul(args...));
}
template<class T1,class T2>T1 &toadd(T1 &a,const T2 b) {
return a=add(a,b);
}
template<class T1,class T2>T1 &tomul(T1 &a,const T2 b) {
return a=mul(a,b);
}
template<class T0,class T,class...Types>T0 &toadd(T0 &a,const T b,const Types...args) {
return toadd(a,b),toadd(a,args...);
}
template<class T0,class T,class...Types>T0 &tomul(T0 &a,const T b,const Types...args) {
return tomul(a,b),tomul(a,args...);
}
} using namespace Modular;
int n,m,w;
int f[2];
struct Mat {
int n,m;
int a[3][3];
int *operator [](int i) {
return a[i];
}
Mat(int n=0,int m=0,bool unit=0):n(n),m(m) {
RCL(a,0,a,1),a[0][0]=a[1][1]=a[2][2]=unit;
}
friend Mat operator *(Mat A,Mat B) {
Mat C(A.n,B.m);
FOR(k,0,A.m-1)FOR(i,0,C.n-1)FOR(j,0,C.m-1)toadd(C[i][j],mul(A[i][k],B[k][j]));
return C;
}
friend Mat operator ^(Mat A,int b) {
Mat res(A.n,A.n,1);
for(; b; b>>=1,A=A*A)if(b&1)res=res*A;
return res;
}
} A,B;
int main() {
I(n,m),A=Mat(1,3),B=Mat(3,3);
A[0][0]=0,A[0][1]=2,A[0][2]=3;
B[0][0]=0,B[0][1]=0,B[0][2]=1;
B[1][0]=1,B[1][1]=0,B[1][2]=1;
B[2][0]=0,B[2][1]=1,B[2][2]=0;
w=(A*(B^(n-3)))[0][2],f[0]=1;
/**/
// g[1]=0,g[2]=2,g[3]=3;
// FOR(i,4,n)g[i]=add(g[i-2],g[i-3]);
// w=g[n];
/**/
FOR(i,0,lN) {
bool typ(m>>i&1);
toadd(tomul(f[typ^1],w),f[typ]);
}
O(f[0],'\n');
return 0;
}
/*
1. m&1<<i == 0
f_{i+1,0} = f_{i,0}
f_{i+1,1} = w * f_{i,1} + f_{i,0}
2. m&1<<i == 1
f_{i+1,0} = w * f_{i,0} + f_{i,1}
f_{i+1,1} = f_{i,1}
*/
/*
[f_1,f_2,f_3] = {0,2,3}
[f_2,f_3,f_4] = [f_1,f_2,f_3] *
\begin{bmatrix}
0 0 1
1 0 1
0 1 0
\end{bmatrix}
*/
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 3528kb
input:
3 2
output:
4
result:
ok "4"
Test #2:
score: 0
Accepted
time: 0ms
memory: 3600kb
input:
3 0
output:
1
result:
ok "1"
Test #3:
score: 0
Accepted
time: 0ms
memory: 3600kb
input:
100 100
output:
343406454
result:
ok "343406454"
Test #4:
score: 0
Accepted
time: 0ms
memory: 3596kb
input:
686579306 119540831
output:
260602195
result:
ok "260602195"
Test #5:
score: 0
Accepted
time: 0ms
memory: 3604kb
input:
26855095 796233790
output:
492736984
result:
ok "492736984"
Test #6:
score: 0
Accepted
time: 0ms
memory: 3532kb
input:
295310488 262950628
output:
503008241
result:
ok "503008241"
Test #7:
score: 0
Accepted
time: 0ms
memory: 3600kb
input:
239670714 149827706
output:
994702969
result:
ok "994702969"
Test #8:
score: 0
Accepted
time: 0ms
memory: 3648kb
input:
790779949 110053353
output:
898252188
result:
ok "898252188"
Test #9:
score: 0
Accepted
time: 0ms
memory: 3600kb
input:
726600542 795285932
output:
183726777
result:
ok "183726777"
Test #10:
score: 0
Accepted
time: 0ms
memory: 3660kb
input:
957970519 585582861
output:
298814018
result:
ok "298814018"
Test #11:
score: 0
Accepted
time: 0ms
memory: 3480kb
input:
93349859 634036506
output:
110693443
result:
ok "110693443"
Test #12:
score: 0
Accepted
time: 0ms
memory: 3524kb
input:
453035113 34126396
output:
306244220
result:
ok "306244220"
Test #13:
score: 0
Accepted
time: 0ms
memory: 3480kb
input:
31994526 100604502
output:
930620701
result:
ok "930620701"
Test #14:
score: 0
Accepted
time: 0ms
memory: 3652kb
input:
234760741 249817734
output:
440195858
result:
ok "440195858"
Test #15:
score: 0
Accepted
time: 0ms
memory: 3652kb
input:
542621111 646412689
output:
207377992
result:
ok "207377992"
Test #16:
score: 0
Accepted
time: 0ms
memory: 3588kb
input:
28492783 602632297
output:
65234506
result:
ok "65234506"
Test #17:
score: 0
Accepted
time: 0ms
memory: 3600kb
input:
213500301 768820204
output:
540205113
result:
ok "540205113"
Test #18:
score: 0
Accepted
time: 0ms
memory: 3524kb
input:
697808101 753041955
output:
93561295
result:
ok "93561295"
Test #19:
score: 0
Accepted
time: 0ms
memory: 3604kb
input:
585126464 450455977
output:
91109717
result:
ok "91109717"
Test #20:
score: 0
Accepted
time: 0ms
memory: 3588kb
input:
236696315 482334538
output:
925575199
result:
ok "925575199"
Test #21:
score: 0
Accepted
time: 0ms
memory: 3652kb
input:
632719214 298704996
output:
358926097
result:
ok "358926097"
Test #22:
score: 0
Accepted
time: 0ms
memory: 3528kb
input:
869119333 933404114
output:
318501108
result:
ok "318501108"
Test #23:
score: 0
Accepted
time: 0ms
memory: 3524kb
input:
6977994 814763202
output:
585332987
result:
ok "585332987"
Test #24:
score: 0
Accepted
time: 0ms
memory: 3604kb
input:
3 1
output:
3
result:
ok "3"
Test #25:
score: 0
Accepted
time: 0ms
memory: 3584kb
input:
3 1000000000
output:
279679226
result:
ok "279679226"
Test #26:
score: 0
Accepted
time: 0ms
memory: 3592kb
input:
4 0
output:
1
result:
ok "1"
Test #27:
score: 0
Accepted
time: 0ms
memory: 3600kb
input:
4 1
output:
2
result:
ok "2"
Test #28:
score: 0
Accepted
time: 0ms
memory: 3668kb
input:
4 1000000000
output:
1755648
result:
ok "1755648"
Test #29:
score: 0
Accepted
time: 0ms
memory: 3532kb
input:
1000000000 0
output:
1
result:
ok "1"
Test #30:
score: 0
Accepted
time: 0ms
memory: 3484kb
input:
1000000000 1
output:
696798313
result:
ok "696798313"
Test #31:
score: 0
Accepted
time: 0ms
memory: 3680kb
input:
1000000000 1000000000
output:
703786301
result:
ok "703786301"
Extra Test:
score: 0
Extra Test Passed