結果

問題 No.2688 Cell Proliferation (Hard)
ユーザー kotatsugame
提出日時 2024-03-20 22:10:47
言語 C++14
(gcc 13.3.0 + boost 1.87.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

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