PKU 2992 【质因数分解】
2010-04-09 13:28
351 查看
利用勒让德定理求n!因数个数,组合数公式C(n,k) =n! /((n-k)!*k!)
首先打素数表。
代码#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int N = 431+10;
bool isp
;
int p
;
int pcnt;
void init()
{
memset(isp,true,sizeof(isp));
pcnt = 0;
for(int i=2;i<N;i++)
{
if(isp[i])
{
p[pcnt++] = i;
for(int j=i+i;j<N;j+=i)
{
isp[j] = false;
}
}
}
}
int divisor(int n,int p)
{
int ans = 0;
while(n>=p)
{
ans += n/p;
n /= p;
}
return ans;
}
int main()
{
init();
int k,n;
while(scanf("%d%d",&n,&k)!=EOF)
{
int t;
long long ans=1;
for(int i=0;i<pcnt;i++)
{
t = divisor(n,p[i]) - divisor(n-k,p[i]) - divisor(k,p[i]);
ans *= (t+1);
}
printf("%lld\n",ans);
}
return 0;
}
首先打素数表。
代码#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int N = 431+10;
bool isp
;
int p
;
int pcnt;
void init()
{
memset(isp,true,sizeof(isp));
pcnt = 0;
for(int i=2;i<N;i++)
{
if(isp[i])
{
p[pcnt++] = i;
for(int j=i+i;j<N;j+=i)
{
isp[j] = false;
}
}
}
}
int divisor(int n,int p)
{
int ans = 0;
while(n>=p)
{
ans += n/p;
n /= p;
}
return ans;
}
int main()
{
init();
int k,n;
while(scanf("%d%d",&n,&k)!=EOF)
{
int t;
long long ans=1;
for(int i=0;i<pcnt;i++)
{
t = divisor(n,p[i]) - divisor(n-k,p[i]) - divisor(k,p[i]);
ans *= (t+1);
}
printf("%lld\n",ans);
}
return 0;
}
相关文章推荐
- PKU2992
- pku 2992 因子和因子个数的应用
- pku2992 Divisors
- pku2992 Divisors
- pku2992(约数,素因子分解)
- Divisors pku 2992
- Pku3673
- PKU1287裸的Kruskal算法低空飘过
- dp_Pku1887
- 1211: [HNOI2004]树的计数 (prufer编码,排列组合,质因数分解)
- PKU2456二分查找
- PKU的Bridging Signals
- PKU 2828 Buy Tickets 线段树想法题 (第4天)
- 假金币问题-PKUacm1029-ACM
- poj 2992 Divisors 整数分解
- DP PKU 1754
- PKU 3630 - Phone List(Trie)
- Bellman算法PKU 3259
- PKU 2226 - Muddy Fields (最小点覆盖)
- HDU 1695 GCD 欧拉函数+容斥原理+质因数分解