求一个数的因子个数

先将n质因数分解后,如n=x1^a*x2^b*x3^c

个数为 (a+1)*(b+1)*(c+1)

poj 2992

#include<stdio.h>
#include<string.h>
int p[90];
int v[432];
int pn;
int e[432][90];

void pri()
{
  int i,j;
  pn=0;
  memset(v,0,sizeof(v));
  for(i=2;i<=431;i++)
  {
    if(v[i]==0)
    {
      p[pn++]=i;
      for(j=i;j<=431;j+=i)
        v[j]=1;
    }
  }
}

void fac()
{
  int i,j,k;
  for(i=2;i<=431;i++)
  {
    k=i;
    for(j=0;j<pn;j++)
      while(k>1&&k%p[j]==0) //质因数分解
      {
        e[i][j]++;
        k/=p[j];
      }
  }
  for(i=3;i<=431;i++)
    for(j=0;j<pn;j++)
      e[i][j]+=e[i-1][j];
}

int main()
{
  //freopen("FactorNumber.in","r",stdin);
  //freopen("FactorNumber.out","w",stdout);
  int i,n,k;
  long long sum;
  pri();
  fac();
  while(scanf("%d%d",&n,&k)!=EOF)
  {
    sum=1;
    for(i=0;i<pn;i++)
      sum*=(1+e[n][i]-e[n-k][i]-e[k][i]);
    printf("%lld\n",sum);
  }
  return 0;
}

转载于:https://www.cnblogs.com/lmjer/p/8137471.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值