結果

問題 No.2688 Cell Proliferation (Hard)
ユーザー kotatsugamekotatsugame
提出日時 2024-03-20 22:10:47
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,472 ms / 4,000 ms
コード長 925 bytes
コンパイル時間 2,174 ms
コンパイル使用メモリ 112,804 KB
実行使用メモリ 38,244 KB
最終ジャッジ日時 2024-09-30 08:05:06
合計ジャッジ時間 24,627 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 4 ms
11,648 KB
testcase_01 AC 4 ms
11,588 KB
testcase_02 AC 4 ms
11,520 KB
testcase_03 AC 70 ms
13,412 KB
testcase_04 AC 1,391 ms
36,708 KB
testcase_05 AC 670 ms
24,812 KB
testcase_06 AC 305 ms
18,020 KB
testcase_07 AC 297 ms
18,028 KB
testcase_08 AC 1,368 ms
36,096 KB
testcase_09 AC 714 ms
27,236 KB
testcase_10 AC 1,409 ms
37,348 KB
testcase_11 AC 1,376 ms
36,412 KB
testcase_12 AC 323 ms
18,536 KB
testcase_13 AC 1,456 ms
38,116 KB
testcase_14 AC 688 ms
26,128 KB
testcase_15 AC 1,464 ms
38,244 KB
testcase_16 AC 679 ms
25,960 KB
testcase_17 AC 678 ms
25,328 KB
testcase_18 AC 1,424 ms
37,492 KB
testcase_19 AC 680 ms
25,656 KB
testcase_20 AC 319 ms
19,044 KB
testcase_21 AC 321 ms
18,920 KB
testcase_22 AC 684 ms
25,188 KB
testcase_23 AC 333 ms
19,428 KB
testcase_24 AC 1,413 ms
37,228 KB
testcase_25 AC 319 ms
18,660 KB
testcase_26 AC 1,472 ms
37,568 KB
testcase_27 AC 696 ms
25,452 KB
testcase_28 AC 724 ms
25,988 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