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; }