A080784 Duplicate of A072978.
1, 6, 10, 14, 22, 26, 34, 36, 38, 46, 58, 60, 62, 74, 82, 84, 86, 94, 100, 106, 118, 122
Offset: 1
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.
a028260 n = a028260_list !! (n-1) a028260_list = filter (even . a001222) [1..] -- Reinhard Zumkeller, Oct 05 2011
with(numtheory); A028260 := proc(n) option remember: local k: if(n=1)then return 1: fi: for k from procname(n-1)+1 do if(bigomega(k) mod 2=0)then return k: fi: od: end: seq(A028260(n),n=1..63); # Nathaniel Johnston, May 27 2011
Select[Range[200],EvenQ[PrimeOmega[#]]&] (* Harvey P. Dale, Aug 14 2011 *) Select[Range@ 134, LiouvilleLambda@# > 0 &] (* Robert G. Wilson v, Jul 06 2012 *)
is(n)=bigomega(n)%2==0 \\ Charles R Greathouse IV, May 29 2013
from math import isqrt, prod from sympy import primerange, primepi, integer_nthroot def A028260(n): def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1))) def f(x): return int(n+x-1-sum(sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,m)) for m in range(2,x.bit_length()+1,2))) m, k = n, f(n) while m != k: m, k = k, f(k) return m # Chai Wah Wu, Apr 10 2025
The terms together with their prime indices begin: 1: {} 6: {1,2} 14: {1,4} 15: {2,3} 26: {1,6} 33: {2,5} 35: {3,4} 36: {1,1,2,2} 38: {1,8} 51: {2,7} 58: {1,10} 60: {1,1,2,3}
prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]; Select[Range[100],!IntegerQ[Median[prix[#]]]&]
The 0-prepended prime indices of 1617 are {0,2,4,4,5}, with sorted differences {0,1,2,2}, with median 3/2, so 1617 is not in the sequence.
prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]; Select[Range[2,100],IntegerQ[Median[Differences[Prepend[prix[#],0]]]]&]
The prime indices of 900 are {1,1,2,2,3,3}, with distinct parts {1,2,3}, with median 2, so 900 is in the sequence. The prime indices of 330 are {1,2,3,5}, with distinct parts {1,2,3,5}, with median 5/2, so 330 is not in the sequence.
Select[Range[2,100],IntegerQ[Median[PrimePi/@First/@FactorInteger[#]]]&]
The definition specifies that we count repeated prime factors. 6 = 2 * 3; the sum of these prime factors is 2 + 3 = 5, an odd number; so 6 is not in the sequence. 50 = 2 * 5 * 5 has 3 prime factors and 3 is an odd number; so 50 is not in the sequence. 60 = 2 * 2 * 3 * 5 has 4 prime factors and 4 is an even number; the sum of these factors is 2 + 2 + 3 + 5 = 12, also an even number; so 60 is in the sequence. 1 has 0 prime factors, which sum to 0 (the empty sum). 0 is even, so 1 is in the sequence.
{1}~Join~Select[Range@1000,(s=Flatten[Table@@@FactorInteger[#]];And@@EvenQ@{Length@s,Total@s})&] (* Giorgos Kalogeropoulos, Jun 24 2021 *)
iseven(x) = ((x%2) == 0); isok(m) = my(f=factor(m)); iseven(sum(k=1, #f~, f[k,1]*f[k,2])) && iseven(sum(k=1, #f~, f[k,2])); \\ Michel Marcus, Jun 24 2021
is(n) = bigomega(n)%2 == 0 && valuation(n, 2)%2 == 0 \\ David A. Corneth, Jun 24 2021
from sympy import factorint def ok(n): f = factorint(n) return sum(f.values())%2 == 0 and sum(p*f[p] for p in f)%2 == 0 print(list(filter(ok, range(1, 216)))) # Michael S. Branicky, Jun 24 2021
The prime factors of 60 are {2,2,3,5}, with median 5/2, so a(60) = 5.
Table[2*Median[Join@@ConstantArray@@@FactorInteger[n]],{n,100}]
The prime indices of 900 are {1,1,2,2,3,3}, with distinct parts {1,2,3}, with median 2, so 900 is not in the sequence. The prime indices of 462 are {1,2,4,5}, with distinct parts {1,2,4,5}, with median 3, so 462 is not in the sequence.
Select[Range[2,100],!IntegerQ[Median[PrimePi/@First/@FactorInteger[#]]]&]
The unordered prime signature of 2520 is {3,2,1,1}, with median 3/2, so 2520 is in the sequence. The unordered prime signature of 12600 is {3,2,2,1}, with median 2, so 12600 is not in the sequence.
Select[Range[2,100],!IntegerQ[Median[Last/@FactorInteger[#]]]&]
The 0-prepended prime indices of 1617 are {0,2,4,4,5}, with sorted differences {0,1,2,2}, with median 3/2, so 1617 is in the sequence.
prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]; Select[Range[2,100],!IntegerQ[Median[Differences[Prepend[prix[#],0]]]]&]
Comments