PE362

Squarefree Factors

对于545454,考虑他的因数分解,如果因子必须是squarefree number,则只有两种方法:
54=3∗3∗6=2∗3∗3∗3 \begin{aligned} 54 &= 3*3*6\\ &= 2*3*3*3 \end{aligned} 54=336=2333
我们记作Fst(54)=2Fst(54)=2Fst(54)=2
我们定义S(n)=∑k=2nFst(k)S(n)=\sum_{k=2}^{n}{Fst(k)}S(n)=k=2nFst(k)
求S(1010)求S(10^{10})S(1010)

思路

我们不考虑计算某一个数的FstFstFst的值,而是直接计算SSS的值
我们定义F(N,start)F(N, start)F(N,start)为一个非递降一个序列,其中都是squarefree number,且第一项≥\gestart且,序列的乘积≤\leN的序列的总数.

形式化的表达如下:
考虑如下序列的个数为F(N,start)F(N, start)F(N,start):
a1,a2,a3...ana_1, a_2, a_3...a_na1,a2,a3...an, 其中 a1≥a2≥a3...≥ana_1 \ge a_2 \ge a_3... \ge a_na1a2a3...ana1⋅a2⋅a3⋅⋅⋅⋅an≤Na_1 \cdot a_2 \cdot a_3 \cdot\cdot\cdot\cdot a_n \le Na1a2a3anN

那么最终答案表达式为F(1010,2)F(10^{10}, 2)F(1010,2)

如何计算

考虑序列的长度为1和大于的情况

  • 如果序列长度为111, 那么F(N,start)=∑i=startNF(N, start) = \sum_{i=start}^{N}F(N,start)=i=startN, iii是squarefree number
  • 如果序列长度大于111,显然a1≤Na_1 \le \sqrt{N}a1N. 考虑枚举a1a_1a1的值,注意a1a_1a1的是一个squarefree number
    F(N,start)=∑i=startNF(Ni,i)F(N, start)=\sum_{i=start}^{\sqrt{N}}F(\frac{N}{i}, i)F(N,start)=i=startNF(iN,i)
    对于是一个squarefree number的xxx,我们有μ2(x)=1\mu^2(x)=1μ2(x)=1,因此可以把式子写为
    F(N,start)=∑i=startN⋅μ2(i)+∑i=startNF(Ni,,i)⋅μ2(i)F(N, start)=\sum_{i=start}^{N}\cdot\mu^2(i) + \sum_{i=start}^{\sqrt{N}}F(\frac{N}{i,}, i)\cdot\mu^2(i)F(N,start)=i=startNμ2(i)+i=startNF(i,N,i)μ2(i)
vasp.5.4.4.18Apr17-6-g9f103f2a35 (build Feb 20 2024 19:41:19) complex executed on LinuxIFC date 2025.07.22 10:07:24 running on 48 total cores distrk: each k-point on 48 cores, 1 groups distr: one band on NCORES_PER_BAND= 4 cores, 12 groups -------------------------------------------------------------------------------------------------------- INCAR: POTCAR: PAW_PBE O 08Apr2002 POTCAR: PAW_PBE In 08Apr2002 POTCAR: PAW_PBE H 15Jun2001 POTCAR: PAW_PBE O 08Apr2002 VRHFIN =O: s2p4 LEXCH = PE EATOM = 432.3788 eV, 31.7789 Ry TITEL = PAW_PBE O 08Apr2002 LULTRA = F use ultrasoft PP ? IUNSCR = 0 unscreen: 0-lin 1-nonlin 2-no RPACOR = .000 partial core radius POMASS = 16.000; ZVAL = 6.000 mass and valenz RCORE = 1.520 outmost cutoff radius RWIGS = 1.550; RWIGS = .820 wigner-seitz radius (au A) ENMAX = 400.000; ENMIN = 300.000 eV ICORE = 2 local potential LCOR = T correct aug charges LPAW = T paw PP EAUG = 605.392 DEXC = .000 RMAX = 2.264 core radius for proj-oper RAUG = 1.300 factor for augmentation sphere RDEP = 1.550 radius for radial grids QCUT = -5.520; QGAM = 11.041 optimization parameters Description l E TYP RCUT TYP RCUT 0 .000 23 1.200 0 -.700 23 1.200 1 .000 23 1.520 1 .600 23 1.520 2 .000 7 1.500 local pseudopotential read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 4 number of lm-projection operators is LMMAX = 8 POTCAR: PAW_PBE In 08Apr2002 VRHFIN =In: s2p1 LEXCH = PE EATOM = 53.7908 eV, 3.9535 Ry TITEL = PAW_PBE In 08Apr2002 LULTRA = F use ultrasoft PP ? IUNSCR = 1 unscreen: 0-lin 1-nonlin 2-no RPACOR = 2.400 partial core radius POMASS = 114.820; ZVAL = 3.000 mass and valenz RCORE = 3.100 outmost cutoff radius RWIGS = 3.170; RWIGS = 1.677 wigner-seitz radius (au A) ENMAX = 95.934; ENMIN = 71.951 eV ICORE = 3 local potential LCOR = T correct aug charges LPAW = T paw PP EAUG = 237.125 DEXC = -.110 RMAX = 4.707 core radius for proj-oper RAUG = 1.300 factor for augmentation sphere RDEP = 3.187 radius for radial grids QCUT = -2.655; QGAM = 5.311 optimization parameters Description l E TYP RCUT TYP RCUT 0 .000 23 3.100 0 .000 23 3.100 1 .000 23 3.100 1 .000 23 3.100 2 -.100 23 3.100 3 -.100 23 3.100 local pseudopotential read in partial core-charges read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in non local Contribution for L= 2 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 5 number of lm-projection operators is LMMAX = 13 POTCAR: PAW_PBE H 15Jun2001 VRHFIN =H: ultrasoft test LEXCH = PE EATOM = 12.4884 eV, .9179 Ry TITEL = PAW_PBE H 15Jun2001 LULTRA = F use ultrasoft PP ? IUNSCR = 0 unscreen: 0-lin 1-nonlin 2-no RPACOR = .000 partial core radius POMASS = 1.000; ZVAL = 1.000 mass and valenz RCORE = 1.100 outmost cutoff radius RWIGS = .700; RWIGS = .370 wigner-seitz radius (au A) ENMAX = 250.000; ENMIN = 200.000 eV RCLOC = .701 cutoff for local pot LCOR = T correct aug charges LPAW = T paw PP EAUG = 400.000 RMAX = 2.174 core radius for proj-oper RAUG = 1.200 factor for augmentation sphere RDEP = 1.112 radius for radial grids QCUT = -5.749; QGAM = 11.498 optimization parameters Description l E TYP RCUT TYP RCUT 0 .000 23 1.100 0 .500 23 1.100 1 -.300 23 1.100 local pseudopotential read in atomic valenz-charges read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 0 read in real space projection operators read in non local Contribution for L= 1 read in real space projection operators read in PAW grid and wavefunctions read in number of l-projection operators is LMAX = 3 number of lm-projection operators is LMMAX = 5 Optimization of the real space projectors (new method) maximal supplied QI-value = 24.76 optimisation between [QCUT,QGAM] = [ 10.64, 21.54] = [ 31.73,129.89] Ry Optimized for a Real-space Cutoff 1.23 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 10.645 4.192 0.69E-04 0.93E-04 0.58E-07 0 8 10.645 8.473 0.40E-03 0.65E-03 0.17E-06 1 7 10.645 2.474 0.80E-04 0.12E-03 0.30E-07 1 7 10.645 3.912 0.24E-03 0.53E-03 0.16E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 12.09 optimisation between [QCUT,QGAM] = [ 10.76, 21.65] = [ 32.44,131.24] Ry Optimized for a Real-space Cutoff 1.65 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 11 10.764 9.989 0.15E-04 0.68E-05 0.17E-07 0 11 10.764 40.533 0.18E-03 0.28E-03 0.11E-05 1 10 10.764 8.252 0.26E-04 0.21E-04 0.71E-07 1 10 10.764 40.362 0.39E-03 0.21E-03 0.95E-06 2 10 10.764 8.820 0.14E-03 0.25E-03 0.18E-06 Optimization of the real space projectors (new method) maximal supplied QI-value = 34.20 optimisation between [QCUT,QGAM] = [ 10.60, 21.55] = [ 31.48,130.03] Ry Optimized for a Real-space Cutoff 1.23 Angstroem l n(q) QCUT max X(q) W(low)/X(q) W(high)/X(q) e(spline) 0 8 10.603 4.436 0.43E-04 0.79E-04 0.42E-07 0 8 10.603 16.131 0.35E-03 0.56E-03 0.22E-06 1 7 10.603 4.654 0.70E-04 0.28E-03 0.29E-06 PAW_PBE O 08Apr2002 : energy of atom 1 EATOM= -432.3788 kinetic energy error for atom= 0.0389 (will be added to EATOM!!) PAW_PBE In 08Apr2002 : energy of atom 2 EATOM= -53.7908 kinetic energy error for atom= 0.0000 (will be added to EATOM!!) PAW_PBE H 15Jun2001 : energy of atom 3 EATOM= -12.4884 kinetic energy error for atom= 0.0045 (will be added to EATOM!!) POSCAR: In2O3-accuracy positions in direct lattice velocities in cartesian coordinates Using predictor-coordinates on the file exchange correlation table for LEXCH = 8 RHO(1)= 0.500 N(1) = 2000 RHO(2)= 100.500 N(2) = 4000 VTST: version 3.2, (02/03/18) CHAIN: initializing optimizer OPT: Using VASP Dynamics algorithm CHAIN: Read ICHAIN 0 POSCAR: In2O3-accuracy positions in direct lattice velocities in cartesian coordinates There are predictor-coordinates on the file. we can't use them due to change of POTIM! -------------------------------------------------------------------------------------------------------- ion position nearest neighbor table 1 0.758 0.297 0.233- 134 2.20 97 2.26 113 2.28 147 2.31 2 0.496 0.716 0.569- 161 0.98 148 2.21 3 0.410 0.621 0.113- 131 2.15 119 2.18 145 2.28 105 2.30 4 0.096 0.964 0.347- 132 2.10 106 2.21 120 2.28 146 2.34 5 0.503 0.759 0.211- 131 2.18 150 2.20 99 2.22 115 2.27 6 0.221 0.114 0.443- 163 0.98 100 2.27 132 2.29 116 2.42 7 0.168 0.089 0.089- 133 2.15 117 2.25 107 2.33 151 2.43 8 0.842 0.439 0.337- 118 2.23 134 2.28 108 2.29 152 2.49 9 0.152 0.373 0.104- 101 2.18 117 2.21 135 2.21 145 2.30 10 0.828 0.707 0.337- 102 2.18 136 2.25 118 2.27 146 2.30 根据上述INCAR运行出的OUTCAR内容来看,参数设置是否合理
07-23
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值