OFFSET
1,13
COMMENTS
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
The prime indices of 2257 are {12,18}, with common prime factors {2,3}, so a(2257) = 2.
MAPLE
f:= proc(n) local g;
g:= igcd(op(map(NumberTheory:-pi, NumberTheory:-PrimeFactors(n))));
nops(NumberTheory:-PrimeFactors(g))
end proc: f(1):= 0:
map(f, [$1..100]); # Robert Israel, Jan 13 2026
MATHEMATICA
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Table[If[n==1, 0, PrimeNu[GCD@@prix[n]]], {n, 100}]
CROSSREFS
For initial interval instead of prime factors we have A055396.
For LCM (or product) instead of GCD we have A303975.
For all divisors instead of just prime factors we have A387114.
A000005 counts divisors.
A003963 multiplies together prime indices.
A120383 lists numbers divisible by all of their prime indices.
A289508 gives greatest common divisor of prime indices.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Sep 11 2025
STATUS
approved
