結果

問題 No.2613 Sum of Combination
ユーザー inksamuraiinksamurai
提出日時 2024-01-20 05:56:38
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 991 ms / 4,500 ms
コード長 2,862 bytes
コンパイル時間 5,141 ms
コンパイル使用メモリ 275,264 KB
実行使用メモリ 19,304 KB
最終ジャッジ日時 2024-09-28 05:27:26
合計ジャッジ時間 28,334 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 21 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 4 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 2 ms
5,376 KB
testcase_13 AC 23 ms
5,376 KB
testcase_14 AC 26 ms
5,376 KB
testcase_15 AC 18 ms
5,376 KB
testcase_16 AC 20 ms
5,376 KB
testcase_17 AC 27 ms
5,376 KB
testcase_18 AC 19 ms
5,376 KB
testcase_19 AC 32 ms
5,376 KB
testcase_20 AC 4 ms
5,376 KB
testcase_21 AC 2 ms
5,376 KB
testcase_22 AC 40 ms
5,376 KB
testcase_23 AC 834 ms
17,224 KB
testcase_24 AC 822 ms
17,056 KB
testcase_25 AC 662 ms
14,716 KB
testcase_26 AC 954 ms
19,076 KB
testcase_27 AC 465 ms
11,636 KB
testcase_28 AC 945 ms
18,844 KB
testcase_29 AC 881 ms
18,028 KB
testcase_30 AC 961 ms
19,100 KB
testcase_31 AC 868 ms
17,884 KB
testcase_32 AC 864 ms
17,320 KB
testcase_33 AC 963 ms
19,304 KB
testcase_34 AC 948 ms
19,304 KB
testcase_35 AC 979 ms
19,304 KB
testcase_36 AC 965 ms
19,176 KB
testcase_37 AC 955 ms
19,300 KB
testcase_38 AC 921 ms
18,896 KB
testcase_39 AC 944 ms
18,956 KB
testcase_40 AC 918 ms
18,680 KB
testcase_41 AC 956 ms
18,908 KB
testcase_42 AC 991 ms
18,884 KB
testcase_43 AC 953 ms
19,176 KB
testcase_44 AC 945 ms
19,304 KB
testcase_45 AC 2 ms
5,376 KB
testcase_46 AC 2 ms
5,376 KB
testcase_47 AC 2 ms
5,376 KB
testcase_48 AC 2 ms
5,376 KB
testcase_49 AC 2 ms
5,376 KB
testcase_50 AC 971 ms
19,176 KB
testcase_51 AC 937 ms
19,176 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
#define int ll
using namespace std;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define rng(i,c,n) for(int i=c;i<n;i++)
#define fi first
#define se second
#define pb push_back
#define sz(a) (int)a.size()
#define vec(...) vector<__VA_ARGS__>
#define _48SMghu ios::sync_with_stdio(0),cin.tie(0)
using ll=long long;
using pii=pair<int,int>;
using vi=vector<int>;
void print(){cout<<'\n';}
template<class h,class...t>
void print(const h&v,const t&...u){cout<<v<<' ',print(u...);}

// snuke's mod int
template <ll mod>
struct modint{
	ll x;
	modint(ll x=0):x((x%mod+mod)%mod){}
	modint operator-()const{return modint(-x);}
	modint& operator+=(const modint a){if((x+=a.x)>=mod) x-=mod; return *this;}
	modint& operator-=(const modint a){if((x+=mod-a.x)>=mod) x-=mod; return *this;}
	modint& operator*=(const modint a){(x*=a.x)%=mod; return *this;}
	modint operator+(const modint a)const{modint res(*this); return res+=a;}
	modint operator-(const modint a)const{modint res(*this); return res-=a;}
	modint operator*(const modint a)const{modint res(*this); return res*=a;}
	modint pow(ll n)const{
		modint res=1,x(*this);
		while(n){
			if(n&1)res*=x;
			x*=x;
			n>>=1;
		}
		return res;
	}
	modint inv()const{return pow(mod-2);}
};

using mint=modint<998244353>;

using vm=vec(mint);

int binpow(int a,int b,int mod){
	int res=1;
	while(b){
		if(b%2){
			res=res*a%mod;
			res%=mod;
		}
		a=a*a%mod;
		a%=mod;
		b/=2;
	}
	return res;
}

void slv(){
	int n,p;
	cin>>n>>p;

	vi a;
	int x=n;
	while(x){
		a.pb(x%p);
		x/=p;
	}
	reverse(a.begin(), a.end());

	// edgy cas
	if(p==2){
		int c=0;
		rep(i,sz(a)){
			c+=a[i];
		}
		print(mint(1ll<<c).x);
		return;
	}

	vi pri_fact;
	{
		int li=p-1;
		int x=li;
		for(int i=2;i<=sqrt(li);i++){
			int c=0;
			while(x%i==0){
				x/=i;
				c+=1;
			}
			if(c){
				pri_fact.pb(i);
			}
		}
	}

	int r=2;
	{
		while(1){
			int pok=1;
			for(auto x:pri_fact){
				pok=pok and binpow(r,(p-1)/x,p)!=1;
			}
			if(pok){
				break;
			}
			r+=1;
		}
	}

	vi ids(p);
	{
		int x=1;
		for(int i=0;i<p;i++){
			ids[x]=i;
			x=(x*r)%p;
		}
	}

	vi fact(p),invfact(p);
	{
		fact[0]=1;
		rng(i,1,p){
			fact[i]=fact[i-1]*i%p;
			fact[i]%=p;
		}
		rep(i,p){
			invfact[i]=binpow(fact[i],p-2,p);
		}
	}

	auto cnk=[&](int k,int n){
		int e=fact[n]*invfact[k]%p;
		e%=p;
		e*=invfact[n-k]%p;
		return e%p;
	};

	int si=sz(a);
	vi dp(p);
	dp[0]=1;
	rep(i,si){
		vi b(p);
		int x=a[i];
		rep(j,x+1){
			int now=cnk(j,x);
			b[ids[now]]+=1;
		}
		vi c=atcoder::convolution<998244353>(dp,b);
		vi ndp(p);
		rep(i,sz(c)){
			int y=binpow(r,i,p);
			ndp[ids[y]]+=c[i];
		}
		dp=ndp;
		// rep(j,p){
		// 	print(j,dp[j]);
		// }
		// print("~ ~");
	}
	mint ans=0;
	rep(i,p){
		ans+=mint(dp[i])*mint(binpow(r,i,p));
	}
	print(ans.x);
}

signed main(){
_48SMghu;
	slv();
}
0