/*    ∫ ∫ ∫    ノヽ   (_  )  (_    ) (______ )  ヽ(´・ω・)ノ     |  /    UU */ #pragma region macro #include typedef long long int64; using namespace std; typedef vector vi; const int MOD = (int)1e9 + 7; const int64 INF = 1LL << 62; const int inf = 1<<30; const char bn = '\n'; templatebool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (b ostream& operator<<(ostream& os, const vector &V){ int N = V.size(); REP(i,N){ os << V[i]; if (i!=N-1) os << " "; } os << "\n"; return os; } template ostream& operator<<(ostream& os, pair const&P){ os << "("; os << P.first; os << " , "; os << P.second; os << ")"; return os; } template ostream& operator<<(ostream& os, set &S){ auto it=S.begin(); while(it!=S.end()){ os << *it; os << " "; it++; } os << "\n"; return os; } template ostream& operator<<(ostream& os, deque &q){ for(auto it=q.begin();it ostream& operator<<(ostream& os, map const&M){ for(auto e:M){ os<> dxdy = {mp(0,1),mp(1,0),mp(-1,0),mp(0,-1)}; #pragma endregion //fixed< bit; for(b=b;b>0;b>>=1){ bit.push_back(b&1); } vector fac(bit.size()); fac[0] = a; int64 res = 1; for(int i=1;i= MOD) x -= MOD; return *this; } mint& operator-=(const mint a) { if ((x += MOD-a.x) >= MOD) x -= MOD; return *this; } mint& operator*=(const mint a) { (x *= a.x) %= MOD; return *this; } mint operator+(const mint a) const { mint res(*this); return res+=a; } mint operator-(const mint a) const { mint res(*this); return res-=a; } mint operator*(const mint a) const { mint res(*this); return res*=a; } mint pow(int64 t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } // for prime MOD mint inv() const { return pow(MOD-2); } mint& operator/=(const mint a) { return (*this) *= a.inv(); } mint operator/(const mint a) const { mint res(*this); return res/=a; } }; ostream& operator<<(ostream& os, mint a){ os << a.x; return os; } class BigCombination{ private: const int MAX = 510'000; mint fac[510000], inv[510000]; public: BigCombination(){ fac[0] = fac[1] = 1; for(int i=2; i0;i--){ inv[i-1] = inv[i] * i; } } mint combination(int n,int r){ if(n> S(n+1,vector(r+1,0)); //S[i][j]:=i人を区別されないj個のグループに分ける場合の数 for(int i=0;i<=min(n,r);i++) S[i][i]=1; for(int i=0;i<=n;i++) S[i][1]=1; for(int i=2;i<=n;i++){ for(int j=2;j<=r;j++){ S[i][j] = S[i-1][j-1] + S[i-1][j] * j; } } mint res = 0; for(int j=0;j<=r;j++) res+=S[n][j]; return res; } //区別できるボール/区別できない箱/1個以下 mint ball_8(int n,int r){ return n<=r; } // //区別できるボール/区別できない箱/1個以上 スターリング数 // mint ball_9(int n,int r){ // BigCombination B; // mint res = 0; // for(int i=0;i> P(n+1,vector(r+1,0)); REP(j,r+1) P[0][j]=1; for(int i=1;i<=n;i++){ for(int j=1;j<=r;j++){ P[i][j] = P[i][j-1]; //0を含む if (i-j>=0) P[i][j] += P[i-j][j]; //0を含まない } } return P[n][r]; } //区別できないボール/区別できない箱/1個以下 mint ball_11(int n,int r){ return n<=r; } //区別できないボール/区別できない箱/1個以上 分割数 mint ball_12(int n,int r){ if (n-r<0) return 0; return ball_10(n-r,r); } int main(){ cin.tie(nullptr); ios::sync_with_stdio(false); int N,K; cin >> N >> K; auto ans = ball_3(N,K); cout << ans << endl; }