|
%I #16 May 03 2019 22:31:18
%S 1,3,4,5,7,8,9,11,12,13,15,16,17,19,20,21,23,24,25,27,28,29,30,31,32,
%T 33,35,36,37,39,40,41,43,44,45,47,48,49,51,52,53,55,56,57,59,60,61,63,
%U 64,65,67,68,69,70,71,72,73,75,76,77,79,80,81,83,84,85,87,88,89,90,91,92,93,95,96,97,99,100
%N 1 together with the numbers with property that at least two successive of their divisors are odd or even.
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Table_of_divisors">Table of divisors</a>
%e 12 is in the sequence because its divisors are 1, 2, 3, 4, 6, 12 and then we can see that the successive divisors 4, 6, and 12 are even numbers.
%p isA183147 := proc(n)
%p local di,i ;
%p if n = 1 then
%p return true;
%p end if;
%p di := sort(convert(numtheory[divisors](n),list)) ;
%p for i from 1 to nops(di)-1 do
%p if (op(i+1,di)-op(i,di)) mod 2 = 0 then
%p return true;
%p end if;
%p end do;
%p return false;
%p end proc:
%p for n from 1 to 100 do
%p if isA183147(n) then
%p printf("%d,",n);
%p end if;
%p end do; # _R. J. Mathar_, Sep 30 2011
%Y Complement of A182991, for n > 1.
%Y Cf. A182996.
%K nonn,easy
%O 1,2
%A _Omar E. Pol_, Aug 04 2011
|
