StarsAndBars1

#include<iostream>
#define NMax 10000
using namespace std;
typedef short Huge[NMax+3];
int n, k;
void AtribValue(Huge H,short X)
{
    H[0] = 0;
    while (X)
    {
      ++H[0];
      H[H[0]] = X % 10;
      X /= 10;
    }
}
void Mult(Huge H,short X)
{
    short i,T=0;
    for (i=1;i<=H[0];i++)
    {
        H[i]=H[i]*X+T;
        T=H[i]/10;
        H[i]=H[i]%10;
    }
    while (T)
    { H[++H[0]]=T%10;
      T/=10;
    }
}

void Divide(Huge A,short X)
{
    short i,R=0;
    for (i=A[0];i;i--)
    {
        A[i]=(R=10*R+A[i])/X;
        R%=X;
    }
    while (!A[A[0]] && A[0]>1) A[0]--;
}

void Afisez(Huge H)
{
    for(short i=H[0];i>0;--i)
        cout<<H[i];
    cout<<'\n';
}

void BinComb(unsigned int n,unsigned int k)
{
    Huge res;
    AtribValue(res,1);
    if(k>n-k)
        k=n-k;
    for(unsigned int i=0;i<k;++i)
    {
        Mult(res,n-i);
        Divide(res,i+1);
    }
    Afisez(res);
}

int main()
{
    int n;
    cin>>n>>k;
    n=n+k-1;
    k--;
    BinComb(n,k);
    return 0;
}