A001026 -id:A001026 - cp's OEIS frontend

A091045 Partial sums of powers of 17 (A001026).

1, 18, 307, 5220, 88741, 1508598, 25646167, 435984840, 7411742281, 125999618778, 2141993519227, 36413889826860, 619036127056621, 10523614159962558, 178901440719363487, 3041324492229179280, 51702516367896047761

Offset: 1

Wolfdieter Lang, Jan 23 2004

17^a(n) is largest power of 17 dividing (17^n)!.

Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=17, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n)=det(A). - Milan Janjic, Feb 21 2010

Robert Israel, Table of n, a(n) for n = 1..812
Index entries for linear recurrences with constant coefficients, signature (18, -17).

Cf. similar sequences of the form (k^n-1)/(k-1) with k prime: A000225 (k=2), A003462 (k=3), A003463 (k=5), A023000 (k=7), A016123 (k=11), A091030 (k=13), this sequence (k=17), A218722 (k=19), A218726 (k=23), A218732 (k=29), A218734 (k=31), A218740 (k=37), A218744 (k=41), A218746 (k=43), A218750 (k=47).

Cf. A001026.

[&+[17^i: i in [0..n]]: n in [0..20]]; // Vincenzo Librandi, Feb 19 2018

ListTools:-PartialSums([seq(17^k,k=0..30)]); # Robert Israel, Feb 18 2018

Table[17^n, {n, 0, 16}] // Accumulate (* Jean-François Alcover, Jul 05 2013 *)

makelist(sum(17^k, k, 0, n), n, 0, 30); /* Martin Ettl, Nov 05 2012 */

[gaussian_binomial(n,1,17) for n in range(1,18)] # Zerinvary Lajos, May 28 2009

a(n) = Sum_{k=0..n-1} 17^k = (17^n - 1)/16.
G.f.: x/((1 - 17*x)*(1 - x))= (1/(1 - 17*x) - 1/(1 - x))/16.
a(n) = 17*a(n-1)+1 (with a(1)=1). - Vincenzo Librandi, Nov 16 2010
E.g.f.: exp(9*x)*sinh(8*x)/8. - Stefano Spezia, Mar 11 2023

A060220 Number of orbits of length n under the full 17-shift (whose periodic points are counted by A001026).

Original entry on oeis.org

17, 136, 1632, 20808, 283968, 4022064, 58619808, 871959240, 13176430176, 201599248032, 3115626937056, 48551851084080, 761890617915840, 12026987582075856, 190828203433892736, 3041324491793194440, 48661191875666868480, 781282469552728498992, 12582759772902701307744

Offset: 1

Thomas Ward, Mar 21 2001

Number of monic irreducible polynomials of degree n over GF(17). - Andrew Howroyd, Dec 10 2017

			a(2)=136 since there are 289 points of period 2 in the full 17-shift and 17 fixed points, so there must be (289-17)/2 = 136 orbits of length 2.

G. C. Greubel, Table of n, a(n) for n = 1..810
Y. Puri and T. Ward, Arithmetic and growth of periodic orbits, J. Integer Seqs., Vol. 4 (2001), #01.2.1.
Yash Puri and Thomas Ward, A dynamical property unique to the Lucas sequence, Fibonacci Quarterly, Volume 39, Number 5 (November 2001), pp. 398-402.
T. Ward, Exactly realizable sequences

Column 17 of A074650.

Cf. A001026, A008683.

A060220:= func< n | (1/n)*(&+[MoebiusMu(d)*(17)^Floor(n/d): d in Divisors(n)]) >;
[A060220(n): n in [1..40]]; // G. C. Greubel, Sep 13 2024

A060220[n_]:= DivisorSum[n, (17)^(n/#)*MoebiusMu[#] &]/n;
Table[A060220[n], {n,40}] (* G. C. Greubel, Sep 13 2024 *)

a001024(n) = 17^n;
a(n) = (1/n)*sumdiv(n, d, moebius(d)*a001024(n/d)); \\ Michel Marcus, Sep 11 2017

def A060220(n): return (1/n)*sum(moebius(k)*(17)^(n/k) for k in (1..n) if (k).divides(n))
[A060220(n) for n in range(1,41)] # G. C. Greubel, Sep 13 2024

a(n) = (1/n)* Sum_{d|n} mu(d)*A001026(n/d).

G.f.: Sum_{k>=1} mu(k)*log(1/(1 - 17*x^k))/k. - Ilya Gutkovskiy, May 20 2019

More terms from Michel Marcus, Sep 11 2017

A008472 Sum of the distinct primes dividing n.

Original entry on oeis.org

0, 2, 3, 2, 5, 5, 7, 2, 3, 7, 11, 5, 13, 9, 8, 2, 17, 5, 19, 7, 10, 13, 23, 5, 5, 15, 3, 9, 29, 10, 31, 2, 14, 19, 12, 5, 37, 21, 16, 7, 41, 12, 43, 13, 8, 25, 47, 5, 7, 7, 20, 15, 53, 5, 16, 9, 22, 31, 59, 10, 61, 33, 10, 2, 18, 16, 67, 19, 26, 14, 71, 5, 73

Offset: 1

Olivier Gérard

Sometimes called sopf(n).
Sum of primes dividing n (without repetition) (compare A001414).
Equals A051731 * A061397 = inverse Mobius transform of [0, 2, 3, 0, 5, 0, 7, ...]. - Gary W. Adamson, Feb 14 2008
Equals row sums of triangle A143535. - Gary W. Adamson, Aug 23 2008
a(n) = n if and only if n is prime. - Daniel Forgues, Mar 24 2009
a(n) = n is a new record if and only if n is prime. - Zak Seidov, Jun 27 2009
a(A001043(n)) = A191583(n);
For n > 0: a(A000079(n)) = 2, a(A000244(n)) = 3, a(A000351(n)) = 5, a(A000420(n)) = 7;
a(A006899(n)) <= 3; a(A003586(n)) = 5; a(A033846(n)) = 7; a(A033849(n)) = 8; a(A033847(n)) = 9; a(A033850(n)) = 10; a(A143207(n)) = 10. - Reinhard Zumkeller, Jun 28 2011
For n > 1: a(n) = Sum(A027748(n,k): 1 <= k <= A001221(n)). - Reinhard Zumkeller, Aug 27 2011
If n is the product of twin primes (A037074), a(n) = 2*sqrt(n+1) = sqrt(4n+4). - Wesley Ivan Hurt, Sep 07 2013
From Wilf A. Wilson, Jul 21 2017: (Start)
a(n) + 2, n > 2, is the number of maximal subsemigroups of the monoid of orientation-preserving or -reversing mappings on a set with n elements.
a(n) + 3, n > 2, is the number of maximal subsemigroups of the monoid of orientation-preserving or -reversing partial mappings on a set with n elements.
(End)
The smallest m such that a(m) = n, or 0 if no such number m exists is A064502(n). The only integers that are not in the sequence are 1, 4 and 6. - Bernard Schott, Feb 07 2022

			a(18) = 5 because 18 = 2 * 3^2 and 2 + 3 = 5.
a(19) = 19 because 19 is prime.
a(20) = 7 because 20 = 2^2 * 5 and 2 + 5 = 7.

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 10000 terms from Franklin T. Adams-Watters)
Johann Bartel, R. K. Bhaduri, Matthias Brack, and M. V. N. Murthy, On the asymptotic prime partitions of integers, arXiv:1609.06497 [math-ph], 2017.
James East, Jitender Kumar, James D. Mitchell, and Wilf A. Wilson, Maximal subsemigroups of finite transformation and partition monoids, arXiv:1706.04967 [math.GR], 2017. [_Wilf A. Wilson_, Jul 21 2017]

First difference of A024924.
Cf. A001414 (sopfr), A001222, A051731, A061397, A064502, A085020, A143535.
Sum of the k-th powers of the primes dividing n for k=0..10 : A001221 (k=0), this sequence (k=1), A005063 (k=2), A005064 (k=3), A005065 (k=4), A351193 (k=5), A351194 (k=6), A351195 (k=7), this sequence (k=8), A351197 (k=9), A351198 (k=10).
Cf. A010051.

a008472 = sum . a027748_row  -- Reinhard Zumkeller, Mar 29 2012

[n eq 1 select 0 else &+[p[1]: p in Factorization(n)]: n in [1..100]]; // Vincenzo Librandi, Jun 24 2017

A008472 := n -> add(d, d = select(isprime, numtheory[divisors](n))):
seq(A008472(i), i = 1..40); # Peter Luschny, Jan 31 2012
A008472 := proc(n)
        add( d, d= numtheory[factorset](n)) ;
end proc: # R. J. Mathar, Jul 08 2012

Prepend[Array[Plus @@ First[Transpose[FactorInteger[#]]] &, 100, 2], 0]
Join[{0}, Rest[Total[Transpose[FactorInteger[#]][[1]]]&/@Range[100]]] (* Harvey P. Dale, Jun 18 2012 *)
(* Requires version 7.0+ *) Table[DivisorSum[n, # &, PrimeQ[#] &], {n, 75}] (* Alonso del Arte, Dec 13 2014 *)
Table[Sum[p, {p, Select[Divisors[n], PrimeQ]}], {n, 1, 100}] (* Vaclav Kotesovec, May 20 2020 *)

sopf(n) = local(fac=factor(n)); sum(i=1,matsize(fac)[1],fac[i,1])

vector(100,n,vecsum(factor(n)[,1]~)) \\ Derek Orr, May 13 2015

A008472(n)=vecsum(factor(n)[,1]) \\ M. F. Hasler, Jul 18 2015

from sympy import primefactors
def A008472(n): return sum(primefactors(n)) # Chai Wah Wu, Feb 03 2022

def A008472(n):
    return add(d for d in divisors(n) if is_prime(d))
print([A008472(i) for i in (1..40)]) # Peter Luschny, Jan 31 2012

[sum(prime_factors(n)) for n in range(1,74)] # Giuseppe Coppoletta, Jan 19 2015

Let n = Product_j prime(j)^k(j) where k(j) >= 1, then a(n) = Sum_j prime(j).
Additive with a(p^e) = p.
G.f.: Sum_{k >= 1} prime(k)*x^prime(k)/(1-x^prime(k)). - Franklin T. Adams-Watters, Sep 01 2009
L.g.f.: -log(Product_{k>=1} (1 - x^prime(k))) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, May 06 2017
Dirichlet g.f.: primezeta(s-1)*zeta(s). - Benedict W. J. Irwin, Jul 11 2018
a(n) = Sum_{p|n, p prime} p. - Wesley Ivan Hurt, Feb 04 2022
From Bernard Schott, Feb 07 2022: (Start)
For n > 0: a(A001020(n)) = 11, a(A001022(n)) = 13, a(A001026(n)) = 17, a(A001029(n)) = 19, a(A009967(n)) = 23, a(A009973(n)) = 29, a(A009975(n)) = 31, a(A009981(n)) = 37, a(A009985(n)) = 41, a(A009987(n)) = 43, a(A009991(n)) = 47.
For p odd prime, a(2*p) = p+2 <==> a(A100484(n)) = A052147(n) for n > 1. (End)
a(n) = Sum_{d|n} d * c(d), where c = A010051. - Wesley Ivan Hurt, Jun 22 2024

A329697 a(n) is the number of iterations needed to reach a power of 2 starting at n and using the map k -> k-(k/p), where p is the largest prime factor of k.

Original entry on oeis.org

0, 0, 1, 0, 1, 1, 2, 0, 2, 1, 2, 1, 2, 2, 2, 0, 1, 2, 3, 1, 3, 2, 3, 1, 2, 2, 3, 2, 3, 2, 3, 0, 3, 1, 3, 2, 3, 3, 3, 1, 2, 3, 4, 2, 3, 3, 4, 1, 4, 2, 2, 2, 3, 3, 3, 2, 4, 3, 4, 2, 3, 3, 4, 0, 3, 3, 4, 1, 4, 3, 4, 2, 3, 3, 3, 3, 4, 3, 4, 1, 4, 2, 3, 3, 2, 4, 4, 2, 3, 3, 4, 3, 4, 4, 4, 1, 2, 4, 4, 2

Offset: 1

Ali Sada and Robert G. Wilson v, Feb 28 2020

From Antti Karttunen, Apr 07 2020: (Start)
Also the least number of iterations of nondeterministic map k -> k-(k/p) needed to reach a power of 2, when any prime factor p of k can be used. The minimal length path to the nearest power of 2 (= 2^A064415(n)) is realized whenever one uses any of the A005087(k) distinct odd prime factors of the current k, at any step of the process. For example, this could be done by iterating with the map k -> k-(k/A078701(k)), i.e., by using the least odd prime factor of k (instead of the largest prime).
Proof: Viewing the prime factorization of changing k as a multiset ("bag") of primes, we see that liquefying any odd prime p with step p -> (p-1) brings at least one more 2 to the bag, while applying p -> (p-1) to any 2 just removes it from the bag, but gives nothing back. Thus the largest (and thus also the nearest) power of 2 is reached by eliminating - step by step - all odd primes from the bag, but none of 2's, and it doesn't matter in which order this is done.
The above implies also that the sequence is totally additive, which also follows because both A064097 and A064415 are. That A064097(n) = A329697(n) + A054725(n) for all n > 1 can be also seen by comparing the initial conditions and the recursion formulas of these three sequences.
For any n, A333787(n) is either the nearest power of 2 reached (= 2^A064415(n)), or occurs on some of the paths from n to there.
(End)
A003401 gives the numbers k where a(k) = A005087(k). See also A336477. - Antti Karttunen, Mar 16 2021

			The trajectory of 15 is {12, 8}, taking 2 iterations to reach 8 = 2^3. So a(15) is 2.
From _Antti Karttunen_, Apr 07 2020: (Start)
Considering all possible paths from 15 to 1 nondeterministic map k -> k-(k/p), where p can be any prime factor of k, we obtain the following graph:
        15
       / \
      /   \
    10     12
    / \   / \
   /   \ /   \
  5     8     6
   \__  |  __/|
      \_|_/   |
        4     3
         \   /
          \ /
           2
           |
           1.
It can be seen that there's also alternative route to 8 via 10 (with 10 = 15-(15/3), where 3 is not the largest prime factor of 15), but it's not any shorter than the route via 12.
(End)

Antti Karttunen, Table of n, a(n) for n = 1..65537
Michael De Vlieger, Annotated fan style binary tree labeling the index n, with a color code where black represents a(n) = 0, red a(n) = 1, and magenta the largest value in a(n) for n = 1..16383.

Cf. A000265, A003401, A005087, A052126, A053575, A054725, A064097, A064415, A078701, A087436, A147545, A171462, A209229, A333123, A333787, A333790, A334107, A334109, A335875, A334204, A335880, A335881, A336396, A336466, A336469 [= a(phi(n))], A336928 [= a(sigma(n))], A336470, A336477, A339970.
Cf. A000079, A334101, A334102, A334103, A334104, A334105, A334106 for positions of 0 .. 6 in this sequence, and also array A334100.
Cf. A334099 (a right inverse, positions of the first occurrence of each n).
Cf. A334091 (first differences), A335429 (partial sums).
Cf. also A331410 (analogous sequence when using the map k -> k + k/p), A334861, A335877 (their sums and differences), see also A335878 and A335884, A335885.

a[n_] := Length@ NestWhileList[# - #/FactorInteger[#][[-1, 1]] &, n, # != 2^IntegerExponent[#, 2] &] -1; Array[a, 100]

A329697(n) = if(!bitand(n,n-1),0,1+A329697(n-(n/vecmax(factor(n)[, 1])))); \\ Antti Karttunen, Apr 07 2020

up_to = 2^24;
A329697list(up_to) = { my(v=vector(up_to)); v[1] = 0; for(n=2, up_to, v[n] = if(!bitand(n,n-1),0,1+vecmin(apply(p -> v[n-n/p], factor(n)[, 1]~)))); (v); };
v329697 = A329697list(up_to);
A329697(n) = v329697[n]; \\ Antti Karttunen, Apr 07 2020

A329697(n) = if(n<=2,0, if(isprime(n), A329697(n-1)+1, my(f=factor(n)); (apply(A329697, f[, 1])~ * f[, 2]))); \\ Antti Karttunen, Apr 19 2020

From Antti Karttunen, Apr 07-19 2020: (Start)
a(1) = a(2) = 0; and for n > 2, a(p) = 1 + a(p-1) if p is an odd prime and a(n*m) = a(n) + a(m) if m,n > 1. [This is otherwise equal to the definition of A064097, except here we have a different initial condition, with a(2) = 0].
a(2n) = a(A000265(n)) = a(n).
a(p) = 1+a(p-1), for all odd primes p.
If A209229(n) == 1 [when n is a power of 2], a(n) = 0,
  otherwise a(n) = 1 + a(n-A052126(n)) = 1 + a(A171462(n)).
Equivalently, for non-powers of 2, a(n) = 1 + a(n-(n/A078701(n))),
or equivalently, for non-powers of 2, a(n) = 1 + Min a(n - n/p), for p prime and dividing n.
a(n) = A064097(n) - A064415(n), or equally, a(n) = A064097(n) - A054725(n), for n > 1.
a(A019434(n)) = 1, a(A334092(n)) = 2, a(A334093(n)) = 3, etc. for all applicable n.
For all n >= 0, a(A334099(n)) = a(A000244(n)) = a(A000351(n)) = a(A001026(n)) = a(257^n) = a(65537^n) = n.
a(A122111(n)) = A334107(n), a(A225546(n)) = A334109(n).
(End)
From Antti Karttunen, Mar 16 2021: (Start)
a(n) = a(A336466(n)) + A087436(n) = A336396(n) + A087436(n).
a(A053575(n)) = A336469(n) = a(n) - A005087(n).
a(A147545(n)) = A000120(A147545(n)) - 1.
(End)

A003992 Square array read by upwards antidiagonals: T(n,k) = n^k for n >= 0, k >= 0.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 3, 4, 1, 0, 1, 4, 9, 8, 1, 0, 1, 5, 16, 27, 16, 1, 0, 1, 6, 25, 64, 81, 32, 1, 0, 1, 7, 36, 125, 256, 243, 64, 1, 0, 1, 8, 49, 216, 625, 1024, 729, 128, 1, 0, 1, 9, 64, 343, 1296, 3125, 4096, 2187, 256, 1, 0, 1, 10, 81, 512, 2401, 7776, 15625, 16384, 6561, 512, 1, 0

Offset: 0

Marc LeBrun

If the array is transposed, T(n,k) is the number of oriented rows of n colors using up to k different colors. The formula would be T(n,k) = [n==0] + [n>0]*k^n. The generating function for column k would be 1/(1-k*x). For T(3,2)=8, the rows are AAA, AAB, ABA, ABB, BAA, BAB, BBA, and BBB. - Robert A. Russell, Nov 08 2018

T(n,k) is the number of multichains of length n from {} to [k] in the Boolean lattice B_k. - Geoffrey Critzer, Apr 03 2020

			Rows begin:
[1, 0,  0,   0,    0,     0,      0,      0, ...],
[1, 1,  1,   1,    1,     1,      1,      1, ...],
[1, 2,  4,   8,   16,    32,     64,    128, ...],
[1, 3,  9,  27,   81,   243,    729,   2187, ...],
[1, 4, 16,  64,  256,  1024,   4096,  16384, ...],
[1, 5, 25, 125,  625,  3125,  15625,  78125, ...],
[1, 6, 36, 216, 1296,  7776,  46656, 279936, ...],
[1, 7, 49, 343, 2401, 16807, 117649, 823543, ...], ...

T. D. Noe, Rows n = 0..50 of triangle, flattened

Rows 0-49 are A000007, A000012, A000079, A000244, A000302, A000351, A000400, A000420, A001018, A001019, A011557, A001020, A001021, A001022, A001023, A001024, A001025, A001026, A001027, A001029, A009964-A009992, A087752.
Columns 0-26 are A000012, A001477, A000290, A000578, A000583, A000584, A001014, A001015, A001016, A001017, A008454, A008455, A008456, A010801-A010813, A089081.
Main diagonal is A000312. Other diagonals include A000169, A007778, A000272, A008788. Antidiagonal sums are in A026898.
Cf. A099555.
Transpose is A004248. See A051128, A095884, A009999 for other versions.
Cf. A277504 (unoriented), A293500 (chiral).

[[(n-k)^k: k in [0..n]]: n in [0..10]]; // G. C. Greubel, Nov 08 2018

Table[If[k == 0, 1, (n - k)^k], {n, 0, 11}, {k, 0, n}]//Flatten

T(n,k) = (n-k)^k \\ Charles R Greathouse IV, Feb 07 2017

E.g.f.: Sum T(n,k)*x^n*y^k/k! = 1/(1-x*exp(y)). - Paul D. Hanna, Oct 22 2004

E.g.f.: Sum T(n,k)*x^n/n!*y^k/k! = e^(x*e^y). - Franklin T. Adams-Watters, Jun 23 2006

More terms from David W. Wilson

Edited by Paul D. Hanna, Oct 22 2004

A329332 Table of powers of squarefree numbers, powers of A019565(n) in increasing order in row n. Square array A(n,k) n >= 0, k >= 0 read by descending antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 8, 9, 6, 1, 1, 16, 27, 36, 5, 1, 1, 32, 81, 216, 25, 10, 1, 1, 64, 243, 1296, 125, 100, 15, 1, 1, 128, 729, 7776, 625, 1000, 225, 30, 1, 1, 256, 2187, 46656, 3125, 10000, 3375, 900, 7, 1, 1, 512, 6561, 279936, 15625, 100000, 50625, 27000, 49, 14

Offset: 0

Peter Munn, Nov 10 2019

The A019565 row order gives the table neat relationships with A003961, A003987, A059897, A225546, A319075 and A329050. See the formula section.

Transposition of this table, that is reflection about its main diagonal, has subtle symmetries. For example, consider the unique factorization of a number into powers of distinct primes. This can be restated as factorization into numbers from rows 2^n (n >= 0) with no more than one from each row. Reflecting about the main diagonal, this factorization becomes factorization (of a related number) into numbers from columns 2^k (k >= 0) with no more than one from each column. This is also unique and is factorization into powers of squarefree numbers with distinct exponents that are powers of two. See the example section.

			Square array A(n,k) begins:
n\k |  0   1     2      3        4          5           6             7
----+------------------------------------------------------------------
   0|  1   1     1      1        1          1           1             1
   1|  1   2     4      8       16         32          64           128
   2|  1   3     9     27       81        243         729          2187
   3|  1   6    36    216     1296       7776       46656        279936
   4|  1   5    25    125      625       3125       15625         78125
   5|  1  10   100   1000    10000     100000     1000000      10000000
   6|  1  15   225   3375    50625     759375    11390625     170859375
   7|  1  30   900  27000   810000   24300000   729000000   21870000000
   8|  1   7    49    343     2401      16807      117649        823543
   9|  1  14   196   2744    38416     537824     7529536     105413504
  10|  1  21   441   9261   194481    4084101    85766121    1801088541
  11|  1  42  1764  74088  3111696  130691232  5489031744  230539333248
  12|  1  35  1225  42875  1500625   52521875  1838265625   64339296875
Reflection of factorization about the main diagonal: (Start)
The canonical (prime power) factorization of 864 is 2^5 * 3^3 = 32 * 27. Reflecting the factors about the main diagonal of the table gives us 10 * 36 = 10^1 * 6^2 = 360. This is the unique factorization of 360 into powers of squarefree numbers with distinct exponents that are powers of two.
Reflection about the main diagonal is given by the self-inverse function A225546(.). Clearly, all positive integers are in the domain of A225546, whether or not they appear in the table. It is valid to start from 360, observe that A225546(360) = 864, then use 864 to derive 360's factorization into appropriate powers of squarefree numbers as above.
(End)

The range of values is A072774.
Rows (abbreviated list): A000079(1), A000244(2), A000400(3), A000351(4), A011557(5), A001024(6), A009974(7), A000420(8), A001023(9), A009965(10), A001020(16), A001022(32), A001026(64).
A019565 is column 1, A334110 is column 2, and columns that are sorted in increasing order (some without the 1) are: A005117(1), A062503(2), A062838(3), A113849(4), A113850(5), A113851(6), A113852(7).
Other subtables: A182944, A319075, A329050.
Re-ordered subtable of A297845, A306697, A329329.
A000290, A003961, A003987, A059897 and A225546 are used to express relationships between terms of this sequence.
Cf. A285322.

A(n,k) = A019565(n)^k.
A(k,n) = A225546(A(n,k)).
A(n,2k) = A000290(A(n,k)) = A(n,k)^2.
A(2n,k) = A003961(A(n,k)).
A(n,2k+1) = A(n,2k) * A(n,1).
A(2n+1,k) = A(2n,k) * A(1,k).
A(A003987(n,m), k) = A059897(A(n,k), A(m,k)).
A(n, A003987(m,k)) = A059897(A(n,m), A(n,k)).
A(2^n,k) = A319075(k,n+1).
A(2^n, 2^k) = A329050(n,k).
A(n,k) = A297845(A(n,1), A(1,k)) = A306697(A(n,1), A(1,k)), = A329329(A(n,1), A(1,k)).
Sum_{n>=0} 1/A(n,k) = zeta(k)/zeta(2*k), for k >= 2. - Amiram Eldar, Dec 03 2022

A128358 Numbers k such that k divides 18^k - 1.

Original entry on oeis.org

1, 17, 289, 4913, 83521, 1419857, 5846759, 24054337, 24137569, 99394903, 408923729, 410338673, 1689713351, 6951703393, 6975757441, 28725126967, 118178957681, 118285781329, 118587876497, 486643291847, 488327158439

Offset: 1

Alexander Adamchuk, Mar 01 2007

17 divides a(n) for n > 1.
a(n) = 17^(n-1) for n < 7.
a(7) = 20231*17^2.

Cf. A001026 (powers of 17), A015969 (k divides 16^k + 1).

is(n)=Mod(18,n)^n==1 \\ Charles R Greathouse IV, Nov 04 2016

a(8)-a(10) from Farideh Firoozbakht, Mar 05 2007
a(11)-a(19) from Ryan Propper, Dec 31 2007
a(20)-a(24) from Alexander Adamchuk, May 17 2010

A153652 Triangle T(n, k, j) = T(n-1, k, j) + T(n-1, k-1, j) + (2j + 1)prime(j)*T(n-2, k-1, j) with T(2, k, j) = prime(j) and j = 7, read by rows.

Original entry on oeis.org

2, 17, 17, 2, 574, 2, 2, 4911, 4911, 2, 2, 5423, 156192, 5423, 2, 2, 5935, 1413920, 1413920, 5935, 2, 2, 6447, 2802720, 42656800, 2802720, 6447, 2, 2, 6959, 4322592, 406009120, 406009120, 4322592, 6959, 2, 2, 7471, 5973536, 1125025312, 11689502240, 1125025312, 5973536, 7471, 2

Offset: 1

Roger L. Bagula, Dec 30 2008

			Triangle begins as:
   2;
  17,   17;
   2,  574,       2;
   2, 4911,    4911,          2;
   2, 5423,  156192,       5423,           2;
   2, 5935, 1413920,    1413920,        5935,          2;
   2, 6447, 2802720,   42656800,     2802720,       6447,       2;
   2, 6959, 4322592,  406009120,   406009120,    4322592,    6959,    2;
   2, 7471, 5973536, 1125025312, 11689502240, 1125025312, 5973536, 7471, 2;

G. C. Greubel, Rows n = 1..50 of the triangle, flattened

Cf. this sequence (j=7), A153653 (j=8), A153654 (j=9), A153655 (j=10).
Cf. A153516, A153518, A153520, A153521, A153648, A153649, A153650, A153651, A153656, A153657.
Cf. A001026 (powers of 17).

f:= func< n,j | Round(((3-(-1)^n)/2)*NthPrime(j)^(n-1) - 2^((3-(-1)^n)/2)) >;
function T(n,k,j)
  if n eq 2 then return NthPrime(j);
  elif (n eq 3 and k eq 2 or n eq 4 and k eq 2 or n eq 4 and k eq 3) then return f(n,j);
  elif (k eq 1 or k eq n) then return 2;
  else return T(n-1,k,j) + T(n-1,k-1,j) + (2*j+1)*NthPrime(j)*T(n-2,k-1,j);
  end if; return T;
end function;
[T(n,k,7): k in [1..n], n in [1..12]]; // G. C. Greubel, Mar 02 2021

T[n_, k_, j_]:= T[n,k,j]= If[n==2, Prime[j], If[n==3 && k==2 || n==4 && 2<=k<=3, ((3-(-1)^n)/2)*Prime[j]^(n-1) -2^((3-(-1)^n)/2), If[k==1 || k==n, 2, T[n-1,k,j] + T[n-1,k-1,j] + (2*j+1)*Prime[j]*T[n-2,k-1,j] ]]];
Table[T[n,k,7], {n,12}, {k,n}]//Flatten (* modified by G. C. Greubel, Mar 02 2021 *)

@CachedFunction
def f(n,j): return ((3-(-1)^n)/2)*nth_prime(j)^(n-1) - 2^((3-(-1)^n)/2)
def T(n,k,j):
    if (n==2): return nth_prime(j)
    elif (n==3 and k==2 or n==4 and 2<=k<=3): return f(n,j)
    elif (k==1 or k==n): return 2
    else: return T(n-1,k,j) + T(n-1,k-1,j) + (2*j+1)*nth_prime(j)*T(n-2,k-1,j)
flatten([[T(n,k,7) for k in (1..n)] for n in (1..12)]) # G. C. Greubel, Mar 02 2021

T(n, k, j) = T(n-1, k, j) + T(n-1, k-1, j) + (2*j + 1)*prime(j)*T(n-2, k-1, j) with T(2, k, j) = prime(j), T(3, 2, j) = 2*prime(j)^2 - 4, T(4, 2, j) = T(4, 3, j) = prime(j)^2 - 2, T(n, 1, j) = T(n, n, j) = 2 and j = 7.

Sum_{k=0..n} T(n, k, j) = 2*prime(j)^(n-1) for j=7 = 2*A001026(n-1).

Edited by G. C. Greubel, Mar 02 2021

A230622 Values of x such that x^2 + y^2 = 17^n with x and y coprime and 0 < x < y.

Original entry on oeis.org

1, 8, 47, 161, 404, 495, 2908, 31679, 203996, 905768, 2529647, 4839120, 4291039, 116593352, 859799153, 4896306240, 14978251196, 35359140465, 28242853388, 375162560801, 3481428994004, 21473668418415, 85367854683953, 228866364286560

Offset: 1

Colin Barker, Oct 26 2013

nonn

The corresponding y-values are in A230623.

			a(2)=8 because 8^2+15^2=289=17^2.

Vincenzo Librandi, Table of n, a(n) for n = 1..200
Chris Busenhart, Lorenz Halbeisen, Norbert Hungerbühler, Oliver Riesen, On primitive solutions of the Diophantine equation x^2+ y^2= M, Eidgenössische Technische Hochschule (ETH Zürich, Switzerland, 2020).

Cf. A001026, A188948, A230623.

Table[Select[PowersRepresentations[17^n, 2, 2], CoprimeQ@@#&][[1, 1]], {n, 1, 40}] (* Vincenzo Librandi, Mar 02 2014 *)

A230623 Values of y such that x^2 + y^2 = 17^n with x and y coprime and 0 < x < y.

Original entry on oeis.org

4, 15, 52, 240, 1121, 4888, 20047, 77280, 277441, 1093425, 5279468, 23647519, 99429196, 393425745, 1457109628, 4968639359, 24553864319, 113193708472, 488133974353, 1980778750800, 7547952442399, 26710380775592, 112605054449252

Offset: 1

Colin Barker, Oct 26 2013

nonn

The corresponding x-values are in A230622.

			a(2)=15 because 8^2 + 15^2 = 289 = 17^2.

Vincenzo Librandi, Table of n, a(n) for n = 1..200
Chris Busenhart, Lorenz Halbeisen, Norbert Hungerbühler, Oliver Riesen, On primitive solutions of the Diophantine equation x^2+ y^2= M, Eidgenössische Technische Hochschule (ETH Zürich, Switzerland, 2020).

Cf. A001026, A188949, A230622.

Table[Select[PowersRepresentations[17^n, 2, 2], CoprimeQ@@#&][[1, 2]], {n, 1, 40}] (* Vincenzo Librandi, Mar 02 2014 *)

cp's OEIS Frontend

A091045 Partial sums of powers of 17 (A001026).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

Sage

Formula

A060220 Number of orbits of length n under the full 17-shift (whose periodic points are counted by A001026).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

SageMath

Formula

Extensions

A008472 Sum of the distinct primes dividing n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Sage

Sage

Formula

A329697 a(n) is the number of iterations needed to reach a power of 2 starting at n and using the map k -> k-(k/p), where p is the largest prime factor of k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Formula

A003992 Square array read by upwards antidiagonals: T(n,k) = n^k for n >= 0, k >= 0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

Extensions

A329332 Table of powers of squarefree numbers, powers of A019565(n) in increasing order in row n. Square array A(n,k) n >= 0, k >= 0 read by descending antidiagonals.

Views

A153652 Triangle T(n, k, j) = T(n-1, k, j) + T(n-1, k-1, j) + (2j + 1)prime(j)*T(n-2, k-1, j) with T(2, k, j) = prime(j) and j = 7, read by rows.