結果

問題 No.1907 DETERMINATION
ユーザー Thangtozy6Thangtozy6
提出日時 2023-10-12 12:01:04
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 5,063 bytes
コンパイル時間 2,261 ms
コンパイル使用メモリ 174,592 KB
実行使用メモリ 10,732 KB
最終ジャッジ日時 2023-10-12 12:01:56
合計ジャッジ時間 41,678 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
7,552 KB
testcase_01 AC 3 ms
7,556 KB
testcase_02 AC 2 ms
7,436 KB
testcase_03 AC 3 ms
7,544 KB
testcase_04 AC 3 ms
7,560 KB
testcase_05 AC 2 ms
7,628 KB
testcase_06 AC 3 ms
7,472 KB
testcase_07 AC 360 ms
9,640 KB
testcase_08 AC 148 ms
8,988 KB
testcase_09 AC 251 ms
9,512 KB
testcase_10 AC 761 ms
10,396 KB
testcase_11 AC 563 ms
10,596 KB
testcase_12 AC 836 ms
10,424 KB
testcase_13 AC 791 ms
10,472 KB
testcase_14 AC 754 ms
10,452 KB
testcase_15 WA -
testcase_16 AC 52 ms
8,492 KB
testcase_17 AC 721 ms
10,372 KB
testcase_18 AC 511 ms
10,000 KB
testcase_19 AC 16 ms
8,036 KB
testcase_20 AC 777 ms
10,408 KB
testcase_21 AC 76 ms
8,816 KB
testcase_22 AC 742 ms
10,608 KB
testcase_23 AC 798 ms
10,404 KB
testcase_24 AC 250 ms
9,408 KB
testcase_25 AC 2 ms
7,540 KB
testcase_26 AC 872 ms
10,568 KB
testcase_27 AC 868 ms
10,524 KB
testcase_28 AC 867 ms
10,500 KB
testcase_29 AC 873 ms
10,456 KB
testcase_30 AC 2 ms
7,592 KB
testcase_31 AC 871 ms
10,460 KB
testcase_32 AC 868 ms
10,540 KB
testcase_33 AC 873 ms
10,580 KB
testcase_34 AC 895 ms
10,444 KB
testcase_35 AC 3 ms
7,556 KB
testcase_36 AC 3 ms
7,564 KB
testcase_37 AC 3 ms
7,556 KB
testcase_38 AC 874 ms
10,484 KB
testcase_39 AC 875 ms
10,500 KB
testcase_40 AC 871 ms
10,524 KB
testcase_41 AC 873 ms
10,568 KB
testcase_42 AC 873 ms
10,448 KB
testcase_43 AC 872 ms
10,640 KB
testcase_44 AC 872 ms
10,684 KB
testcase_45 AC 871 ms
10,716 KB
testcase_46 AC 863 ms
10,440 KB
testcase_47 AC 871 ms
10,468 KB
testcase_48 AC 873 ms
10,568 KB
testcase_49 AC 867 ms
10,552 KB
testcase_50 AC 872 ms
10,476 KB
testcase_51 AC 872 ms
10,552 KB
testcase_52 AC 3 ms
7,512 KB
testcase_53 AC 873 ms
10,556 KB
testcase_54 AC 873 ms
10,488 KB
testcase_55 AC 2 ms
7,496 KB
testcase_56 AC 873 ms
10,500 KB
testcase_57 AC 872 ms
10,460 KB
testcase_58 AC 660 ms
10,444 KB
testcase_59 AC 673 ms
10,732 KB
testcase_60 AC 678 ms
10,468 KB
testcase_61 AC 733 ms
10,496 KB
testcase_62 AC 678 ms
10,480 KB
testcase_63 AC 875 ms
10,444 KB
testcase_64 AC 3 ms
7,444 KB
testcase_65 AC 3 ms
7,580 KB
testcase_66 AC 2 ms
7,408 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std; 
const int mod=998244353;
class mint{
	private:
		int x;
		long long ksm(long long w,int y)
		{
			if(!y)return 1;
			if(y%2==0)return ksm(w*w%mod,y>>1);
			return w*ksm(w*w%mod,y>>1)%mod;
		}
	public:
		void read()
		{
			int fl=0;x=0;
			char ch=getchar();
			while(ch<'0'||ch>'9')
			{
				if(ch=='-')fl=1;
				ch=getchar();
			}
			while('0'<=ch&&ch<='9')
			{
				x=(10ll*x+ch-'0')%mod;
				ch=getchar();
			}
			if(fl&&x)ch=mod-x;
		}
		void print(){cout<<x;}
		void operator=(long long y)
		{
			x=y%mod;
			if(x<0)x+=mod;
		}
		mint operator+(mint y)const{mint a;a.x=x+y.x;if(a.x>=mod)a.x-=mod;return a;}
		mint operator~(void){mint a;a.x=x;if(a.x)a.x=mod-a.x;return a;}
		mint operator-(mint y)const{mint a;a=x;return a+(~y);}
		mint operator*(mint y)const{mint a;a.x=1ll*x*y.x%mod;return a;}
		mint inv(){mint a;a.x=ksm(x,mod-2);return a;}
		mint operator/(mint y)const{mint a;a=x;return a*(y.inv());}
		int val(){return x;}
		mint &operator++(void){x++;if(x==mod)x=0;return *this;}
		mint &operator--(void){x--;if(x==-1)x=mod-1;return *this;}
		mint &operator+=(mint y){x+=y.x;if(x>=mod)x-=mod;return *this;}
		mint &operator-=(mint y){x-=y.x;if(x<0)x+=mod;return *this;}
		mint &operator*=(mint y){x=1ll*x*y.x%mod;return *this;}
		mint &operator/=(mint y){x=1ll*x*(y.inv()).x%mod;return *this;}
		friend istream &operator >>(istream&,mint&);
		friend ostream &operator <<(ostream&,mint);
};
istream &operator >>(istream& input,mint& x){
	x.read();
	return input;
}
ostream &operator <<(ostream& output,mint x){
	x.print();
	return output;
}
mint m0[505][505],m1[505][505];
mint a[505][505];
mint pl[505][505];
mint ans[505];
//求a的特征多项式 
void char_poly(int n)
{
	//hessenberg
	for(int i=1;i<=n;i++)
	{
		int p=-1;
		for(int j=i+1;j<=n;j++)
		{
			if(a[j][i].val())
			{
				p=j;
				break;
			}
		}
		if(p==-1)continue;
		swap(a[i+1],a[p]);
		for(int j=1;j<=n;j++)swap(a[j][i+1],a[j][p]);
		for(int j=i+2;j<=n;j++)
		{
			mint v=a[j][i]/a[i+1][i];
			for(int k=1;k<=n;k++)a[j][k]-=a[i+1][k]*v;
			for(int k=1;k<=n;k++)a[k][i+1]+=a[k][j]*v;
		}
	}
	pl[0][0]=1;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<i;j++)
		{
			mint v=~a[j][i];
			for(int k=j+1;k<=i;k++)v*=a[k][k-1];
			for(int k=0;k<=n;k++)pl[i][k]+=pl[j-1][k]*v;
		}
		for(int j=1;j<=i;j++)pl[i][j]+=pl[i-1][j-1];
		for(int j=0;j<=i;j++)pl[i][j]-=pl[i-1][j]*a[i][i];
	}
	for(int i=0;i<=n;i++)ans[i]=pl[n][i];
}
mint g[505][505];
mint h[505][505];
//h=g^-1
void matriv(int n)
{
	for(int i=1;i<=n;i++)h[i][i]=1;
	for(int i=1;i<=n;i++)
	{
		int p=-1;
		for(int j=i;j<=n;j++)
		{
			if(g[i][j].val())
			{
				p=j;
				break;
			}
		}
		for(int j=1;j<=n;j++)
		{
			swap(g[i][j],g[p][j]);
			swap(h[i][j],h[p][j]); 
		}
		mint iv=g[i][i].inv();
		for(int j=1;j<=n;j++)
		{
			g[i][j]=g[i][j]*iv;
			h[i][j]=h[i][j]*iv;
		}
		if(p==-1)continue;
		for(int j=i+1;j<=n;j++)
		{
			mint v=g[j][i];
			for(int k=1;k<=n;k++)
			{
				g[j][k]-=g[i][k]*v;
				h[j][k]-=h[i][k]*v;
			}
		}
	}
	for(int i=n;i>=1;i--)
	{
		for(int j=i-1;j>=1;j--)
		{
			mint v=g[j][i];
			for(int k=1;k<=n;k++)
			{
				g[j][k]-=g[i][k]*v;
				h[j][k]-=h[i][k]*v;
			}
		}
	}
}
mint r[505][505];
//r=g*h
void product(int n)
{
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			r[i][j]=0;
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			for(int k=1;k<=n;k++)
			{
				r[i][j]+=g[i][k]*h[k][j];
			}
		}
	}
}
//det of h
mint det(int n)
{
	mint ret;
	ret=1;
	for(int i=1;i<=n;i++)
	{
		int p=-1;
		for(int j=1;j<=n;j++)
		{
			if(h[i][j].val())
			{
				p=j;
				break;
			}
		}
		if(p==-1)continue;
		if(p!=i)
		{
			swap(h[i],h[p]);
			ret=~ret;
		}
		for(int j=i+1;j<=n;j++)
		{
			mint v=h[j][i]/h[i][i];
			for(int k=1;k<=n;k++)
			{
				h[j][k]-=h[i][k]*v;
			}
		}
	}
	for(int i=1;i<=n;i++)ret*=h[i][i];
	return ret;
}
mt19937 rnd(time(0));
mint sum[505];
mint cho[505][505];
mint mi[505];
int main()
{
//	freopen("ex_5.in","r",stdin);
	int n;
	std::ios::sync_with_stdio(false);
	std::cin.tie(nullptr);
	cin>>n;
	int x;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			cin>>x;
			m0[i][j]=x;
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			cin>>x;
			m1[i][j]=x;
		}
	}
	mint w;
	w=rnd()%(mod-1)+1;
	w=63504649;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			g[i][j]=m0[i][j]-w*m1[i][j];
		}
	}
	matriv(n);
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			g[i][j]=m1[i][j];
		}
	}
	product(n);
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			a[i][j]=~r[i][j];
		}
	}
	char_poly(n);
	mint _v=det(n).inv();
	for(int i=0;i<=n;i++)sum[i]=ans[n-i]*_v;
//	for(int i=0;i<=n;i++)cout<<sum[i]<<" ";
//	cout<<endl;
	cho[0][0]=1;
	for(int i=1;i<=n;i++)
	{
		cho[i][0]=1;
		for(int j=1;j<=i;j++)
		{
			cho[i][j]=cho[i-1][j-1]+cho[i-1][j];
		}
	}
	mi[0]=1;
	for(int i=1;i<=n;i++)mi[i]=mi[i-1]*w;
	for(int i=1;i<=n;i++)ans[i]=0;
	ans[0]=sum[0];
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<=i;j++)
		{
			ans[j]+=mi[i-j]*cho[i][j]*sum[i];
		}
	}
	for(int i=0;i<=n;i++)cout<<ans[i].val()<<"\n";
	return 0;
}
0