%I #2 Mar 30 2012 17:34:20
%S -2,7,3,-26,-13,-7,96,48,24,12,-362,-181,-91,-46,-23,1350,675,337,168,
%T 84,42,-5042,-2521,-1261,-631,-316,-158,-79,18816,9408,4704,2352,1176,
%U 588,294,147,-70226,-35113,-17557,-8779,-4390,-2195,-1098,-549,-275,262086,131043,65521,32760,16380,8190,4095
%N a(n,m) =Floor[N[(-2 + Sqrt[3])^n + (-2 - Sqrt[3])^n]/2^m].
%C A triangular prime alternating sign Binet like function.
%e -2
%e 7, 3
%e -26, -13,-7
%e 96, 48, 24, 12
%e -362,-181, -91, -46, -23
%e 1350, 675, 337, 168, 84, 42
%t f[n_, m_] = N[(-2 + Sqrt[3])^n + (-2 - Sqrt[3])^n]/2^m a = Table[Table[Floor[f[n, m]], {m, 1, n}], {n, 1, 10}] aa = Flatten[a] pp = Flatten[Table[If[PrimeQ[aa[[n]]], aa[[n]], {}], {n, 1, Length[aa]}]] Length[pp]/Length[aa]
%K sign,tabl
%O 0,1
%A _Roger L. Bagula_, Apr 29 2006
|