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

2721: [Violet 5]樱花|约数个数

时间:2016-03-20 08:10:25      阅读:197      评论:0      收藏:0      [点我收藏+]

标签:

先跪一发题目背景QAQ
显然x,y>n!,然后可以设y=n!+d
原式子可以化简成

x=n!2d+n!

那么解的个数也就是n!的因子个数,然后线性筛随便搞一搞

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
#define N 1000008
#define mod  1000000007
using namespace std;
int sc()
{
    int i=0,f=1; char c=getchar();
    while(c>‘9‘||c<‘0‘){if(c==‘-‘)f=-1;c=getchar();}
    while(c>=‘0‘&&c<=‘9‘)i=i*10+c-‘0‘,c=getchar();
    return i*f;
}
long long ans=1;
int lo[N],low[N],a[N],prime[N],s[N],top;
int sum[N],n;
void pre(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!a[i])
            s[prime[++top]=low[i]=lo[i]=i]=1;
        for(int j=1;prime[j]*i<=n;j++)
        {
            a[i*prime[j]]=1;
            lo[i*prime[j]]=prime[j];
            if(i%prime[j]==0)
            {
                low[i*prime[j]]=low[i]*prime[j];
                s[i*prime[j]]=s[i]+1;
                break;
            }
            low[i*prime[j]]=prime[j];
            s[i*prime[j]]=1;
        }
    }
}
int main()
{
    pre(n=sc());
    for(int i=2;i<=n;i++)
    {
        int now=i;
        while(now!=1)
            sum[lo[now]]+=2*s[now],now/=low[now];
    }
    for(int i=1;i<=n;i++)
        ans=ans*(sum[i]+1)%mod;
    cout<<ans;
    return 0;
}

2721: [Violet 5]樱花|约数个数

标签:

原文地址:http://blog.csdn.net/ws_yzy/article/details/50935095

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