https://ac.nowcoder.com/acm/contest/9982/I
https://codeforces.com/problemset/problem/546/D
思路:总的思路就是要o(n)拿出每个数的最大/最小质因数。然后就可以logn处理。
欧拉筛本身的处理就是可以把v[i]存成i的最小质因数来筛,详见F. 1.小W 的质数(prime)[欧拉筛再理解]
如果我要处理出每个数的最大质因数呢?类似一个dp的思路
Maxp[x]:表示x的最大质因子
枚举的i是x的最大因子
Maxp[i*primes[j]]=max(primes[j],Maxp[i]);不断更新能整除的小的质因子得出最大的/
n=read();
for(int i=2;i<=10000000;++i){
if(!vis[i]){p[++cnt]=i;Maxp[i]=i;}
for(int j=1;p[j]*i<=10000000;++j){
vis[p[j]*i]=1;//线性筛
Maxp[p[j]*i]=max(p[j],Maxp[i]);//求最大的质因子
if(i%p[j]==0)break;
}
}
for(int i=1;i<=n;++i){
a[i]=read();
int tmp=a[i];
while(tmp!=1){
int t=Maxp[tmp];
while(tmp%t==0)tmp/=t,printf("%d ",t);
}
puts("");
}
#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<cstdio>
#include<algorithm>
#define debug(a) cout<<#a<<"="<<a<<endl;
using namespace std;
const int maxn=4e6+100;
typedef int LL;
const int mod=1e9+7;
inline LL read(){LL x=0,f=1;char ch=getchar(); while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;}
LL primes[maxn];
LL v[maxn];
LL m=0;
LL ksm(LL a,LL b){
LL ans=1;
while(b){
if(b&1)ans=ans*a%mod;
a=a*a%mod;
b/=2;
}
return ans%mod;
}
void getprimes(int n)
{
for(LL i=2;i<=n;i++)
{
if(v[i]==0) {v[i]=i;primes[++m]=i;}
for(LL j=1;j<=m;j++)
{
if(primes[j]>v[i]||primes[j]>n/i) break;
v[i*primes[j]]=primes[j];
}
}
}
LL get(LL x){
LL sum=0;
while(x>0) x/=10,sum++;
return sum;
}
LL solve(LL x){
long long sum=0;
if(x==1) return 0;
LL tmp=x;
while(tmp!=1){
LL t=v[tmp];
while(tmp%t==0){
tmp/=t;
sum=(sum*ksm(10,get(t))%mod+t%mod)%mod;
}
}
if(tmp>1) sum=(sum*ksm(10,get(tmp))%mod+tmp%mod)%mod;
return sum%mod;
}
int main(void)
{
getprimes(4000000);
LL n;n=read();
long long ans=0;
for(LL i=1;i<=n;i++){
ans=(ans%mod+solve(i)%mod)%mod;
}
printf("%lld\n",ans);
return 0;
}