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.

A083722 Product of primes greater than the greatest prime factor of n but not greater than n.

This page as a plain text file.
%I A083722 #18 Nov 02 2017 15:35:50
%S A083722 1,1,1,3,1,5,1,105,35,7,1,385,1,143,1001,15015,1,85085,1,323323,46189,
%T A083722 4199,1,37182145,7436429,7429,37182145,1062347,1,215656441,1,
%U A083722 100280245065,86822723,392863,955049953,33426748355,1,765049,247110827,247357937827,1,1448810778701,1
%N A083722 Product of primes greater than the greatest prime factor of n but not greater than n.
%C A083722 a(n) = 1 iff n is prime or n is 1.
%C A083722 Apart from 1-terms, the other duplicated terms are: a(24) = a(27), a(120) = a(125), a(140) = a(147), a(528) = a(539), etc, whose positions are listed by A293893 and A293894. - _Antti Karttunen_, Nov 01 2017
%H A083722 Antti Karttunen, <a href="/A083722/b083722.txt">Table of n, a(n) for n = 1..2001</a>
%F A083722 a(n) = A002110(A049084(A007917(n)))/A002110(A049084(A006530(n))).
%t A083722 Array[Times @@ Select[Prime@ Range[#1, #1 + #2], Function[p, p <= #3]] & @@ {PrimePi@ NextPrime[FactorInteger[#][[-1, 1]]], PrimePi@ #, #} &, 43] (* _Michael De Vlieger_, Nov 01 2017 *)
%o A083722 (PARI) a(n) = {if (n==1, return (1)); my(gpf = vecmax(factor(n)[,1])); my(pp = 1); forprime(p=gpf+1, n, pp *= p;); pp;} \\ _Michel Marcus_, Jun 26 2016
%Y A083722 Cf. A006530, A083721, A083720.
%Y A083722 Cf. A293892 (restricted growth sequence transform), A293893, A293894.
%K A083722 nonn
%O A083722 1,4
%A A083722 _Reinhard Zumkeller_, May 04 2003
%E A083722 More terms from _Michel Marcus_, Jun 26 2016