FOJ 1851 分解素数因子

Problem 1851 组合数

Accept: 121    Submit: 409
Time Limit: 1000 mSec    Memory Limit : 32768 KB

Problem Description

组合数 C(N,M)=N!/(M!*(N-M)!). 问题是求 C(N,M)中不同的素因子的个数

Input

输入包含多组数据(数据组数<= 10)，每组数据只有一行 2个正整数表示N和M (1 <= M <= N <= 1000000)

N=M=0表示数据结束.

Output

输出一行答案,表示C(N,M)中不同素因子个数.

Sample Input

7 4 0 0

Sample Output

2

Hint

C(7,4)=35=5*7,所以有2种不同的素因子.

Source

FOJ有奖月赛-2009年10月——稚鹰翱翔

 

#include<stdio.h> int a[78500]; int b[78500]; int prime[78500]; int pp[1000050]={0}; int max(int n,int m) { return m>n-m ?m:n; } void fun() { int i,j,k; for(i=2,k=0;i<1000050;i++) { if(!pp[i]) { for(j=2*i;j<1000050;j+=i) pp[j]=1; prime[k++]=i; } } } void merge(int n,int *c) { int i; int tt=n; for(i=0;n>prime[i];i++) { while(tt) { c[i]+=(tt/=prime[i]); } tt=n; } } int main() { int n,m; int i; int temp; int ans; fun(); while(scanf("%d%d",&n,&m),m+n) { for(i=0;i<78498;i++) { a[i]=b[i]=0; } merge(n,a); merge(m,b); merge(n-m,b); temp=max(n,m); for(i=0;temp>prime[i];i++) { a[i]-=b[i]; } for(i=ans=0;n>prime[i];i++) if(a[i]) ans++; printf("%d/n",ans); } return 0; }