cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A288097 Square array read by antidiagonals downwards: A(n, 1) = smallest base-prime(n) Wieferich prime and A(n, k) = smallest base-A(n, k-1) Wieferich prime for k > 1.

This page as a plain text file.
%I A288097 #9 Jan 21 2023 13:21:27
%S A288097 1093,2,11,1093,71,2,2,3,1093,5,1093,11,2,2,71,2,71,1093,1093,3,2,
%T A288097 1093,3,2,2,11,1093,2,2,11,1093,1093,71,2,1093,3,1093,71,2,2,3,1093,2,
%U A288097 11,13,2,3,1093,1093,11,2,1093,71,2,2,1093,11,2,2,71,1093,2,3,1093,1093,7
%N A288097 Square array read by antidiagonals downwards: A(n, 1) = smallest base-prime(n) Wieferich prime and A(n, k) = smallest base-A(n, k-1) Wieferich prime for k > 1.
%e A288097 Array starts
%e A288097 1093,    2, 1093,    2, 1093,    2, 1093,    2, 1093,    2
%e A288097   11,   71,    3,   11,   71,    3,   11,   71,    3,   11
%e A288097    2, 1093,    2, 1093,    2, 1093,    2, 1093,    2, 1093
%e A288097    5,    2, 1093,    2, 1093,    2, 1093,    2, 1093,    2
%e A288097   71,    3,   11,   71,    3,   11,   71,    3,   11,   71
%e A288097    2, 1093,    2, 1093,    2, 1093,    2, 1093,    2, 1093
%e A288097    2, 1093,    2, 1093,    2, 1093,    2, 1093,    2, 1093
%e A288097    3,   11,   71,    3,   11,   71,    3,   11,   71,    3
%e A288097   13,    2, 1093,    2, 1093,    2, 1093,    2, 1093,    2
%e A288097    2, 1093,    2, 1093,    2, 1093,    2, 1093,    2, 1093
%t A288097 f[n_] := Block[{p = 2}, While[! Divisible[n^(p - 1) - 1, p^2], p = NextPrime@ p]; p]; T[n_, k_] := T[n, k] = If[k == 1, f@ Prime@ n, f@ T[n, k - 1]]; Table[Function[n, T[n, k]][m - k + 1], {m, 12}, {k, m, 1, -1}] // Flatten (* _Michael De Vlieger_, Jun 06 2017 *)
%o A288097 (PARI) a039951(n) = forprime(p=1, , if(Mod(n, p^2)^(p-1)==1, return(p)))
%o A288097 table(rows, cols) = forprime(p=1, prime(rows), my(i=0, w=a039951(p)); while(i < cols, print1(w, ", "); w=a039951(w); i++); print(""))
%o A288097 table(10, 10) \\ print initial 10 rows and 10 columns of table
%Y A288097 Cf. A039951, A252801, A252802, A252812, A269111, A281001, A281002, A268479.
%K A288097 nonn,tabl
%O A288097 1,1
%A A288097 _Felix Fröhlich_, Jun 05 2017
%E A288097 More terms from _Michael De Vlieger_, Jun 06 2017