回复：欧拉计划的解密计划的交流计划

214题，快了一点。
Euler[1]=1;Euler[n_]:=Euler[n]=If[Euler[EulerPhi@n]>=24,Infinity,Euler[EulerPhi@n]+1];STP=Select[Table[Prime[i]-1,{i,636606,PrimePi[40000000]}],Euler[#]==24&];Length@STP+Total@STP

用CompilationTarget产生C代码，运行快，但是为什么我运行虚拟全书的例子，它输出 一个库函数不能从编译函数中产生，你们也是这样吗

@草红样
C编译器怎么自己安装，欧拉计划后面的题目如不产生C代码，很多题估计一分钟内都很搞不定，比如：221题
cf=Compile[{},Total@Select[Range[64000000],FractionalPart@Sqrt[DivisorSigma[2.,#]]==0&],CompilationOptions->{"InlineExternalDefinitions"->True},RuntimeOptions->"Speed"];cf[]//AbsoluteTiming
I

还需要以下的设置（注意要改成自己的路径），不知道重启Mathematica后是否会失效，不想每次都运行的话，
可以写到init.m里面（可以执行 SystemOpen[DirectoryName@FindFile["init.m"]] 来打开）
Needs["CCompilerDriver`"]
$CCompiler = {"Compiler" ->
CCompilerDriver`MinGWCompiler`MinGWCompiler,
"CompilerInstallation" -> "C:/MinGW64",
"CompilerName" -> Automatic};I

203题，又找到个简单的。
Total[Select[DeleteDuplicates[Flatten[Table[Binomial[n,k],{n,0,50},{k,0,n}]]],Times@@(Transpose[FactorInteger[#]])[[2]]==1&]]

第51题：
Clear["`*"];
f[x_] := Module[{m = (x - 1)/2},
Count[Join[7 - 10 #1 + 4 #1^2 & /@ Range[2, m + 1],
3 - 6 #1 + 4 #1^2 & /@ Range[2, m + 1],
5 - 8 #1 + 4 #1^2 & /@ Range[2, m + 1]], _?PrimeQ]/(4 m + 1)];
f[26589] // N
f[26591] // N得到0.100006和0.0999981
为什么26591错的？好纠结啊I

179题：
TD=Table[DivisorSigma[0,n],{n,2,10^7}];Count[TD-RotateRight@TD,0]

138题：
Reduce[(b/2)^2+h^2==L^2&&(h==b+1||h==b-1)&&b>1&&L>0,Integers];
Expand@Total@Table[%[[1,4,2]]+%[[2,4,2]]/.C[1]->i,{i,6}]

51题：
k=1;While[!Or@@Table[Block[{ID=IntegerDigits@Prime@k},If[MemberQ[ID,j],Count[PrimeQ@FromDigits@#&/@Table[ID/.{j->i},{i,j+1,9}],True]>6,False]],{j,0,2}],k++];Prime@k

140题：
Reduce[Solve[Sum[RSolve[{f[n]==f[n-1]+f[n-2],f[1]==1&&f[2]==4},f[n],n][[1,1,2]]x^n,{n,Infinity}]==k,x][[1,1,2,3,3,2,1]]==n^2&&k>0&&n>0,Integers];
Expand@Total@Table[(Total[Table[%[[i,3,2]],{i,4}]]/.C[1]->j)+(Total[Table[%[[i,3,2]],{i,5,6}]]/.C[1]->j+1),{j,0,4}]

137题：
Expand[Reduce[Solve[Sum[Fibonacci[n]x^n,{n,Infinity}]==k,x][[1,1,2,3,3,1]]==n^2&&k>0&&n>0,Integers][[2,3,2]]/.C[1]->5]

太厉害了，支持，学习了

76题：
SeriesCoefficient[1/QPochhammer[x], {x, 0, 100}] - 1

80题：
Total[Total@RealDigits[Sqrt[#],10,100][[1]]&/@Complement[Range[100],Range[10]^2]]