20170426, 16:43  #1 
"Matthew Anderson"
Dec 2010
Oregon, USA
2×11×29 Posts 
OEIS sequence A027750
Here is my mathematical doodle for today.
I noticed that OEIS.org/A027750/ has no Maple code in it. I know that there are not many Maple people who frequent mersenneforum.org, but maybe someone else will follow through with this, if I don't. The Maple code goes like this  for a from 1 to 20 do if isprime(a) = false then print(ifactor(a)); else print(a); end if; end do; This code functions, but it does not produce the same sequence as A27750. For what its worth, Matt 
20170426, 17:18  #2 
Romulan Interpreter
Jun 2011
Thailand
2×3×5^{2}×61 Posts 
you forgot to print the 1's (every number is divisible by 1) and all the composite divisors. In Pari that would be "n=100; for(i=1,n,print(divisors(i)))" or so, with a trick to avoid printing all the brackets (like adding them to a list, or so, instead of printing).

20170426, 22:19  #3  
"Forget I exist"
Jul 2009
Dumbassville
2^{6}·131 Posts 
Quote:
Code:
v=List(); for(n=1, 20, fordiv(n, d, listput(v, d))); Vec(v) 

Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Primes in A048788 OEIS Sequence  carpetpool  Miscellaneous Math  9  20170317 22:57 
OEIS  A057732 : 2^n + 3  T.Rex  Miscellaneous Math  40  20150915 14:01 
OEIS  2·A261539: (2^n+5)/3 , n even  T.Rex  Miscellaneous Math  38  20150905 16:14 
OEIS  A059242 : 2^n+5 , n odd  T.Rex  Miscellaneous Math  7  20150828 18:04 
my misunderstanding of the OEIS  science_man_88  Miscellaneous Math  11  20110518 15:04 