Вычисление (2n)!/(n!*n!)

   long long KPL = 2;  // Кол-во позиций
   for(i=N+1; i<2*N; i++) KPL *= i;
   for(i=2; i< N; i++) KPL /= i;
   long KP = (LU)KPL;

int nod(int a, int b)
{  int t;
   while(a) {
     if (a < b) {
       t = a;
       a = b;
       b = t;
     }
     a = a % b;
   }
   return b;
}
// ---------
unsigned long N2N(int n)  // (2n)! / (n!*n!)
{ int chis[40], znam[40];  // Сомножители числителя и знаменателя
  int i, j;
  unsigned long res = 1;
   for(i=0; i<n; i++) {
     chis[i] = 2*n - i;
     znam[i] = n - i;
   }
   for(i=0; i<n; i++) {
     for (j=0; j<n; j++) {
       int nd = nod(znam[i], chis[j]);
       znam[i] /= nd;
       chis[j] /= nd;
       if (znam[i]==1) break;
     }
   }
   for(j=0; j<n; j++)  res *= chis[j];
   return res;
}
// ---------

#include <iostream>
#include <list>
using namespace std;
 
unsigned c_n_m(unsigned n, unsigned m){
if(!n || !m || (m>=n)) return 1;// шара)))
unsigned res=1;
unsigned mm=m, np1=n-m+1, nmax=n+1;
list<long> dividers;
for( ; mm>0; --mm)dividers.push_back(mm);
 
list<long>::iterator it_to_del, dividers_it=dividers.begin(), dividers_end=dividers.end();
 
for( ; np1<nmax; ++np1){
res*=np1;
dividers_it=dividers.begin();
dividers_end=dividers.end();
for( ; dividers_it != dividers_end; ++dividers_it){
    if(res% *dividers_it == 0 && *dividers_it != 1){
res/=*dividers_it;
dividers.erase(dividers_it); 
dividers_it=dividers.begin();
dividers_end=dividers.end();
 
    }
}
}
 
return res;
}
 
long c(long n, long k)
{
    if ((n==k) || (k==0))
       return 1;
    else
       return c(n-1,k-1)+c(n-1,k);
}        
 
int main(int argc, char* argv[])
{
int m, n;
n=30;
m=15;
//n=44;
//m=22;
long c_nm=c_n_m(n,m);
//long cr_nm=c(n,m);
 
cout<<c_nm<<endl;
//cout<<cr_nm<<endl;
system("pause");
return 0;
}