bzoj2721樱花——质因数分解

Posted narh

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj2721樱花——质因数分解相关的知识,希望对你有一定的参考价值。

题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2721

要推式子!

发现x和y一定都比 n! 大。不妨设 x = n!+k;

  则1/x + 1/y = 1/ n!

  <=> ( n! + k + y ) / ( n! + k ) * y = 1 / n!

  <=> n! * y+ k * y= (n!)^2 + n! * k + n! * y

  <=> y = (n!)^2 / k + n!

所以( x, y ) 的个数 =  y 的个数 = k的个数 = (n!)^2 的约数个数。

通过质因数分解求出约数个数。

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
const int N=1e6+5;
const ll mod=1e9+7;
int n,pri[N],mindiv[N],cnt[N],xnt;
ll ans;
void init()
{
    ans=1;
    for(int i=2;i<=n;i++)
    {
        if(!mindiv[i])pri[++xnt]=i,mindiv[i]=xnt;
        for(int j=1;j<=xnt&&i*pri[j]<=n;j++)// i*pri[j]<=n!!!
        {
            mindiv[i*pri[j]]=j;
            if((i%pri[j])==0)break;
        }
    }
}
ll calc(int a)
{
    while(a!=1)
    {
        cnt[mindiv[a]]++;
        a/=pri[mindiv[a]];
    }
}
int main()
{
    scanf("%d",&n);
    init();
    while(n)calc(n),n--;
    for(int i=1;i<=xnt;i++)
        (ans*=(cnt[i]*2+1))%=mod;
    printf("%lld",ans);
    return 0;
}

 

以上是关于bzoj2721樱花——质因数分解的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ_2721_[Violet 5]樱花_数学

P1445 [Violet]樱花(质因数分解)

BZOJ 2721 樱花

bzoj 2721[Violet 5]樱花 数论

bzoj2721[Violet 5]樱花

bzoj2721: [Violet 5]樱花