A066092 Erroneous version of A063427.
2, 6, 4, 20, 3, 42, 8, 18, 10, 110, 4, 156, 14, 10, 16, 272, 18, 342, 5, 28, 22, 506, 8, 100
Offset: 2
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.
a(10) = 4 since 1/10 = 1/5 - 1/10 = 1/6 - 1/15 = 1/8 - 1/40 = 1/9 - 1/90. a(12) = 7: the divisors of 12 are 1, 2, 3, 4, 6 and 12 and the decompositions are (1, 2), (1, 3), (1, 4), (1, 6), (1, 12), (2, 3), (3, 4).
[(NumberOfDivisors(n^2)-1)/2 : n in [1..100]]; // Vincenzo Librandi, Apr 18 2018
Table[(Length[Divisors[n^2]] - 1)/2, {n, 1, 100}] (DivisorSigma[0,Range[100]^2]-1)/2 (* Harvey P. Dale, Apr 15 2013 *)
for(n=1,100,print1(sum(i=1,n^2,if((n*i)%(i+n),0,1)),","))
a(n)=numdiv(n^2)\2 \\ Charles R Greathouse IV, Oct 03 2016
Row 6 is (3,6,12,30) because 6+3 = 9 divides 6*3 = 18, 6+6 = 12 divides 6*6 = 36, 6+12 = 18 divides 6*12 = 72 and 6+30 = 36 divides 6*30 = 180.
for n from 2 to 20 do for m from 1 to n*(n-1) do if(m*n mod (m+n) = 0)then printf("%d, ",m): fi: od: od: # Nathaniel Johnston, Jun 22 2011
f[n_] := Select[Range[n^2], Mod[n*#, n + # ] == 0 &];Table[f[n], {n, 2, 24}] // Flatten (* Ray Chandler, Feb 13 2007 *)
arow(n)=local(d,m);d=divisors(n^2);vector(#d\2,k,m=d[ #d\2-k+1];n*(n-m)/m) \\ Franklin T. Adams-Watters, Aug 07 2009
a(6) = 2 because 6*3/(6+3) = 2 is the smallest integer of the form 6*k/(6+k). a(10) = 5 since 1/10 + 1/10 = 1/5, 1/10 + 1/15 = 1/6, 1/10 + 1/40 = 1/8, 1/10 + 1/90 = 1/9 and so the first sum provides the value.
spi[n_]:=Module[{k=1},While[!IntegerQ[(n*k)/(n+k)],k++];(n*k)/(n+k)]; Array[ spi,80,2] (* Harvey P. Dale, May 05 2022 *)
a(n)={my(k=1); if(n>1, while (n*k%(n + k), k++); n*k/(n + k))} \\ Harry J. Smith, Aug 20 2009
a(2) = 2 as 2*2 = 4 which contains 2 + 2 = 4 as a substring. a(4) = 68 as 4*68 = 272 which contains 4+68 = 72 as a substring. a(69) = 16961 as 69*16961 = 1170309 which contains 69+16961 = 17030 as a substring. a(501000) = 1002 as 501000*1002 = 502002000 which contains 501000+1002 = 502002 as a substring. This is the first of 500 consecutive terms with a(n) = 1002. a(554635) = 879948670 as 554635*879948670 = 488050330585450 which contain 554635+879948670 = 880503305 as a substring. This is the largest value of a(n) for the first one million terms.
isok(n, k) = #strsplit(Str(n*k), Str(n+k)) > 1; a(n) = {if (vecsearch([1, 3, 5, 6, 7, 9, 26], n), return (-1)); my(k=1); while (! isok(k, n), k++); k;} \\ Michel Marcus, Dec 02 2020 and Jan 23 2021
a(3) = 24 as 3*24 = 72 which contains reverse(3+24) = reverse(27) = 72 as a substring. a(6) = 34 as 6*34 = 204 which contains reverse(6+34) = reverse(40) = 04 as a substring. Note the leading zero is included. a(29) = 46716 as 29*46716 = 1354764 which contains reverse(29+4671) = reverse(46745) = 54764 as a substring. a(110) = 11 as 110*11 = 1210 which contains reverse(110+11) = reverse(121) = 121 as a substring. This is the first of two consecutive terms with a(n) = 11. a(20000) = 666843331 as 20000*666843331 = 13336866620000 which contains reverse(20000+666843331) = reverse(666863331) = 133368666 as a substring. This is the largest value in the first 100000 terms.
isok(n, k) = #strsplit(Str(n*k), concat(Vecrev(Str(n+k)))) > 1; ispt(n) = my(t); ispower(n,,&t) && (t==10); a(n) = {if ((n==1) || (n==10) || ispt(n), return (-1)); my(k=0); while (! isok(n, k), k++); k;} \\ Michel Marcus, Jan 22 2021
a(10)=15 since 1/10=1/20+1/20=1/30+1/15=1/35+1/14=1/60+1/12=1/110+1/11, but the first sum does not have c>b, leaving the second sum to provide the value.
f:= proc(n) local b; for b from 2*n-1 by -1 do if n*b mod (b-n) = 0 then return b fi od end proc: map(f, [$2..100]); # Robert Israel, Dec 01 2019
a[n_] := n + SelectFirst[Divisors[n^2] // Reverse, #Jean-François Alcover, Jun 07 2020 *)
a(2) = 6. The only candidate for a(2) that satisfies 3*a(2) being divisible by 3+a(2) is a(2) = 6. a(3) = 12. The candidates for a(3) given a(2) = 6 are 3,6,12,30, all of which satisfy 6*a(3) being divisible by 6+a(3). 3 and 6 have already appeared so the next smallest candidate is chosen, being 12. a(4) = 24. The candidates for a(4) given a(3) = 12 are 4,6,12,24,36,60,132, all of which satisfy 12*a(4) being divisible by 12+a(4). 4 is not allowed as there would be no candidates for a(5), and 6 and 12 have already appeared, so the next smallest candidate is chosen, being 24. a(5) = 8. There are 10 candidates for a(5) given a(4) = 24, the smallest that has not appeared is 8.
R:= 3: a:= 3: S:= {4,3}: for i from 2 to 100 do Cands:= remove (t -> t < 1 or isprime(t), map(`-`,numtheory:-divisors(a^2),a) minus S); a:= min(Cands); R:= R, a; S:= S union {a}; od: R; # Robert Israel, Mar 23 2023
a(10)=30 since 1/10=1/20+1/20=1/30+1/15=1/35+1/14=1/60+1/12=1/110+1/11, but the first sum does not have c>b, leaving the second sum to provide the value.
a(45)=75 because divisors of 45^2 are {1, 3, 5, 9, 15, 25, 27, 45, 75, 81, 135, 225, 405, 675, 2025} and the smallest element greater than 45 is 75.
with(numtheory): for n from 2 to 200 do a := divisors(n^2): b := a[(tau(n^2)-1)/2]: printf(`%d,`,n^2/b); od:
sdgn[n_]:=Select[Divisors[n^2],#>n&,1]; Flatten[Array[sdgn,70]] (* Harvey P. Dale, Jun 18 2012 *)
{ for (n=2, 1000, d=divisors(n^2); write("b063718.txt", n, " ", d[2 + length(d)\2]) ) } \\ Harry J. Smith, Aug 28 2009
Comments