传送门
题
意
:
求
∑
i
=
2
n
F
(
i
)
,
其
中
F
(
i
)
表
示
将
i
的
每
一
个
质
因
数
从
大
到
小
拼
接
成
题意:求\sum_{i=2}^{n}F(i),其中F(i)表示将i的每一个质因数从大到小拼接成
题意:求∑i=2nF(i),其中F(i)表示将i的每一个质因数从大到小拼接成
的
数
,
如
12
=
2
∗
2
∗
3
,
那
么
F
(
12
)
=
223
的数,如12=2*2*3,那么F(12)=223
的数,如12=2∗2∗3,那么F(12)=223
分
析
:
欧
拉
筛
的
妙
用
,
如
果
i
为
素
数
,
那
么
F
(
i
)
=
i
.
分析:欧拉筛的妙用,如果i为素数,那么F(i)=i.
分析:欧拉筛的妙用,如果i为素数,那么F(i)=i.
在
分
解
素
数
的
同
时
可
以
用
d
p
来
预
处
理
出
每
一
个
数
在分解素数的同时可以用dp来预处理出每一个数
在分解素数的同时可以用dp来预处理出每一个数
关
键
在
于
,
状
态
是
怎
么
转
移
的
:
关键在于,状态是怎么转移的:
关键在于,状态是怎么转移的:
d
p
[
j
]
=
d
p
[
j
/
i
]
∗
q
p
o
w
(
10
,
c
n
t
)
,
其
中
j
为
合
数
,
i
为
j
的
dp[j]=dp[j/i]*qpow(10,cnt),其中j为合数,i为j的
dp[j]=dp[j/i]∗qpow(10,cnt),其中j为合数,i为j的
其
中
一
个
素
因
子
,
c
n
t
代
表
i
为
多
少
位
数
其中一个素因子,cnt代表i为多少位数
其中一个素因子,cnt代表i为多少位数
对
于
12
=
2
∗
2
∗
3
,
当
素
因
子
为
2
时
,
可
以
知
道
:
对于12=2*2*3,当素因子为2时,可以知道:
对于12=2∗2∗3,当素因子为2时,可以知道:
d
p
[
2
]
=
2
,
dp[2]=2,
dp[2]=2,
d
p
[
4
]
=
d
p
[
4
/
2
]
∗
10
+
2
=
22
,
dp[4]=dp[4/2]*10+2=22,
dp[4]=dp[4/2]∗10+2=22,
d
p
[
6
]
=
d
p
[
6
/
2
]
∗
10
+
2
=
2
dp[6]=dp[6/2]*10+2=2
dp[6]=dp[6/2]∗10+2=2
.
.
.
.
....
....
d
p
[
12
]
=
d
p
[
12
/
2
]
∗
10
+
2
=
22
dp[12]=dp[12/2]*10+2=22
dp[12]=dp[12/2]∗10+2=22
当
素
因
子
为
3
时
,
可
以
知
道
:
当素因子为3时,可以知道:
当素因子为3时,可以知道:
d
p
[
3
]
=
3
dp[3]=3
dp[3]=3
d
p
[
6
]
=
d
p
[
6
/
3
]
∗
10
+
3
=
23
,
dp[6]=dp[6/3]*10+3=23,
dp[6]=dp[6/3]∗10+3=23,
d
p
[
12
]
=
d
p
[
12
/
3
]
∗
10
+
3
=
223
dp[12]=dp[12/3]*10+3=223
dp[12]=dp[12/3]∗10+3=223
.
.
.
.
.
.
......
......
#include <set>
#include <map>
#include <list>
#include <cmath>
#include <stack>
#include <queue>
#include <string>
#include <bitset>
#include <vector>
#include<cstring>
#include <stdio.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define lson rt<<1,l,mid
#define rson rt<<1|1,mid+1,r
#define lowbit(x) ((x)&-(x))
typedef unsigned long long ull;
#define freedom() ios::sync_with_stdio(false);
inline int read(){int s=0,w=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();return s*w;}
const int maxn = 4e6+5;
int pri[maxn+5],tot;
bool vis[maxn+5];
ll mod=1e9+7,dp[maxn+5];
ll qpow(ll a,ll b) // 快速幂
{
ll ans=1;
while (b){
if(b&1) ans=ans*a%mod;
a=a*a%mod;
b>>=1;
}
return ans;
}
void init()
{
for(int i=2;i<maxn;i++)
if(!vis[i])
{
pri[++tot]=i;
dp[i]=i;
for(int j=2*i;j<maxn;j+=i)
{
vis[j]=true;
int k=i,cnt=0;
while (k) k/=10,cnt++;
dp[j]=(dp[j/i]*qpow(10,cnt)%mod+(ll)i)%mod;
}
}
}
int main()
{
int n=read();
ll res=0; init();
for(int i=2;i<=n;i++)
res=(res+dp[i])%mod;
printf("%lld\n",res);
return 0;
}