|
%I #21 Apr 20 2019 11:52:04
%S 6,1,3,4,7,5,2,6,9,2,0,2,2,3,4,4,1,6,0,1,8,0,4,1,6,6,3,8
%N Decimal expansion of estimate of the strongly triple-free set constant.
%H Steven R. Finch, <a href="/FinchTriple.html">Triple-Free Sets of Integers</a> [From Steven Finch, Apr 20 2019]
%H Julien Cassaigne and Paul Zimmermann, <a href="http://iml.univ-mrs.fr/~cassaign/publis/cz.ps.gz">Numerical Evaluation of the Strongly Triple-Free Constant</a> (1996). [From _Steven Finch_, Feb 25 2009]
%H Julien Cassaigne and Paul Zimmermann, <a href="https://oeis.org/A086316/a086316_3.pdf">Numerical Evaluation of the Strongly Triple-Free Constant</a> (pdf file, 1996).
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Triple-FreeSet.html">Triple-Free Set</a>
%e 0.613475269...
%e 0.6134752692022344160180416638... - _Steven Finch_, Feb 25 2009
%t f[k_,n_]:=1+Floor[FullSimplify[Log[n/3^k]/Log[2]]]; g[n_]:=Floor[FullSimplify[Log[n]/Log[3]]]; peven[n_]:=Sum[Quotient[f[k,n]+Mod[k+1,2],2],{k,0,g[n]}]; podd[n_]:=Sum[Quotient[f[k,n]+Mod[k,2],2],{k,0,g[n]}]; p[n_]:=Max[peven[n],podd[n]]; v[1]=1;j=1;k=1;n=4001; For[k=2,k=n,k++,If[2*v[k-j]<3^j,v[k]=2*v[k-j],{v[k]=3^j,j++}]]; Sum[p[v[n]]*(1/v[n]-1/v[n+1]),{n,1,4000}]/3 (* _Steven Finch_, Feb 25 2009 *)
%Y Cf. A050295, A050296.
%K nonn,cons,more
%O 0,1
%A _Eric W. Weisstein_, Jul 15 2003
%E More terms from _Steven Finch_, Feb 25 2009
|
