login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A268756
Number of triples x, y, r such that r divides x*y, r divides at least one of x or y, and x*y <= n.
1
1, 5, 9, 17, 21, 35, 39, 53, 61, 75, 79, 107, 111, 125, 139, 160, 164, 192, 196, 224, 238, 252, 256, 304, 312, 326, 340, 368, 372, 418, 422, 452, 466, 480, 494, 550, 554, 568, 582, 630, 634, 680, 684, 712, 740, 754, 758, 830, 838, 866, 880, 908, 912, 960, 974
OFFSET
1,2
COMMENTS
a(n) - a(n - 1) only depends on the prime signature of n. - David A. Corneth, Aug 30 2018
LINKS
David A. Corneth, PARI prog
Adrian W. Dudek, On the Success of Mishandling Euclid's Lemma, arXiv:1602.03555 [math.HO], 2016. See A(n) p. 2.
Adrian W. Dudek, On the Success of Mishandling Euclid's Lemma, The American Mathematical Monthly, Vol. 123, No. 9 (2016), 924-927.
PROG
(PARI) a(n) = {s = 0; for (x=1, n, for (y = 1, n, if (x*y <= n, s += sum(r = 1, x*y, !(x*y % r) && (!(x % r) || !(y % r))); ); ); ); s; }
(PARI) \\ See PARI link \\ David A. Corneth, Aug 30 2018
CROSSREFS
Sequence in context: A255651 A216877 A350989 * A095725 A005006 A369318
KEYWORD
nonn
AUTHOR
Michel Marcus, Feb 13 2016
STATUS
approved