結果

問題 No.2668 Trees on Graph Paper
ユーザー hamo21
提出日時 2024-03-08 23:06:24
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
MLE  
実行時間 -
コード長 2,895 bytes
コンパイル時間 4,416 ms
コンパイル使用メモリ 256,848 KB
最終ジャッジ日時 2025-02-20 03:00:57
ジャッジサーバーID
(参考情報)
judge1 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 23 MLE * 10
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include<bits/stdc++.h>
#include<atcoder/all>
using namespace std;
using namespace atcoder;
typedef int64_t ll;
using mint=modint998244353;
typedef long double ld;
const ll MOD=1000000007;
const ll MODA=998244353;
ll vx[4]={0,1,0,-1};
ll vy[4]={1,0,-1,0};
#define rep(i,n) for(ll i=0;i<(ll)(n);i++)
long long gcd(long long a,long long b){
a=abs(a);
b=abs(b);
ll gcdmax=max(a,b);
ll gcdmin=min(a,b);
if(gcdmin==0)return gcdmax;
while(true){
if(gcdmax%gcdmin==0)break;
else gcdmax%=gcdmin;
swap(gcdmin,gcdmax);
}
return gcdmin;
}
ll pow(ll N,ll P,ll M){
if(P==0)return 1;
else if(P%2==0){
ll t=pow(N,P/2,M);
return t*t%M;
}
else return N*pow(N,P-1,M)%M;
}
ll pow(ll N,ll P){
if(P==0)return 1;
else if(P%2==0){
ll t=pow(N,P/2);
return t*t;
}
else return N*pow(N,P-1);
}
vector<ll> fac;
vector<ll> finv;
vector<ll> inv;
void COMinit(ll N,ll P){
rep(i,N+1){
if(i==0){
fac.push_back(1);
finv.push_back(1);
inv.push_back(1);
}
else if(i==1){
fac.push_back(1);
finv.push_back(1);
inv.push_back(1);
}
else{
fac.push_back(fac.at(i-1)*i%P);
inv.push_back(P-inv.at(P%i)*(P/i)%P);
finv.push_back(finv.at(i-1)*inv.at(i)%P);
}
}
}
ll COM(ll n,ll k,ll P){
if(n<k)return 0;
if(n<0||k<0)return 0;
return fac.at(n)*(finv.at(k)*finv.at(n-k)%P)%P;
}
struct UnionFind {
vector<ll> par; // par[i]:i () par[3] = 2 : 32
UnionFind(ll N) : par(N) { //
for(ll i = 0; i < N; i++) par[i] = i;
}
ll root(ll x) { // xroot(x) = {x}
if (par[x] == x) return x;
return par[x] = root(par[x]);
}
void unite(ll x, ll y) { // xy
ll rx = root(x); //xrx
ll ry = root(y); //yry
if (rx == ry) return; //xy(=)
par[rx] = ry; //xy(=)xrxyry
}
bool same(ll x, ll y) { // 2x, ytrue
ll rx = root(x);
ll ry = root(y);
return rx == ry;
}
};
int main(){
ll n,m;
cin>>n>>m;
COMinit(2*n,m);
ll v=1%m;
vector<vector<ll>> dp(2*n,vector<ll>(3));
dp[0][0]=1%m;
rep(i,2*n-1){
dp[i+1][0]=dp[i][0]+(i+1)*dp[i][1];
dp[i+1][1]=dp[i][0]+dp[i][1]+dp[i][2]*(i+1);
dp[i+1][2]=dp[i][0]+dp[i][1]+dp[i][2];
rep(j,3)dp[i+1][j]%=m;
}
rep(i,n-2)v=(v*dp[2*i+1][2])%m;
cout<<(dp[2*n-3][2]*v%m)*fac[2*n-1]%m<<endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0