|
%I #39 Jan 16 2023 08:07:31
%S 3,4,6,8,10,18,22,24,30,32,40,42,46,54,56,66,70,72,88,96,98,102,114,
%T 118,126,128,130,136,150,152,168,182,200,224,234,238,246,250,266,270,
%U 294,312,318,328,330,350,354,360,370,392,402,406,416,424,434,440,442,450,472,480,486,510,536,546,594,600,630,640,646,648,650,654,666,680,690,722
%N Numbers k such that A048675(k) is a prime.
%C After 3 and 4 each term is an even number with an odd exponent of 2. - _David A. Corneth_ and _Antti Karttunen_, Oct 11 2016
%H Antti Karttunen (terms 1..4994) & Hans Havermann, <a href="/A277319/b277319.txt">Table of n, a(n) for n = 1..25000</a>
%H Hans Havermann, <a href="http://chesswanks.com/num/(a277319,a277321).txt">70000 terms with their associated primes</a>
%o (PARI)
%o allocatemem(2^30);
%o A048675(n) = my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; \\ From _Michel Marcus_, Oct 10 2016
%o isA277319 = n -> isprime(A048675(n));
%o i=0; n=1; while(i < 10000, n++; if(isA277319(n), i++; write("b277319.txt", i, " ", n)));
%o (Python)
%o from sympy import factorint, primepi, isprime
%o def a048675(n):
%o if n==1: return 0
%o f=factorint(n)
%o return sum([f[i]*2**(primepi(i) - 1) for i in f])
%o print([n for n in range(1, 1001) if isprime(a048675(n))]) # _Indranil Ghosh_, Jun 19 2017
%Y Row 1 of A277898. Positions of ones in A277892.
%Y Cf. A048675 and A277321 for the primes themselves.
%Y Cf. A277317 (a subsequence).
%Y After two initial terms a subsequence of A036554.
%K nonn
%O 1,1
%A _Antti Karttunen_, Oct 11 2016
|
