|
%I #12 Nov 02 2016 12:54:46
%S 0,6,37,43,74,87,118,155,186,192,199,230,236,267,280,304,311,317,348,
%T 385,392,416,429,460,466,497,504,510,541,578,622,659,690,696,703,734,
%U 740,771,784,808,815,852,889,896,920,933,964,970,1001,1008,1014,1045,1082,1126,1163,1194,1200,1207,1238,1244,1275,1288,1312,1319,1356,1387,1393,1400,1424
%N Intersection of A277722 and A277723.
%C See A277728 for discussion.
%p Digits := 120;
%p isA277722 := proc(n)
%p a276800 := 3.3829757679062374941227085364550345869493820437485761820195626772353718960099402922235933340043661396041006 ;
%p for x from floor((n-3)/a276800) to (n+3)/a276800 do
%p if floor(x*a276800) = n then
%p return true;
%p end if;
%p end do:
%p return false;
%p end proc:
%p isA277723 := proc(n)
%p a276801 := 6.2222625231203986266745611011083211873735607898461684287983213166395751180919067179620287534326731537460804;
%p for x from floor((n-3)/a276801) to (n+3)/a276801 do
%p if floor(x*a276801) = n then
%p return true;
%p end if;
%p end do:
%p return false;
%p end proc:
%p isA277726 := proc(n)
%p isA277722(n) and isA277723(n) ;
%p end proc:
%p for n from 0 to 8000 do
%p if isA277726(n) then
%p printf("%d,",n) ;
%p end if;
%p end do: # _R. J. Mathar_, Nov 02 2016
%Y Cf. A003144, A003145, A003146, A058265, A158919, A275926, A276799, A276800, A277722, A277723, A277724, A277725, A277726, A277727, A277728.
%K nonn
%O 1,2
%A _N. J. A. Sloane_, Oct 30 2016
%E Corrected by _R. J. Mathar_, Nov 01 2016
|
