码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj2721: [Violet 5]樱花

时间:2018-09-06 14:30:55      阅读:154      评论:0      收藏:0      [点我收藏+]

标签:个数   sizeof   ons   ++   namespace   typedef   get   ret   als   

本来以为这题比较难的。。。

推下柿子

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

(x+y)*n!=xy

设y=n!+d

(x+n!+d)*n!=x(n!+d)

n!^2+d*n!=x*d

(n!^2+d*n!)/d=x 只要令n!^2/d是整数就行了

那么就是算n!的约数个数

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const LL mod=1e9+7;

int pr,prime[1100000];
bool v[1100000];
void getprime(int n)
{
    pr=0;
    memset(v,true,sizeof(v));
    for(int i=2;i<=n;i++)
    {
        if(v[i]==true)prime[++pr]=i;
        for(int j=1;j<=pr&&i*prime[j]<=n;j++)
        {
            v[i*prime[j]]=false;
            if(i%prime[j]==0)break;
        }
    }
}


int main()
{
    int n;
    scanf("%d",&n);
    getprime(n);
    
    LL ans=1;
    for(int i=1;i<=pr;i++)
    {
        LL c=0;
        for(LL j=prime[i];j<=n;j*=prime[i])
            c=(c+n/j)%mod;
            
        ans=(ans*(c*2+1))%mod;
    }
    printf("%lld\n",ans);
    return 0;
}

 

bzoj2721: [Violet 5]樱花

标签:个数   sizeof   ons   ++   namespace   typedef   get   ret   als   

原文地址:https://www.cnblogs.com/AKCqhzdy/p/9597736.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!