A383727 Main diagonal of A383726.
30, 6279, 8970, 2532915, 3327870
Offset: 3
Crossrefs
Cf. A383726.
Programs
Formula
a(n) = A383726(n,n-2).
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.
32085 is a term because it has 4 distinct prime factors (3, 5, 23 and 31) and the largest one is the sum of the others (3 + 5 + 23 = 31).
N:= 10^5: # for terms <= N P:= select(isprime,[2,seq(i,i=3..N/(3*5*7),2)]): V:= NULL: for j from 1 while P[j]^3*(3*P[j]) < N do for k from j+1 while P[j]*P[k]^2*(P[j]+2*P[k]) < N do for l from k+1 while P[j]*P[k]*P[l] * (P[j]+P[k]+P[l]) <= N do p4:= P[j]+P[k]+P[l]; if not isprime(p4) then next fi; for d1 from 1 while P[j]^d1 * P[k] * P[l] * p4 <= N do for d2 from 1 while P[j]^d1 * P[k]^d2 * P[l] * p4 <= N do for d3 from 1 while P[j]^d1 * P[k]^d2 * P[l]^d3 * p4 <= N do for d4 from 1 while P[j]^d1 * P[k]^d2 * P[l]^d3 * p4^d4 <= N do V:= V,P[j]^d1 * P[k]^d2 * P[l]^d3 * p4^d4 od od od od od od od: sort([V]); # Robert Israel, Jun 09 2025
A383728Q[k_] := Length[#] == 4 && Total[Most[#]] == Last[#] & [FactorInteger[k][[All, 1]]]; Select[Range[10^5], A383728Q]
10710 is a term because it has 5 distinct prime factors (2, 3, 5, 7 and 17) and the largest one is the sum of the others (2 + 3 + 5 + 7 = 17).
N:= 10^5: # for terms <= N P:= select(isprime,[2,seq(i,i=3..N/(2*3*5*7),2)]): V:= NULL: i:= 1: for j from i+1 while P[i]*P[j]^3*(P[i]+3*P[j]) < N do for k from j+1 while P[i]*P[j]*P[k]^2*(P[i]+P[j]+2*P[k]) < N do for l from k+1 while P[i]*P[j]*P[k]*P[l] * (P[i]+P[j]+P[k]+P[l]) <= N do p5:= P[i]+P[j]+P[k]+P[l]; if not isprime(p5) then next fi; for d1 from 1 while P[i]^d1 * P[j] * P[k] * P[l] * p5 <= N do for d2 from 1 while P[i]^d1 * P[j]^d2 * P[k] * P[l] * p5 <= N do for d3 from 1 while P[i]^d1 * P[j]^d2 * P[k]^d3 * P[l] * p5 <= N do for d4 from 1 while P[i]^d1 * P[j]^d2 * P[k]^d3 * P[l]^d4 * p5 <= N do for d5 from 1 while P[i]^d1 * P[j]^d2 * P[k]^d3 * P[l]^d4 * p5^d5 <= N do V:= V,P[i]^d1 * P[j]^d2 * P[k]^d3 * P[l]^d4 * p5^d5 od od od od od od od od: sort([V]); # Robert Israel, Jun 09 2025
A383729Q[k_] := Length[#] == 5 && Total[Most[#]] == Last[#] & [FactorInteger[k][[All, 1]]]; Select[Range[10^5], A383729Q]
a(3) = 30 is the smallest number having 3 distinct prime factors (namely 2, 3, and 5) such that the largest one is the sum of the others (2 + 3 = 5). a(4) = 3135 is the smallest number having 4 distinct prime factors (namely 3, 5, 11 and 19) such that the largest one is the sum of the others (3 + 5 + 11 = 19).
isok(k, n) = my(f=factor(k)); (omega(f)==n) && (vecsum(f[,1]) == 2*vecmax(f[,1])); a(n) = my(k=1); while (!isok(k, n), k++); k; \\ Michel Marcus, May 08 2025
T(4,3) = 1080 is the smallest 4-digit number having 3 distinct prime factors (namely 2, 3, and 5) such that the largest one is the sum of the others (2 + 3 = 5). T(5,4) = 10695 is the smallest 5-digit number having 4 distinct prime factors (namely 3, 5, 23 and 31) such that the largest one is the sum of the others (3 + 5 + 23 = 31). Triangle begins: 30; 120, -1; 1080, 3135, 3570; 10336, 10695, 10626, -1; 100672, 102695, 103730, 844305, -1; 1003520, 1005039, 1003450, 1218945, 1231230, -1; ...
T(4,4) = 2145 = 3*5*11*13 is the smallest four-digit number with four distinct prime factors (3, 5, 11, and 13), where the prime factors can be partitioned into two subsets of at least two elements each, both summing to the same value: 3+13 = 5+11. T(5,4) = 2310 = 2*3*5*7*11 is the smallest five-digit number with four distinct prime factors (2, 3, 5, 7 and 11), where the prime factors can be partitioned into two subsets of at least two elements each, both summing to the same value: 3+11 = 2+5+7. The lower triangle begins at T(4,4): [ 2145, 2310]; [ 10725, 10374]; [ 101065, 100050, 255255, 510510]; [ 1005993, 1000350, 1036035, 1009470]; [ 10006081, 10000130, 10012065, 10004610]; [100010225, 100001300, 100001195, 100009910, 111546435, 223092870]; ...
Comments