LibreOJ #6221. 幂数 ！（数论+dfs+剪枝）

　　写新题然后艹翻标程的感觉真是舒爽啊...

　　这题就是个dfs。。。先预处理出sqrt(n)范围内的素数，然后dfs构造合法的数就行了。

　　直接暴搜会TLE，需要剪一剪枝，不需要跑到最后一层再计算答案，边构造边更新答案，发现下一层无法出现新答案直接退出就行了

　　效率大概O(答案个数)

　　#6222也就是这题加强版需要高精度，出题人直接用python了= =   欺负C++高精度得手写...

　　而且这种写法是过不了的QAQ，加强版就不是很会了

#include<bits/stdc++.h>
#define ull unsigned long long
using namespace std;
const int maxn=;
ull n,cnt,ans;
int p;
ull prime[maxn];
bool v[maxn];
void read(ull &k)
{
    int f=;k=;char c=getchar();
    while(c<''||c>'')c=='-'&&(f=-),c=getchar();
    while(c<=''&&c>='')k=k*+c-'',c=getchar();
    k*=f;
}
void dfs(int dep,ull sum)
{
    if(dep==p)return;
    if(sum*prime[dep+]*prime[dep+]<=n)dfs(dep+,sum);
    ull now=prime[dep];
    for(int i=;;i++)
    {
        now*=prime[dep];
        if(sum*now>n)return;
        dfs(dep+,sum*now);
        cnt++;ans+=sum*now;
    }
}
int main()
{
    read(n);int m=(int)sqrt(n);v[]=;
    for(int i=;i*i<=m;i++)
    if(!v[i])for(int j=i<<;j<=m;j+=i)v[j]=;
    for(int i=;i<=m;i++)if(!v[i])prime[p++]=i;
    cnt=ans=;dfs(,);
    printf("%llu\n%llu",cnt,ans);
    return ;
}