OFFSET
2,1
COMMENTS
The entries of the sequence are the solution of the differential equation m'=n/n'*m.
If n=n' then m'=m (A051674) and m=p^p, with p prime. Taking the minimum prime p=2, m=4.
If n'|n then m'=b*m, where b=n/n', and m=2^(2*b).
In general n and n’ can have same common factor. Let us consider the fraction k/h obtained by reduction of n/n' [GCD(k,h)=1]. If h=p, p prime, then if k>h m=4*p^(k-h) else if k<h m=p^k.
If h is composite let us make the variable substitution m=h*t. Now, m'=h'*t+h*t'=(k/h)*h*t=k*t that leads to t'=1/h*(k-h')*t. Let us reduce the fraction (k-h')/h, if necessary. If we get a prime denominator we can solve the equation by the previous formulas. Otherwise we must iterate the variable substitutions until we get:
1) a prime denominator.
2) a differential equation of the type y'=0, whose solution is y=1.
3) a differential equation of the type y'=x*y, with x<0. In this case there is no solution. This corresponds to the zeros in the sequence (12, 14, 15, 28, 35, 39, 44, 46, 50, 51, 55, 65, etc.).
In cases 1) and 2), moving backward through all the substitutions, we reach the final solution.
LINKS
Paolo P. Lava, Table of n, a(n) for n = 2..100
EXAMPLE
n=18; n'=11. m'=18/11*m. Here n>n', n'=p and therefore m=4*11^7.
n=24; n'=44. m'=24/44*m. Simplifying the fraction we have m'=6/11*m. Here k<k' e k'=p and therefore m=116.
n=44; n'=48. m'=44/48*m. Simplifying the fraction we have m'=11/12*m. Let m=12t, m’=16t+12t'=11t. We have t'=-5/12*t that has no solutions. Therefore a(44)=0.
n=62; n'=33. m'=62/33*m. Let m=33*t, m'=14*t+33*t'=62*t. We have t’=48/33*t=16/11*t. therefore t=4*11^5 and m=33*4*11^5.
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Apr 20 2012
STATUS
approved