結果

問題 No.2668 Trees on Graph Paper
ユーザー hamo21hamo21
提出日時 2024-03-08 23:11:35
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,529 ms / 3,000 ms
コード長 2,766 bytes
コンパイル時間 3,925 ms
コンパイル使用メモリ 265,232 KB
実行使用メモリ 268,516 KB
最終ジャッジ日時 2024-09-29 20:39:05
合計ジャッジ時間 19,863 ms
ジャッジサーバーID
(参考情報)
judge1 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 33
権限があれば一括ダウンロードができます

ソースコード

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;
ll a=1%m,b=0,c=0;
rep(i,2*n-3){
ll na=(a+(i+1)*b)%m;
ll nb=(a+b+c*(i+1))%m;
ll nc=(a+b+c)%m;
if(i%2==0)v=(v*nc)%m;
a=na,b=nb,c=nc;
}
cout<<(v*fac[2*n-1])%m<<endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0