login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A181461 Numbers n such that 53 is the largest prime factor of n^2-1. 3
52, 54, 105, 107, 160, 211, 319, 370, 476, 529, 531, 584, 637, 741, 743, 847, 849, 900, 902, 953, 1059, 1220, 1273, 1324, 1377, 1379, 1483, 1538, 1644, 1695, 1801, 1803, 2015, 2174, 2278, 2386, 2437, 2543, 2651, 2755, 2861, 2969, 3073, 3181, 3497, 3499 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Sequence is finite, for proof see A175607.
Search for terms can be restricted to the range from 2 to A175607(16) = 2907159732049; primepi(53) = 16.
LINKS
MATHEMATICA
jj = 2^36*3^23*5^15*7^13*11^10*13^9*17^8*19^8*23^8*29^7*31^7*37^7*41^6 *43^6*47^6*53^6*59^6*61^6*67^6*71^5*73^5*79^5*83^5*89^5*97^5; rr = {}; n = 2; While[n < 3222617400, If[GCD[jj, n^2 - 1] == n^2 - 1, k = FactorInteger[n^2 - 1]; kk = Last[k][[1]]; If[kk == 53, AppendTo[rr, n]]]; n++ ]; rr (* Artur Jasinski *)
Select[Range[300000], FactorInteger[#^2-1][[-1, 1]]==53&]
PROG
(Magma) [ n: n in [2..300000] | m eq 53 where m is D[#D] where D is PrimeDivisors(n^2-1) ]; // Klaus Brockhaus, Feb 19 2011
(Magma) p:=(97*89*83*79*73*71)^5 *(67*61*59*53*47*43*41)^6 *(37*31*29)^7 *(23*19*17)^8 *13^9 *11^10 *7^13 *5^15 *3^23 *2^36; [ n: n in [2..50000000] | p mod (n^2-1) eq 0 and (D[#D] eq 53 where D is PrimeDivisors(n^2-1)) ]; // Klaus Brockhaus, Feb 20 2011
(PARI) is(n)=n=n^2-1; forprime(p=2, 47, n/=p^valuation(n, p)); n>1 && 53^valuation(n, 53)==n \\ Charles R Greathouse IV, Jul 01 2013
CROSSREFS
Sequence in context: A327109 A327108 A295156 * A346455 A050422 A350232
KEYWORD
fini,nonn
AUTHOR
Artur Jasinski, Oct 21 2010
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 09:04 EDT 2024. Contains 371240 sequences. (Running on oeis4.)