a(n) = q^2*(q^2+2*q-1)/2, where q = n-th prime power A000961(n).

Terms

a(0) =1a(1) =14a(2) =63a(3) =184a(4) =425a(5) =1519a(6) =2528a(7) =3969a(8) =8591a(9) =16393a(10) =36736a(11) =46529a(12) =71839a(13) =151823a(14) =210625a(15) =285039a(16) =377609a(17) =491071a(18) =556544a(19) =987049a(20) =1480961a(21) =1787983a(22) =2542559a(23) =2998849a(24) =4092713a(25) =6262319a(26) =7148041a(27) =8648704a(28) =10374079a(29) =13061231