結果

問題 No.2688 Cell Proliferation (Hard)
ユーザー kotatsugamekotatsugame
提出日時 2024-03-20 22:10:47
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,718 ms / 4,000 ms
コード長 925 bytes
コンパイル時間 2,249 ms
コンパイル使用メモリ 112,564 KB
実行使用メモリ 38,432 KB
最終ジャッジ日時 2024-03-20 22:11:19
合計ジャッジ時間 28,198 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
11,840 KB
testcase_01 AC 5 ms
11,840 KB
testcase_02 AC 5 ms
11,840 KB
testcase_03 AC 83 ms
13,664 KB
testcase_04 AC 1,641 ms
37,008 KB
testcase_05 AC 804 ms
25,036 KB
testcase_06 AC 352 ms
18,144 KB
testcase_07 AC 344 ms
18,148 KB
testcase_08 AC 1,611 ms
36,356 KB
testcase_09 AC 828 ms
27,468 KB
testcase_10 AC 1,658 ms
37,592 KB
testcase_11 AC 1,625 ms
36,660 KB
testcase_12 AC 359 ms
18,660 KB
testcase_13 AC 1,718 ms
38,424 KB
testcase_14 AC 819 ms
26,380 KB
testcase_15 AC 1,718 ms
38,432 KB
testcase_16 AC 788 ms
26,176 KB
testcase_17 AC 782 ms
25,680 KB
testcase_18 AC 1,671 ms
37,744 KB
testcase_19 AC 785 ms
25,908 KB
testcase_20 AC 369 ms
19,168 KB
testcase_21 AC 371 ms
19,040 KB
testcase_22 AC 796 ms
25,476 KB
testcase_23 AC 384 ms
19,552 KB
testcase_24 AC 1,659 ms
37,448 KB
testcase_25 AC 360 ms
18,788 KB
testcase_26 AC 1,697 ms
37,912 KB
testcase_27 AC 784 ms
25,828 KB
testcase_28 AC 816 ms
26,240 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<cassert>
#include<atcoder/modint>
#include<atcoder/convolution>
using namespace std;
using mint=atcoder::modint998244353;
int P1,P2,Q1,Q2,T;
mint P,Q;
mint Qt[1<<20];
mint dp[1<<20];
void solve(int l,int r)
{
	if(r-l<=100)
	{
		for(int i=l;i<r;i++)
		{
			for(int j=i+1;j<=r;j++)
			{
				dp[j]+=dp[i]*P*Qt[j-i-1];
			}
		}
		return;
	}
	int mid=(l+r)/2;
	solve(l,mid);
	vector<mint>L(mid-l+1);
	for(int i=l;i<=mid;i++)L[i-l]=dp[i]*P;
	vector<mint>R(r-l+1);
	for(int i=1;i<=r-l;i++)R[i]=Qt[i-1];
	L=atcoder::convolution(L,R);
	for(int i=mid+1;i<=r;i++)dp[i]+=L[i-l];
	solve(mid+1,r);
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	cin>>P1>>P2>>Q1>>Q2>>T;
	P=mint(P1)/P2,Q=mint(Q1)/Q2;
	{
		Qt[0]=1;
		mint v=1;
		for(int i=1;i<=T;i++)
		{
			v*=Q;
			Qt[i]=Qt[i-1]*v;
		}
	}
	dp[0]=1;
	solve(0,T);
	mint ans=0;
	for(int i=0;i<=T;i++)ans+=dp[i]*Qt[T-i];
	cout<<ans.val()<<endl;
}
0