? for(n=1,10^8,if(Mod(sigma(n),n+eulerphi(n))==0,print(n," ",sigma(n)/(n+eulerphi(n))," ",factor(n))))
2 1 Mat([2, 1])
456 2 [2, 3; 3, 1; 19, 1]
828 2 [2, 2; 3, 2; 23, 1]
7584 2 [2, 5; 3, 1; 79, 1]
33462 2 [2, 1; 3, 2; 11, 1; 13, 2]
1357440 3 [2, 7; 3, 1; 5, 1; 7, 1; 101, 1]
1596048 2 [2, 4; 3, 1; 41, 1; 811, 1]
1964544 2 [2, 9; 3, 1; 1279, 1]