[seqfan] Re: Uphi(m)=UPhi(n)=3/4*(m*n)^(1/2)

Richard Mathar mathar at strw.leidenuniv.nl
Tue Nov 4 21:42:27 CET 2008


yk> From seqfan-bounces at list.seqfan.eu  Tue Nov  4 08:37:51 2008
yk> Date: Tue, 04 Nov 2008 15:20:21 +0900
yk> From: "koh" <zbi74583.boat at orange.zero.jp>
yk> To: seqfan at list.seqfan.eu
yk> Subject: [seqfan]  Uphi(m)=UPhi(n)=3/4*(m*n)^(1/2)
yk> 
yk>     Hi, Seqfans
yk>     Definition of Amicable Number is the following.
yk> 
yk>     Sigma(x) = Sigma(y) = x+y
yk> 
yk>     If "Sigma" is replaced with the other divisor functions or "x+y" is replaced with high degree or rational or irrational formulas then it becomes generalized Amicable Number.
yk> 
yk>     I calculated the examples of the following equtation.    
yk> 
yk>     UnitaryPhi(x) = UnitaryPhi(y) = 3/4*(x*y)^(1/2) , y<=x
yk> 
yk>     x=y=2^2
yk>     x=y=2^3*7
yk>     x=y=2^4*5
yk>     x=y=2^5*5*31
yk>     x=y=2^6*3^2*7
yk>     x=y=2^6*3^3*7*13
yk>     x=y=2^7*3^2*7*127
yk>     x=y=2^7*3^3*7*13*127
yk>     x=y=2^8*5*17          
yk>     x=y=2^9*3^2*7*73
yk>     x=y=2^9*3^3*7*13*73
yk>     x=y=2^10*3^3*5^2*11*13*31
yk>                
yk>     x=5*3^2*11^2,      y=5*31^2 
yk>     x=2^11*5*3^2*89^2, y=2^11*5*11^2*23^2
yk> 
yk>     Could anyone do exhausitive search?
yk>     I would like to know the other examples of the case not{x=y}.
yk> 
yk> 
yk> 
yk> 
yk>     %I A000001
yk>     %S A000001 4, 56, 80, 4032, 4960, 5445, 21760, 157248, 1024128, 39940992, 2354688, 91832832, 729999360
yk>     %N A000001 Numbers m,n such that UnitaryPhi(m)=UnitaryPhi(m)=3/4*(m*n)^(1/2), n<=m

There is an obvious typo in this definition (?!). I am also strongly in favor of including
explicit OEIS A-number references to rare functions such as UnitaryPhi(.) into
submissions of that kind (see the discussion in the seqfan list of Oct 2006).

I ran the following Maple prog and confirm that this list is
complete in the range y<=x<=21760. There are no further values
 <= 21760 in the sequence. A Java version is available on request.

A047994 := proc(n) option remember ;
	local ifacts;
	ifacts := ifactors(n)[2] ;
	mul( op(1,op(i,ifacts))^op(2,op(i,ifacts))-1,i=1..nops(ifacts)) ;
end:
x := 1:
while true do
	uphix := A047994(x) ;
	if uphix mod 3 = 0 then
		for y from 1 to x do
			uphiy := A047994(y) ;
			if uphix^2 = uphiy^2 and uphix^2 = 9/16*x*y then
				print(x,ifactors(x)[2],y,ifactors(y)[2]) ;
			fi ;
		od :
	fi:
	x := x+1 :
od :

yk>     %C A000001 a(6) and a(13) are the cases not{n=m}  
yk>     %e A000001 Factorization of cases not{n=m}
yk>                m,n=5*3^2*11^2,5*31^2
yk>                m.n=2^11*5*3^2*89^2,2^11*5*11^2*23^2
yk>     %Y A000001 A000002
yk>     %K A000001 none
yk>     %O A000001 0,1
yk>     %A A000001 Yasutsohi Kohmoto zbi74583.boat at orange.zero.jp
yk> 
yk>     %I A000002
yk>     %S A000002 4, 56, 80, 4032, 4960, 4805, 21760, 157248, 1024128, 39940992, 2354688, 91832832, 655452160
yk> 
yk>     %N A000002 Numbers m,n such that UnitaryPhi(m)=UnitaryPhi(m)=3/4*(m*n)^(1/2), n<=m

There is the obvious typo in the definition again. Of course it is puzzling that
the definition is the same as A00001 above, although the numbers are different (?!).
(BTW, I know the answers to the questions I am asking.)

yk>     %C A000002 a(6) and a(13) are the cases not{n=m}  
yk>     %Y A000002 A000001
yk>     %K A000002 none
yk>     %O A000002 0,1
yk>     %A A000002 Yasutsohi Kohmoto zbi74583.boat at orange.zero.jp
-- 
Richard J Mathar
The Netherlands                 URL http://www.strw.leidenuniv.nl/~mathar




More information about the SeqFan mailing list