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.

A081761 Least number >n having same type of prime factorization, a(1)=1.

This page as a plain text file.
%I A081761 #14 Dec 26 2020 03:49:27
%S A081761 1,3,5,9,7,10,11,27,25,14,13,18,17,15,21,81,19,20,23,28,22,26,29,40,
%T A081761 49,33,125,44,31,42,37,243,34,35,38,100,41,39,46,54,43,66,47,45,50,51,
%U A081761 53,80,121,52,55,63,59,56,57,88,58,62,61,84,67,65,68,729,69,70,71,75,74,78
%N A081761 Least number >n having same type of prime factorization, a(1)=1.
%H A081761 Amiram Eldar, <a href="/A081761/b081761.txt">Table of n, a(n) for n = 1..10000</a>
%F A081761 A046523(a(n)) = A046523(n).
%F A081761 a(n) = n + A085072(n). - _Amiram Eldar_, Dec 26 2020
%t A081761 sig[n_] := Sort @ FactorInteger[n][[;; , 2]]; a[1] = 1; a[n_] := Module[{sign = sig[n], k = n + 1}, While[sig[k] != sign, k++]; k]; Array[a, 70] (* _Amiram Eldar_, Dec 26 2020 *)
%o A081761 (PARI) a(n) = {if (n==1, 1, my(k=1, s = vecsort(factor(n)[,2]~)); while (vecsort(factor(n+k)[,2]~) != s, k++); n+k;)} \\ _Michel Marcus_, Nov 02 2020
%Y A081761 Cf. A046523, A085072.
%K A081761 nonn
%O A081761 1,2
%A A081761 _Reinhard Zumkeller_, Apr 09 2003