PARI's commands

[CRGreathouse]

Originally Posted by science_man_88

oh wait k*b^b+1 would give k*1+1 for b=1 still all the primes I think.

Indeed. This the first *row* in the array calculated by my first suggestion is precisely the prime numbers, A000040. But the sequence is not its first row. It starts 2, 3, 5, 5, 13, 109, ... and so is clearly different from A000040. The second sequence is about half 1s and so is clearly distinct from the primes; it starts 1, 1, 2, 1, 1, 2, 2, 1, ... if I haven't made any mistakes. The third sequence is 5, 13, 109,163, ...

[science_man_88]

finally a question on Pari commands again:

how do you make alter and print the antidiagonals of a matrix:

[CRGreathouse]

Originally Posted by science_man_88

finally a question on Pari commands again:

how do you make alter an print the antidiagonals of a matrix:

#1216 has details. If you have a matrix T, you flatten it with vector(100,n,T[t1(n)+1, t2(n)+1]). If you have a function T, you flatten it with vector(100,n,T(t1(n)+1, t2(n)+1)). (The +1s are if you want it to be 1-based instead of 0-based; remove otherwise.) Of course you can replace 100 with whatever number you like.

[3.14159]

Originally Posted by CRGreathouse

Indeed. This the first *row* in the array calculated by my first suggestion is precisely the prime numbers, A000040. But the sequence is not its first row. It starts 2, 3, 5, 5, 13, 109, ... and so is clearly different from A000040. The second sequence is about half 1s and so is clearly distinct from the primes; it starts 1, 1, 2, 1, 1, 2, 2, 1, ... if I haven't made any mistakes. The third sequence is 5, 13, 109,163, ...

The sequence, with the restriction is:

2, 3, 5, 109, 163, 257, 271, 379, 433, 487, 541, 769, 3329, 7681, 7937, 9473, 10753, 11777, 12289, 13313, 14081, 14593, 15361, 17921, 18433, 19457, 22273, 23041, 23297, 25601, 26113, 26881, 30977, 31489, 32257, 36097, 36353, 37501, 37633, 37889, 39937, 40193, 40961, 41729, 43777, 45569, 46337, 49409, 49921, 50177, 51713, 57089, 57601, 58369, 59393, 60161, 61441, 62501, 64513, 65537, ...

[science_man_88]

Originally Posted by CRGreathouse

#1216 has details. If you have a matrix T, you flatten it with vector(100,n,T[t1(n)+1, t2(n)+1]). If you have a function T, you flatten it with vector(100,n,T(t1(n)+1, t2(n)+1)). (The +1s are if you want it to be 1-based instead of 0-based; remove otherwise.) Of course you can replace 100 with whatever number you like.

heres a better first question lol how do you originally create the matrix ? or am I missing something do you just map vectors into another vector ?

[CRGreathouse]

Originally Posted by science_man_88

heres a better first question lol how do you originally create the matrix ?

Ah. Well, Pi here has a whole thread devoted to that topic... maybe he'll post instructions there and you can help him on his project.

[CRGreathouse]

Originally Posted by 3.14159

The sequence, with the restriction is:

2, 3, 5, 109, 163, 257, 271, 379, 433, 487, 541, 769, 3329, 7681, 7937, 9473, 10753, 11777, 12289, 13313, 14081, 14593, 15361, 17921, 18433, 19457, 22273, 23041, 23297, 25601, 26113, 26881, 30977, 31489, 32257, 36097, 36353, 37501, 37633, 37889, 39937, 40193, 40961, 41729, 43777, 45569, 46337, 49409, 49921, 50177, 51713, 57089, 57601, 58369, 59393, 60161, 61441, 62501, 64513, 65537, ...

What is this, exactly?

[science_man_88]

Originally Posted by CRGreathouse

Ah. Well, Pi here has a whole thread devoted to that topic... maybe he'll post instructions there and you can help him on his project.

did he start if if not this will be hard for me to find lol.

[science_man_88]

I found how to make a matrix but I don't see how the way they make it differs from a vector.

[3.14159]

Originally Posted by CRGreathouse

What is this, exactly?

1. What I sent to the OEIS.
2. The sequence with the restriction. (k * b^b + 1, k < b^b)

[science_man_88]

I got the matrix formed and vectors to fill it but can you get it to print them as rows not just all at once ?