結果
問題 | No.535 自然数の収納方法 |
ユーザー |
![]() |
提出日時 | 2018-12-04 10:54:53 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 28 ms / 2,000 ms |
コード長 | 1,911 bytes |
コンパイル時間 | 2,245 ms |
コンパイル使用メモリ | 198,324 KB |
最終ジャッジ日時 | 2025-01-06 18:07:37 |
ジャッジサーバーID (参考情報) |
judge3 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 23 |
ソースコード
#include<bits/stdc++.h>using namespace std;using Int = long long;template<typename T1,typename T2> inline void chmin(T1 &a,T2 b){if(a>b) a=b;}template<typename T1,typename T2> inline void chmax(T1 &a,T2 b){if(a<b) a=b;}template<typename T,T MOD = 1000000007>struct Mint{T v;Mint():v(0){}Mint(signed v):v(v){}Mint(long long t){v=t%MOD;if(v<0) v+=MOD;}Mint pow(long long k){Mint res(1),tmp(v);while(k){if(k&1) res*=tmp;tmp*=tmp;k>>=1;}return res;}Mint inv(){return pow(MOD-2);}Mint& operator+=(Mint a){v+=a.v;if(v>=MOD)v-=MOD;return *this;}Mint& operator-=(Mint a){v+=MOD-a.v;if(v>=MOD)v-=MOD;return *this;}Mint& operator*=(Mint a){v=1LL*v*a.v%MOD;return *this;}Mint& operator/=(Mint a){return (*this)*=a.inv();}Mint operator+(Mint a) const{return Mint(v)+=a;};Mint operator-(Mint a) const{return Mint(v)-=a;};Mint operator*(Mint a) const{return Mint(v)*=a;};Mint operator/(Mint a) const{return Mint(v)/=a;};Mint operator-(){return v?MOD-v:v;}bool operator==(const Mint a)const{return v==a.v;}bool operator!=(const Mint a)const{return v!=a.v;}bool operator <(const Mint a)const{return v <a.v;}};//INSERT ABOVE HEREsigned main(){int n;cin>>n;using M = Mint<int>;vector<M> dp(n+1,1);dp[0]=M(0);for(int i=n-2;i>=0;i--){vector<M> sm(n+1,0);for(int j=1;j<=n;j++) sm[j]=sm[j-1]+dp[j];int k=i?i:1;vector<M> nx(n+1,0);for(int j=1;j<=n;j++) nx[j]=sm[n]-sm[max(0,j-k)];swap(dp,nx);}M ans;for(int j=1;j<=n;j++) ans+=dp[j];for(int j=1;j<=n;j++) dp[j]=0;dp[n]=M(1);for(int i=n-2;i>=0;i--){vector<M> sm(n+1,0);for(int j=1;j<=n;j++) sm[j]=sm[j-1]+dp[j];int k=i?i:1;vector<M> nx(n+1,0);for(int j=1;j<=n;j++) nx[j]=sm[n]-sm[max(0,j-k)];swap(dp,nx);}ans-=dp[1];cout<<ans.v<<endl;return 0;}