wrong formula in A104350

[Richard Mathar]

It seems that one of the formulas given in A104350,
a(n)=A076928(n)*A000142(n),
is wrong, and only correct for the first 3 terms. The Maple functions to
compare the sequences are below. With a little bit of effort one could
probably figure out which product of which factors is actually missing...

This was detected while submitting a correction to A104362...

--Richard

A076928 := proc(n)
	local q ;
	if n = 1 then
		RETURN(1) ;
	else
		q := ifactors(n)[2] ;
		RETURN(A076928(n-1)*n/q[nops(q)][1]) ;
	fi ;
end:

A000142 := proc(n)
	RETURN(n!) ;
end:

A006530 := proc(n)
	local i, t1, t2, t3, t4, t5;
	if n = 1 then
		RETURN(1) ;
	else
		t1 := numtheory[divisors](n);
		t2 := convert(t1, list);
		t3 := sort(t2);
		t4 := nops(t3);
		t5 := 1;
		for i from 1 to t4 do
			if isprime(t3[t4+1-i]) then
				RETURN(t3[t4+1-i]);
			fi;
		od;
		RETURN(1);
	fi ;
end:

A104350 := proc(n)
	resul := 1 ;
	for k from 1 to n do
		resul := resul*A006530(k) ;
	od ;
	RETURN(resul) ;
end:

A104350wrong := proc(n)
	A076928(n)*A000142(n) ; # ????
end:

for n from 1 to 30 do
	#A006530(n) ;
	#A076928(n) ;
	print(n,A104350(n),A104350wrong(n)) ;
od ;