poj1306 Combinations

//poj 1306 Combinations 简单题
/*
本来想用高精度，后来看数据好小，就直接模拟平时的算法，
把上面跟下面大的那一块先约掉不管，然后再算剩余的，这里
用了素数分解，可能还有其他方法，不管了，简单题能水就行。
*/
#include <iostream>
#include <algorithm>
using namespace std;
const int inf = 1<<28;

int n,m;
int prime[100],isprime[100],k;
void make_prime()
{
     for (int i=0;i<100;i++) isprime[i]=1;
     k=0;
     for (int i=2;i<100;i++)
     if (isprime[i])
     {
                    prime[++k]=i;
                    for (int j=i+i;j<100;j+=i) isprime[j]=0;
     }
}
int p[2][100];

void doit(int n,int j)
{
     for (int i=1;i<=k;i++)
     {
         if (n==1) break;
         while (n%prime[i]==0)
         {
               n/=prime[i];
               p[j][i]++;
         }
     }
}
int main()
{
    make_prime();
    while (scanf("%d%d",&n,&m) && n+m)
    {
          memset(p,0,sizeof(p));
          for (int i=n;i>max(m,n-m);i--)
          {
              doit(i,0);
          }
          for (int i=min(m,n-m);i>=1;i--)
          {
              doit(i,1);
          }
         // for (int i=1;i<=k;i++)
         // printf("%d:%d %d/n",prime[i],p[0][i],p[1][i]);
          int ans=1;
          for (int i=1;i<=k;i++)
          {
             for (int j=0;j<p[0][i]-p[1][i];j++)
             ans*=prime[i];
          }
         
          printf("%d things taken %d at a time is %d exactly./n",n,m,ans);
         
         
    }
    system("pause");
    return 0;
}