cp's OEIS Frontend

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.

Previous Showing 31-40 of 80 results. Next

A376846 Number of m <= n such that rad(m) | n and Omega(m) > Omega(n), where rad = A007947 and Omega = A001222.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 2, 0, 0, 0, 2, 0, 1, 0, 2, 0, 0, 1, 3, 0, 1, 0, 3, 1, 1, 0, 4, 0, 2, 0, 3, 0, 0, 0, 3, 1, 2, 0, 2, 0, 1, 1, 3, 0, 2, 0, 3, 0, 0, 0, 7, 0, 3, 1, 5, 0, 1, 0, 4, 0, 3, 0, 8, 0, 1, 0, 4, 0, 4, 0, 4, 2
Offset: 1

Views

Author

Michael De Vlieger, Oct 06 2024

Keywords

Comments

Number of m not exceeding n such that the squarefree kernel of m divides n, and m has more prime factors with repetition than does n.
Number of m in row n of A162306 such that Omega(m) > Omega(n).

Examples

			Table of select n such that a(n) > 0:
   n  a(n)  List of m such that Omega(m) > Omega(n).
  -------------------------------------------------
  10   1    {8}
  14   1    {8}
  18   1    {16}
  20   1    {16}
  22   2    {8, 16}
  26   2    {8, 16}
  28   1    {16}
  30   2    {16, 24}
  33   1    {27}
  34   3    {8, 16, 32}
  36   1    {32}
  38   3    {8, 16, 32}
  39   1    {27}
  40   1    {32}
  42   4    {16, 24, 32, 36}
		

Crossrefs

Programs

  • Mathematica
    rad[x_] := rad[x] = Times @@ FactorInteger[x][[All, 1]];
    {0}~Join~Table[With[{k = PrimeOmega[n]}, Count[Range[n], _?(And[Divisible[n, rad[#]], PrimeOmega[#] > k] &)]], {n, 2, 120}]

Formula

a(n) = card({m <= n : rad(m) | n, Omega(m) > Omega(n) }).
a(n) = 0 for prime power n (in A000961).
a(n) < A010846(n).

A381800 a(n) = number of distinct residues r mod n of numbers k such that rad(k) | n, where rad = A007947.

Original entry on oeis.org

1, 2, 2, 3, 2, 5, 2, 4, 3, 7, 2, 8, 2, 6, 8, 5, 2, 12, 2, 9, 9, 13, 2, 11, 3, 15, 4, 9, 2, 19, 2, 6, 9, 11, 12, 16, 2, 21, 6, 12, 2, 24, 2, 16, 15, 14, 2, 16, 3, 28, 20, 17, 2, 31, 8, 12, 21, 31, 2, 28, 2, 8, 13, 7, 10, 32, 2, 13, 15, 35, 2, 20, 2, 39, 29, 24
Offset: 1

Views

Author

Michael De Vlieger, Mar 07 2025

Keywords

Examples

			 n  a(n)  row n of A381801
----------------------------------------------
 1    1   {0}
 2    2   {0,1}
 3    2   {0,1}
 4    3   {0,1,2}
 6    5   {0,1,2,3,4}
 8    4   {0,1,2,4}
10    7   {0,1,2,4,5,6,8}
12    8   {0,1,2,3,4,6,8,9}
14    6   {0,1,2,4,7,8}
15    8   {0,1,3,5,6,9,10,12}
18   12   {0,1,2,3,4,6,8,9,10,12,14,16}
20    9   {0,1,2,4,5,8,10,12,16}
21    9   {0,1,3,6,7,9,12,15,18}
22   13   {0,1,2,4,6,8,10,11,12,14,16,18,20}
24   11   {0,1,2,3,4,6,8,9,12,16,18}
26   15   {0,1,2,4,6,8,10,12,13,14,16,18,20,22,24}
28    9   {0,1,2,4,7,8,14,16,21}
30   19   {0,1,2,3,4,5,6,8,9,10,12,15,16,18,20,21,24,25,27}
36   16   {0,1,2,3,4,6,8,9,12,16,18,20,24,27,28,32}
		

Crossrefs

Programs

  • Mathematica
    Table[CountDistinct@ Flatten@ Mod[TensorProduct @@ Map[(p = #; NestWhileList[Mod[p*#, n] &, 1, UnsameQ, All]) &, FactorInteger[n][[All, 1]]], n], {n, 120}]

Formula

a(n) = length of row n of A381801.
a(1) = 1 since 1 is the empty product.
A010846(n) <= a(n) <= A051953(n).
a(n) >= 2 for n > 1.
For prime p, a(p) = A010846(p^m) = A000005(p^m) = A381798(p) = 2.
For prime power p^m, m > 0, a(p^m) = A010846(p^m) = A000005(p^m) = A381798(p^m) = m+1.
For n in A024619, a(n) > A381798(n).

A299992 Composite n with A001221(n) > 1 for which A243822(n) < A000005(n).

Original entry on oeis.org

6, 10, 12, 14, 15, 18, 20, 21, 22, 24, 26, 28, 33, 35, 36, 39, 40, 44, 45, 48, 51, 52, 55, 56, 57, 63, 65, 68, 69, 72, 75, 76, 77, 80, 85, 87, 88, 91, 92, 93, 95, 96, 99, 100, 104, 108, 111, 112, 115, 116, 117, 119, 123, 124, 129, 133, 135, 136, 141, 143, 144
Offset: 1

Views

Author

Michael De Vlieger, Feb 26 2018

Keywords

Comments

Composite numbers m have nondivisors k in the cototient such that k | n^e with e > 1. These k appear in row n of A272618 and are enumerated by A243822(n). These nondivisors k are a kind of "regular" number along with divisors d of n; both are listed in row n of A162306 and are together enumerated by A045763(n).
Primes p have 2 divisors {1, p}; these two numbers constitute the cototient of p: there are no nondivisors in the cototient.
Prime powers p^i have (i + 1) divisors; all smaller powers of the same prime p, i.e., p^j with 0 <= j <= i, also divide p^i. These numbers constitute the cototient of p^i; there are no nondivisors in the cototient.
Therefore, we can ignore cases where n has no nondivisors in the cototient, since they obviously have more divisors than nondivisors therein.
This sequence lists (composite) numbers n with omega(n) > 1 that have fewer nondivisors k in the cototient of n than divisors d.
The smallest odd term is 15.
The number m = 1001 is the smallest term with A001221(m) = 3. No term less than 36,000,000 has A001221(m) > 3.
The following terms m are the smallest to have A001222(m) = {2, 3, 4, ...}: {6, 12, 24, 48, 96, 192, 384, 1152, 2304, 4608, 13824, 27648, 55296, 110592, 331776, 663552, 1327104, 3981312, 7962624, 15925248, ...}
Number of terms less than 10^k for 0 <= k <= 7: {0, 2, 44, 319, 2171, 15545, 119469, 969749}.

Examples

			6 is the first term since it is the smallest number with more than one distinct prime divisor that has more divisors (4) than numbers in A243822(6) = 1.
		

Crossrefs

Programs

  • Mathematica
    Select[Range@ 144, Function[n, And[PrimeNu[n] > 1, Count[Range[n], _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)] < 2 DivisorSigma[0, n]]]]

A360543 a(n) = number of numbers k < n, gcd(k, n) > 1, such that omega(k) > omega(n) and rad(n) | rad(k), where omega(n) = A001221(n) and rad(n) = A007947(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 6, 0, 0, 0, 0, 11, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 2, 5, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 0, 0, 0, 1, 26, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 2, 0, 0, 0, 0, 3, 23, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 7, 0, 3, 1, 4
Offset: 1

Views

Author

Michael De Vlieger, Mar 06 2023

Keywords

Examples

			a(4) = 0 since k = 1..3 are prime powers.
a(8) = 1 since only k = 6 is such that p = 3, q = 5, but gcd(6, 10) = 2.
a(9) = 1 since the following satisfies definition: {6},
a(16) = 4, i.e., {6, 10, 12, 14},
a(25) = 3, i.e., {10, 15, 20},
a(27) = 6, i.e., {6, 12, 15, 18, 21, 24},
a(32) = 11, i.e., {6, 10, 12, 14, 18, 20, 22, 24, 26, 28, 30},
a(36) = 1, i.e., {30},
a(40) = 1, i.e., {30},
a(45) = 1, i.e., {30}, etc.
		

Crossrefs

Programs

  • Mathematica
    nn = 120; rad[n_] := rad[n] = Times @@ FactorInteger[n][[All, 1]]; c = Select[Range[4, nn], CompositeQ]; Table[Function[{q, r}, Count[TakeWhile[c, # <= n &], _?(And[PrimeNu[#] > q, Divisible[rad[#], r]] &)]] @@ {PrimeNu[n], rad[n]}, {n, nn}]

Formula

a(n) = A243823(n) - A360480(n).
a(n) = A045763(n) - A243822(n) - A360480(n).
a(n) = A051953(n) - A000005(n) - A243822(n) - A360480(n).
a(n) = A051953(n) - A010846(n) - A360480(n).
a(n) = A243823(n) = A045763(n) for n in A246547.
For prime power n = p^e, n > 1, a(n) = p^(e-1) - e.
For n in A360765, a(n) > 0.

A208815 n for which A079277(n) + phi(n) < n.

Original entry on oeis.org

115, 329, 1243, 2119, 2171, 4709, 4777, 4811, 6593, 6631, 6707, 6821, 11707, 11983, 12029, 14597, 15463, 16793, 23809, 23867, 23983, 24041, 24331, 29047, 29171, 29357, 29543, 50357, 50579, 67937, 68183, 68347, 68429, 77873, 78389, 78733, 79421, 83351, 83453, 102413
Offset: 1

Views

Author

Robert Israel, Mar 01 2012

Keywords

Comments

Includes (among other terms, see below) semiprimes pq where p and q are primes with p^k-p+1 < q < p^k for an integer k>1. In particular, by the Prime Number Theorem this sequence is infinite. - clarified by Antti Karttunen, Apr 26 2017
From Antti Karttunen, Apr 26 2017: (Start)
Numbers n for which A051953(n) > A079277(n).
Factorization of terms a(1) .. a(29): 5*23, 7*47, 11*113, 13*163, 13*167, 17*277, 17*281, 17*283, 19*347, 19*349, 19*353, 19*359, 23*509, 23*521, 23*523, 11*1327, 7*47*47, 7*2399, 29*821, 29*823, 29*827, 29*829, 29*839, 31*937, 31*941, 31*947, 31*953, 37*1361, 37*1367. Note that a(17) = 15463 is not a semiprime.
(End)

Examples

			A079277(115) + phi(115) = 25 + 88 = 113 < 115 so 115 is in the sequence, where phi = A000010.
		

Crossrefs

Positions of negative terms in A285709.

Programs

  • Mathematica
    Select[Range[2, 10^4], Function[n, If[n == 2, 1, Module[{k = n - 2, e = Floor@ Log2@ n}, While[PowerMod[n, e, k] != 0, k--]; k]] + EulerPhi@ n < n]] (* or *)
    Do[If[If[n == 2, 1, Module[{k = n - 2, e = Floor@ Log2@ n}, While[PowerMod[n, e, k] != 0, k--]; k]] + EulerPhi@ n < n, Print@ n], {n, 2, 10^5}] (* Michael De Vlieger, Apr 27 2017 *)

Extensions

a(28)-a(29) from Antti Karttunen, Apr 26 2017
a(30)-a(40) from David A. Corneth, Apr 26 2017

A275280 Irregular triangle listing numbers m of n that have prime divisors p that also divide n, in order of appearance in a matrix of products that arranges the powers of prime divisors p of n along independent axes.

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 2, 4, 1, 5, 1, 2, 4, 3, 6, 1, 7, 1, 2, 4, 8, 1, 3, 9, 1, 2, 4, 8, 5, 10, 1, 11, 1, 2, 4, 8, 3, 6, 12, 9, 1, 13, 1, 2, 4, 8, 7, 14, 1, 3, 9, 5, 15, 1, 2, 4, 8, 16, 1, 17, 1, 2, 4, 8, 16, 3, 6, 12, 9, 18, 1, 19, 1, 2, 4, 8, 16, 5, 10, 20, 1, 3, 9, 7, 21, 1, 2, 4, 8, 16, 11, 22, 1, 23
Offset: 1

Views

Author

Michael De Vlieger, Jul 28 2016

Keywords

Comments

Product matrix of tensors T = 1,p,p^2,...,p^e that include the powers 1 <= e of prime divisors p such that p^e <= n.
This sequence is analogous to A275055 but differs from it in that the tensors T include not only powers p^e that divide n but all powers p^e <= n.
The matrix a(n) is bounded by n, thus all products m <= n.
Let omega(n) = A001221(n). The matrix a(n) has omega(n) dimensions and is an omega(n)-dimensional simplex with (omega(n)-1) "right-angle: sides and 1 irregular surface that is bounded by n.
A027750(n) is a subset of A162306(n) and in a(n), the terms of A275055(n) appear in an contiguous omega(n)-dimensional parallelepiped (parallelotope) with 1 at the origin and n at the opposite corner. Thus the omega(n)-dimensional array described by A275055(n) is fully contained in the simplex-like matrix described by a(n). Divisors appear within the parallelepiped while nondivisors appear in the field outside the parallelepiped (see examples). Terms within the parallelepiped appear in A027750(n) while those outside appear in A272618(n).
For a(2^x + 2) there is a term m = (n-2); m != (n-1) except for n=2, since GCD(n, n-1)=1.
a(p^e) = A027750(p^e) = A162306(p^e) = A275055(p^e) for e >= 1.

Examples

			Triangle begins:
1;
1, 2;
1, 3;
1, 2, 4;
1, 5;
1, 2, 4, 3, 6;
1, 7;
1, 2, 4, 8;
1, 3, 9;
1, 2, 4, 8, 5, 10;
1, 11;
1, 2, 4, 8, 3, 6, 12, 9;
1, 13;
1, 2, 4, 8, 7, 14;
1, 3, 9, 5, 15;
1  2, 4, 8, 16;
1, 17;
1, 2, 4, 8, 16, 3, 6, 12, 9, 18;
...
2 prime divisors: n = 96
   1  2  4  8 16 32 64
   3  6 12 24 48 96
   9 18 36 72
  27 54
  81
thus a(96) = {1,2,4,8,16,32,64,3,6,12,24,48,96,9,18,36,72,27,54,81}.
The divisors of 72 (thus the terms of A275055(72)) appear in a rectangle delimited by 1 at top left and 72 at bottom right.
3 prime divisors: n = 60
(the 3 dimensional levels correspond with powers of 5)
   level 5^0:            level 5^1:         level 5^2:
   1  2  4  8 16 32  |    5 10 20 40    |   25 50
   3  6 12 24 48     |   15 30 60       |
   9 18 36           |   45             |
  27 54              |                  |
thus a(60) = {1,2,4,8,16,32,3,6,12,24,48,9,18,36,27,54,5,10,20,40,15,30,60,45,25,50}.
The divisors of 60 (thus the terms of A275055(60)) appear in a parallelepiped delimited by 1 at top left of level 5^0 and 60 at bottom right of level 5^1.
		

Crossrefs

Cf. A162306, A010846 (row length), A243103 (row product), A027750 (divisors of n), A000005 (number of divisors of n), A272618 (nondivisors m <= n that have prime divisors p that also divide n), A243822 (number of such nondivisors of n), A275055 (Product of tensor of prime divisor powers that are also divisors).

Programs

  • Mathematica
    f[n_] := If[n == 1, 1, Function[w, ToExpression@ StringJoin["With[{n=", ToString@ n, "}, Table[", ToString@ InputForm[Times @@ Map[Power @@ # &, w]], ", ", Most@ Flatten@ Map[{#, ", "} &, #], "]]"] &@ MapIndexed[Function[p, StringJoin["{", ToString@ Last@ p, ", 0, Log[", ToString@ First@ p, ", n/(", ToString@ InputForm[Times @@ Map[Power @@ # &, Take[w, First@ #2 - 1]]], ")]}"]]@ w[[First@ #2]] &, w]]@ Map[{#, ToExpression["p" <> ToString@ PrimePi@ #]} &, Reverse[FactorInteger[n][[All, 1]]]] ]; Array[f, 24] // Flatten (* Michael De Vlieger, Mar 08 2017 *)

A279907 Triangle read by rows: T(n,k) is the smallest power of n that is divisible by k, or -1 if no such power exists.

Original entry on oeis.org

0, 0, 1, 0, -1, 1, 0, 1, -1, 1, 0, -1, -1, -1, 1, 0, 1, 1, 2, -1, 1, 0, -1, -1, -1, -1, -1, 1, 0, 1, -1, 1, -1, -1, -1, 1, 0, -1, 1, -1, -1, -1, -1, -1, 1, 0, 1, -1, 2, 1, -1, -1, 3, -1, 1, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 0, 1, 1, 1, -1, 1, -1, 2, 2, -1, -1, 1, 0, -1, -1, -1, -1, -1, -1
Offset: 1

Views

Author

Michael De Vlieger, Dec 26 2016

Keywords

Comments

T(n,1) = 0 since 1 | n^0.
T(n,p) = 1 for prime divisors p of n since p | n^1.
T(n,d) = 1 for divisors d > 1 of n since d | n^1.
Row n for prime p have maximum value 1, since all k < p are coprime to p, and k | p^1 only when k = p.
Values greater than 1 pertain only to composite k of composite n > 4, but not in all cases. T(n,k) = 1 for squarefree kernels k of composite n.
T(n,k) = -1 for numbers k > 1 coprime to n and for numbers that are products of at least one prime q coprime to n and one prime p | n.
T(n,k) is nonnegative for all numbers k for which n^k (mod k) = 0, i.e., all the prime divisors p of k also divide n.
The largest possible value s in row n of T = floor(log_2(n)), since the largest possible multiplicity of any number m <= n pertains to perfect powers of 2, as 2 is the smallest prime. This number s first appears at T(2^s + 2, 2^s) for s > 1.
If T(n,k) is positive, 1/k terminates T(n,k) digits after the radix point in base n. If T(n,k) is negative, 1/k is recurrent in base n.
From Robert Israel, Dec 28 2016: (Start)
T(a*b,c*d) = max(T(a,c),T(b,d)) if GCD(a,b)=1, GCD(b,d)=1,T(a,c)>=0 and T(b,d)>=0.
T(n,a*b) = max(T(n,a),T(n,b)) if GCD(a,b)=1 and T(n,a)>=0 and T(n,b)>=0. (End)

Examples

			The triangle T(n,k) begins (with -1 shown as "." for clarity):
n\k 1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 ...
1:  0
2:  0  1
3:  0  .  1
4:  0  1  .  1
5:  0  .  .  .  1
6:  0  1  1  2  .  1
7:  0  .  .  .  .  .  1
8:  0  1  .  1  .  .  .  1
9:  0  .  1  .  .  .  .  .  1
10: 0  1  .  2  1  .  .  3  .  1
11: 0  .  .  .  .  .  .  .  .  .  1
12: 0  1  1  1  .  1  .  2  2  .  .  1
13: 0  .  .  .  .  .  .  .  .  .  .  .  1
14: 0  1  .  2  .  .  1  3  .  .  .  .  .  1
15: 0  .  1  .  1  .  .  .  2  .  .  .  .  .  1
...
		

Crossrefs

Cf.: A010846 (number of nonnegative k in row n), A162306 (k with nonnegative values in a(n)), A051731 (k with values 0 or 1), A000005 (number of k in row n with values 0 or 1), A272618 (k with values > 1), A243822 (number of k in row n with values > 1), A007947.

Programs

  • Maple
    f:= proc(n,k) local Fk,Fn,i;
       if k = 1 then return 0 fi;
       Fk:= ifactors(k)[2];
       Fn:= map(t -> padic:-ordp(n,t[1]),Fk);
       if min(Fn) = 0 then -1 else max(seq(ceil(Fk[i,2]/Fn[i]),i=1..nops(Fk))) fi
    end proc:
    seq(seq(f(n,k),k=1..n),n=1..20); # Robert Israel, Dec 28 2016
  • Mathematica
    Table[Boole[k == 1] + (Boole[#[[-1, 1]] == 1] (-1 + Length@ #) /. 0 -> -1) &@ NestWhileList[Function[s, {#1/s, s}]@ GCD[#1, #2] & @@ # &, {k, n}, And[First@ # != 1, ! CoprimeQ @@ #] &], {n, 16}, {k, n}] // Flatten (* or *)
    Table[SelectFirst[Range[0, Floor@ Log2@ n], PowerMod[n, #, k] == 0 &] /. k_ /; MissingQ@ k -> -1, {n, 12}, {k, n}] // TableForm (* Version 10.2 *)

A280269 Irregular triangle T(n,m) read by rows: smallest power e of n that is divisible by m = term k in row n of A162306.

Original entry on oeis.org

0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 2, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 2, 1, 3, 1, 0, 1, 0, 1, 1, 1, 1, 2, 2, 1, 0, 1, 0, 1, 2, 1, 3, 1, 0, 1, 1, 2, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 1, 2, 1, 3, 1, 2, 4, 1, 0, 1, 0, 1, 1, 1, 2, 1, 2, 1, 0, 1, 1, 2, 1, 0, 1, 2, 3, 1, 4, 1, 0, 1, 0, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Michael De Vlieger, Dec 30 2016

Keywords

Comments

This table eliminates the negative values in row n of A279907.
Let k = A162306(n,m), i.e., the value in column m of row n.
T(n,1) = 0 since 1 | n^0.
T(n,p) = 1 for prime divisors p of n since p | n^1.
T(n,d) = 1 for divisors d > 1 of n since d | n^1.
Row n for prime p have two terms, {0,1}, the maximum value 1, since all k < p are coprime to p, and k | p^1 only when k = p.
Row n for prime power p^i have (i+1) terms, one zero and i ones, since all k that appear in corresponding row n of A162306 are divisors d of p^i.
Values greater than 1 pertain only to composite k of composite n > 4, but not in all cases. T(n,k) = 1 for squarefree kernels k of composite n.
Numbers k > 1 coprime to n and numbers that are products of at least one prime q coprime to n and one prime p | n do not appear in A162306; these do not divide n^e evenly.
T(n,k) is nonnegative for all numbers k for which n^k (mod k) = 0, i.e., all the prime divisors p of k also divide n.
The largest possible value s in row n of T = floor(log_2(n)), since the largest possible multiplicity of any number m <= n pertains to perfect powers of 2, as 2 is the smallest prime. This number s first appears at T(2^s + 2, 2^s) for s > 1.
1/k terminates T(n,k) digits after the radix point in base n for values of k that appear in row n of A162306.
Originally from Robert Israel at A279907: (Start)
T(a*b,c*d) = max(T(a,c),T(b,d)) if GCD(a,b)=1, GCD(b,d)=1,T(a,c)>=0 and T(b,d)>=0.
T(n,a*b) = max(T(n,a),T(n,b)) if GCD(a,b)=1 and T(n,a)>=0 and T(n,b)>=0.
(End)

Examples

			Triangle T(n,m) begins:  Triangle A162306(n,k):
1:  0                    1
2:  0  1                 1  2
3:  0  1                 1  3
4:  0  1  1              1  2  4
5:  0  1                 1  5
6:  0  1  1  2  1        1  2  3  4  6
7:  0  1                 1  7
8:  0  1  1  1           1  2  4  8
9:  0  1  1              1  3  9
10: 0  1  2  1  3  1     1  2  4  5  8  10
...
		

Crossrefs

Cf. A162306, A279907 (T(n,k) with values for all 1 <= k <= n), A280274 (maximum values in row n), A010846 (number of nonnegative k in row n), A051731 (k with e <= 1), A000005 (number of k in row n with e <= 1), A272618 (k with e > 1), A243822 (number of k in row n with e > 1), A007947.

Programs

  • Mathematica
    Table[SelectFirst[Range[0, #], PowerMod[n, #, k] == 0 &] /. m_ /; MissingQ@ m -> Nothing &@ Floor@ Log2@ n, {n, 24}, {k, n}] // Flatten (* Version 10.2, or *)
    DeleteCases[#, -1] & /@ Table[If[# == {}, -1, First@ #] &@ Select[Range[0, #], PowerMod[n, #, k] == 0 &] &@ Floor@ Log2@ n, {n, 24}, {k, n}] // Flatten (* or *)
    DeleteCases[#, -1] & /@ Table[Boole[k == 1] + (Boole[#[[-1, 1]] == 1] (-1 + Length@ #) /. 0 -> -1) &@ NestWhileList[Function[s, {#1/s, s}]@ GCD[#1, #2] & @@ # &, {k, n}, And[First@# != 1, ! CoprimeQ @@ #] &], {n, 24}, {k, n}] // Flatten

A300155 Numbers n for which A243822(n) = A000005(n).

Original entry on oeis.org

34, 38, 46, 50, 54, 58, 62, 105, 249, 267, 268, 284, 291, 292, 303, 309, 316, 321, 324, 327, 332, 339, 356, 363, 381, 385, 388, 393, 404, 411, 412, 417, 428, 436, 447, 452, 453, 455, 471, 484, 489, 500, 501, 507, 508, 519, 537, 543, 573, 579, 591, 595, 597
Offset: 1

Views

Author

Michael De Vlieger, Feb 26 2018

Keywords

Comments

Indices of zeros in A299990, i.e., A010846(n) - 2*A000005(n) = 0.
Composite numbers m have nondivisors k in the cototient such that k | n^e with e > 1. These k appear in row n of A272618 and are enumerated by A243822(n). These nondivisors k are a kind of "regular" number along with divisors d of n; both are listed in row n of A162306 and are together enumerated by A045763(n). Divisors of n are listed in row n of A027750.
This sequence lists numbers that have an equal number of nondivisors k in the cototient of n as divisors d.
The smallest odd term is 105.

Examples

			34 is the first term since it is the smallest number for which A243822(34) = A000005(34). For n = 34, there are 4 divisors {1, 2, 17, 34} and 4 nondivisors 1 <= m <= n such that m | n^e with e > 1: {4, 8, 16, 32}.
		

Crossrefs

Programs

  • Mathematica
    Select[Range@ 600, Function[n, Count[Range[n], _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)] == 2 DivisorSigma[0, n]]]

A300156 Indices of records in A299990.

Original entry on oeis.org

1, 30, 42, 66, 78, 90, 102, 114, 138, 150, 210, 330, 390, 510, 570, 630, 870, 990, 1050, 1470, 1890, 2100, 2310, 2730, 3570, 3990, 4620, 5460, 6510, 6930, 8190, 9240, 10710, 11550, 13650, 16170, 19110, 20790, 23100, 24570, 25410, 30030, 39270, 43890, 46410
Offset: 1

Views

Author

Michael De Vlieger, Feb 26 2018

Keywords

Comments

A010846(n) = A000005(n) + A243822(n).
Successive terms in this sequence represent increasing differences A243822(n) - A000005(n).
A000079 = records in -1 * A299990, since A243822(p^e)=0 for e>=0, n = 2^k sets records in A000005(n). The corresponding records are in A000027.

Examples

			The first term is 1: A299990(1) = -1. The first term that exceeds -1 in A299990 corresponds to n = 30; A299990(30) = 2, thus 30 is the second term.
		

Crossrefs

Programs

  • Mathematica
    With[{s = Table[Count[Range[n], _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)] - 2 DivisorSigma[0, n], {n, 10^3}]}, Map[FirstPosition[s, #][[1]] &, Union@ FoldList[Max, s]]]
Previous Showing 31-40 of 80 results. Next