mersenneforum.org OEIS sequence A027750
 User Name Remember Me? Password
 Register FAQ Search Today's Posts Mark Forums Read

2017-04-26, 16:43   #1
MattcAnderson

"Matthew Anderson"
Dec 2010
Oregon, USA

58810 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
Attached Files
 A027750.pdf (79.0 KB, 46 views)

 2017-04-26, 17:18 #2 LaurV Romulan Interpreter     Jun 2011 Thailand 2×17×251 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).
2017-04-26, 22:19   #3
science_man_88

"Forget I exist"
Jul 2009
Dumbassville

26×131 Posts

Quote:
 Originally Posted by LaurV 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).
Code:
v=List(); for(n=1, 20, fordiv(n, d, listput(v, d))); Vec(v)
is already listed in the sequence

 Similar Threads Thread Thread Starter Forum Replies Last Post carpetpool Miscellaneous Math 9 2017-03-17 22:57 T.Rex Miscellaneous Math 40 2015-09-15 14:01 T.Rex Miscellaneous Math 38 2015-09-05 16:14 T.Rex Miscellaneous Math 7 2015-08-28 18:04 science_man_88 Miscellaneous Math 11 2011-05-18 15:04

All times are UTC. The time now is 07:33.

Thu May 28 07:33:12 UTC 2020 up 64 days, 5:06, 1 user, load averages: 2.86, 2.55, 2.20