結果
問題 | No.287 場合の数 |
ユーザー |
|
提出日時 | 2015-10-13 14:15:26 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 3 ms / 5,000 ms |
コード長 | 1,841 bytes |
コンパイル時間 | 1,490 ms |
コンパイル使用メモリ | 159,860 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-08 00:01:00 |
合計ジャッジ時間 | 2,208 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define all(c) (c).begin(),(c).end()#define rrep(i,n) for(int i=(int)(n)-1;i>=0;i--)#define REP(i,m,n) for(int i=(int)(m);i<(int)(n);i++)#define rep(i,n) REP(i,0,n)#define iter(c) __typeof((c).begin())#define tr(it,c) for(iter(c) it=(c).begin();it!=(c).end();it++)#define mem(a) memset(a,0,sizeof(a))#define pd(a) printf("%.10f\n",a)#define pb(a) push_back(a)#define in(a) insert(a)#define pi M_PI#define R cin>>#define F first#define S second#define C class#define ll long long#define ln cout<<'\n'#define _(_1,_2,_3,N,...)N#define pr(...) _(__VA_ARGS__,pr3,pr2,pr1)(__VA_ARGS__)template<C T>void pr1(T a){cout<<a;ln;}template<C T,C T2>void pr2(T a,T2 b){cout<<a<<' '<<b;ln;}template<C T,C T2,C T3>void pr3(T a,T2 b,T3 c){cout<<a<<' '<<b<<' '<<c;ln;}template<C T>void PR(T a,int n){rep(i,n){if(i)cout<<' ';cout<<a[i];}ln;}bool check(int n,int m,int x,int y){return x>=0&&x<n&&y>=0&&y<m;}const ll MAX=1000000007,MAXL=1LL<<60,dx[4]={-1,0,1,0},dy[4]={0,1,0,-1};const int N=1<<11;const double th=2*pi/N;typedef complex<double> P;void FFT(P a[],double t) {for(int m=N; m>=2; m>>=1) {int mh=m>>1;for(int i=0; i<mh; i++) {P w=exp(t*i*P(0,1));for(int j=i; j<N; j+=m) {int k=j+mh;P x=a[j]-a[k];a[j]+=a[k];a[k]=w*x;}}t*=2;}int i=0;for(int j=1; j<N-1; j++) {for(int k=N>>1; k>(i^=k); k>>=1);if(j<i) swap(a[i],a[j]);}}P a[N];ll d[N];void solve() {for(int i=0; i<N; i++) a[i]=P(d[i],0);FFT(a,th);for(int i=0; i<N; i++) a[i]*=a[i];FFT(a,-th);for(int i=0; i<N; i++) d[i]=round(a[i].real()/N);}void Main() {int n;R n;rep(i,N) d[i]=i<=n;rep(i,3) solve();pr(d[6*n]);}int main() {ios::sync_with_stdio(0);cin.tie(0);Main();return 0;}