QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#104520 | #59. Determinant of A+Bz | AFewSuns | Judging | / | / | C++14 | 3.5kb | 2023-05-10 22:15:10 | 2024-03-04 04:40:49 |
Judging History
你现在查看的是测评时间为 2024-03-04 04:40:49 的历史记录
- [2024-05-05 11:54:08]
- hack成功,自动添加数据
- (/hack/617)
- [2024-05-05 11:38:15]
- hack成功,自动添加数据
- (/hack/616)
- [2024-03-04 04:35:40]
- hack成功,自动添加数据
- (/hack/552)
- [2023-08-10 23:21:45]
- System Update: QOJ starts to keep a history of the judgings of all the submissions.
- [2023-05-10 22:15:10]
- 提交
answer
#include<bits/stdc++.h>
using namespace std;
namespace my_std{
#define ll long long
#define bl bool
ll my_pow(ll a,ll b,ll mod){
ll res=1;
if(!b) return 1;
while(b){
if(b&1) res=(res*a)%mod;
a=(a*a)%mod;
b>>=1;
}
return res;
}
ll qpow(ll a,ll b){
ll res=1;
if(!b) return 1;
while(b){
if(b&1) res*=a;
a*=a;
b>>=1;
}
return res;
}
#define db double
#define pf printf
#define pc putchar
#define fr(i,x,y) for(register ll i=(x);i<=(y);i++)
#define pfr(i,x,y) for(register ll i=(x);i>=(y);i--)
#define go(u) for(ll i=head[u];i;i=e[i].nxt)
#define enter pc('\n')
#define space pc(' ')
#define fir first
#define sec second
#define MP make_pair
#define il inline
#define inf 8e18
#define random(x) rand()*rand()%(x)
#define inv(a,mod) my_pow((a),(mod-2),(mod))
il ll read(){
ll sum=0,f=1;
char ch=0;
while(!isdigit(ch)){
if(ch=='-') f=-1;
ch=getchar();
}
while(isdigit(ch)){
sum=sum*10+(ch^48);
ch=getchar();
}
return sum*f;
}
il void write(ll x){
if(x<0){
x=-x;
pc('-');
}
if(x>9) write(x/10);
pc(x%10+'0');
}
il void writeln(ll x){
write(x);
enter;
}
il void writesp(ll x){
write(x);
space;
}
}
using namespace my_std;
#define mod 998244353
ll n,A[550][550],B[550][550],f[550][550],nz=0,ans=1;
bl ck=0;
void gauss1(){
fr(i,1,n){
ll tmp=i;
fr(j,i,n) if(B[j][i]) tmp=j;
if(i!=tmp){
swap(A[i],A[tmp]);
swap(B[i],B[tmp]);
ans=(mod-ans)%mod;
}
tmp=i;
fr(j,i,n) if(B[i][j]) tmp=j;
if(i!=tmp){
fr(j,1,n) swap(A[j][tmp],A[j][i]);
fr(j,1,n) swap(B[j][tmp],B[j][i]);
ans=(mod-ans)%mod;
}
if(!B[i][i]){
nz++;
swap(A[i],B[i]);
}
fr(j,1,i-1){
tmp=B[i][j];
fr(k,1,n) B[i][k]=(B[i][k]-tmp*B[j][k]%mod+mod)%mod;
fr(k,1,n) A[i][k]=(A[i][k]-tmp*A[j][k]%mod+mod)%mod;
}
tmp=i;
fr(j,i,n) if(B[i][j]) tmp=j;
if(i!=tmp){
fr(j,1,n) swap(A[j][tmp],A[j][i]);
fr(j,1,n) swap(B[j][tmp],B[j][i]);
ans=(mod-ans)%mod;
}
if(!B[i][i]){
ck=1;
return;
}
tmp=inv(B[i][i],mod);
ans=ans*B[i][i]%mod;
fr(j,1,n) B[i][j]=B[i][j]*tmp%mod;
fr(j,1,n) A[i][j]=A[i][j]*tmp%mod;
fr(j,i+1,n){
tmp=B[j][i];
fr(k,1,n) B[j][k]=(B[j][k]-tmp*B[i][k]%mod+mod)%mod;
fr(k,1,n) A[j][k]=(A[j][k]-tmp*A[i][k]%mod+mod)%mod;
}
}
pfr(i,n,1){
pfr(j,i-1,1){
ll tmp=B[j][i];
fr(k,1,n) B[j][k]=(B[j][k]-tmp*B[i][k]%mod+mod)%mod;
fr(k,1,n) A[j][k]=(A[j][k]-tmp*A[i][k]%mod+mod)%mod;
}
}
}
void gauss2(){
fr(i,2,n){
ll tmp=i;
fr(j,i,n) if(A[j][i-1]) tmp=j;
if(i!=tmp){
swap(A[i],A[tmp]);
fr(k,1,n) swap(A[k][i],A[k][tmp]);
}
if(!A[i][i-1]) continue;
ll invv=inv(A[i][i-1],mod);
fr(j,i+1,n){
tmp=A[j][i-1]*invv%mod;
fr(k,i-1,n) A[j][k]=(A[j][k]-tmp*A[i][k]%mod+mod)%mod;
fr(k,1,n) A[k][i]=(A[k][i]+tmp*A[k][j]%mod)%mod;
}
}
}
int main(){
n=read();
fr(i,1,n) fr(j,1,n) A[i][j]=read();
fr(i,1,n) fr(j,1,n) B[i][j]=read();
gauss1();
if(ck){
fr(i,0,n) writesp(0);
return 0;
}
gauss2();
f[0][0]=1;
fr(i,1,n){
fr(j,0,n) f[i][j+1]=(f[i][j+1]+f[i-1][j])%mod;
fr(j,0,n) f[i][j]=(f[i][j]+A[i][i]*f[i-1][j]%mod)%mod;
ll mul=1;
pfr(j,i,2){
mul=mul*A[j][j-1]%mod;
fr(k,0,n){
if((i-j+1)&1) f[i][k]=(f[i][k]-mul*A[j-1][i]%mod*f[j-2][k]%mod+mod)%mod;
else f[i][k]=(f[i][k]+mul*A[j-1][i]%mod*f[j-2][k]%mod)%mod;
}
}
}
fr(i,nz,n) writesp(f[n][i]*ans%mod);
fr(i,1,nz) writesp(0);
}