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.

A076426 Fixed points of the mapping k -> abs(reverse(lpd(k))-reverse(Lpf(k))). lpd(k) is the largest proper divisor and Lpf(k) is the largest prime factor of k.

This page as a plain text file.
%I A076426 #6 Mar 30 2012 17:34:50
%S A076426 5750,33866,74841,517250,577750,5538710,51414250,51454250,51687250,
%T A076426 51727250,51748250,51858250,52525250,57515750,57535750,57575750,
%U A076426 57757750,67597352,841794296,5120202250,5120802250,5121612250
%N A076426 Fixed points of the mapping k -> abs(reverse(lpd(k))-reverse(Lpf(k))). lpd(k) is the largest proper divisor and Lpf(k) is the largest prime factor of k.
%C A076426 Besides these fixed points (cycles of length 1) there are five cycles of length 2 ([9378, 9739], [518775, 522075], [5170250, 5197250], [5219475, 5249775], [5255750, 5755250]) and one cycle of length 3 ([7285, 7467, 9711]) below 8000000.
%F A076426 abs(reverse(lpd(n))-reverse(Lpf(n))) = n.
%e A076426 lpd(5750) = 2875; Lpf(5750) = 23; 5782 - 32 = 5750.
%o A076426 (PARI) {for(n=1,34000,v=divisors(n); a=matsize(v)[2]; z=if(a>1,v[a-1],1); p=0; while(z>0,d=divrem(z,10); z=d[1]; p=10*p+d[2]); z=if(n==1,1,vecmax(component(factor(n),1))); q=0; while(z>0,d=divrem(z,10); z=d[1]; q=10*q+d[2]); if(abs(p-q)==n,print1(n,",")))}
%Y A076426 Cf. A076423, A076425.
%K A076426 base,nonn
%O A076426 1,1
%A A076426 _Klaus Brockhaus_, Oct 11 2002
%E A076426 Offset corrected and a(7)-a(22) from _Donovan Johnson_, Aug 09 2010