bzoj 2721

题解：首先推一发式子（见csdn https://blog.csdn.net/lleozhang/article/details/83415995）

因为x是整数，所以x的数量显然为能使取得整数的t的个数，也就是求的约数个数

而根据约数个数和公式(设一个数)

可以将前n个数质因子分解，然后将质因子的幂次相乘，最后将所有幂次*2+1后乘在一起即可。

#include <cstdio>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>
#define ll long long
#define mode 1000000007
#define maxn 1000000
using namespace std;
ll fac[];
ll pri[];
bool used[];
int tot=;
int n;
void init()
{
    scanf("%d",&n);
    for(int i=;(ll)i*i<=n;i++)
    {
        if(!used[i])
        {
            pri[++tot]=i;
        }
        for(int j=;j<=tot&&(ll)i*pri[j]<=n;j++)
        {
            used[i*pri[j]]=;
            if(i%pri[j]==)
            {
                break;
            }
        }
    }
    for(int i=;i<=n;i++)
    {
        int t=i;
        for(int j=;(ll)pri[j]*pri[j]<=t&&j<=tot;j++)
        {
            while(t%pri[j]==)
            {
                fac[pri[j]]++;
                t/=pri[j];
                if(fac[pri[j]]>=mode)
                {
                    fac[pri[j]]-=mode;
                }
            }
        }
        if(t!=)
        {
            fac[t]++;
            if(fac[t]>=mode)
            {
                fac[t]-=mode;
            }
        }
    }
    ll ans=;
    for(int i=;i<=n;i++)
    {
        fac[i]<<=;
        fac[i]%=mode;
        ans*=(fac[i]+);
        ans%=mode;
    }
    printf("%lld\n",ans);
}
int main()
{
    init();
    return ;
}