nyoj 56阶乘因式分解(一)

来源:互联网 发布:windows nc命令下载 编辑:程序博客网 时间:2024/06/10 15:15

阶乘因式分解(一)

时间限制:3000 ms  |  内存限制:65535 KB
难度:2
描述

给定两个数m,n,其中m是一个素数。

将n(0<=n<=10000)的阶乘分解质因数,求其中有多少个m。

输入
第一行是一个整数s(0<s<=100),表示测试数据的组数
随后的s行, 每行有两个整数n,m。
输出
输出m的个数。
样例输入
2100 516 2
样例输出
2415
#include<stdio.h>int a[10000];void prim()  //计算10000以内的质数 {int i,j,k=1;for(i=2;i<10000;i++) {for(j=2;j<i;j++)if(i%j==0)break;if(j==i)        a[k++]=i;}}int f(int n, int m){int i,k=0;for(i=1;n>=a[i];)if(n%a[i]==0){if(a[i]==m||n==m)k++;n=n/a[i];continue;}elsei++;return k;}int main(){int m,n,s,i,su,count;prim();scanf("%d",&s);while(s--){su=count=0;scanf("%d%d",&n,&m);for(i=2;i<=n;i++){su=f(i,m);count+=su;}printf("%d\n",count);}return 0;}


0 0