cp's OEIS Frontend

A063198 Dimension of the space of weight 2n cuspidal newforms for Gamma_0( 10 ).

%I A063198 #16 Jul 19 2024 11:10:47
%S A063198 0,1,3,1,3,5,3,5,7,5,7,9,7,9,11,9,11,13,11,13,15,13,15,17,15,17,19,17,
%T A063198 19,21,19,21,23,21,23,25,23,25,27,25,27,29,27,29,31,29,31,33,31,33
%N A063198 Dimension of the space of weight 2n cuspidal newforms for Gamma_0( 10 ).
%C A063198 The dimension of weight n is apparently given by 0, 0, 2, 1, 0, 3, 2, 1, 4,... etc as in A063942. - _R. J. Mathar_, Jul 14 2015
%H A063198 R. J. Mathar, <a href="/A063198/b063198.txt">Table of n, a(n) for n = 1..1000</a>
%H A063198 William A. Stein, <a href="http://wstein.org/Tables/dimskg0new.gp">Dimensions of the spaces S_k^{new}(Gamma_0(N))</a>
%H A063198 William A. Stein, <a href="http://wstein.org/Tables/">The modular forms database</a>
%H A063198 <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (1,0,1,-1).
%F A063198 G.f.: x^2*(1+2*x-2*x^2+x^3) / ( (1+x+x^2)*(x-1)^2 ). - _R. J. Mathar_, Jul 15 2015
%F A063198 For n>1, a(n) = (6*n-3+12*cos(2*n*Pi/3)-4*sqrt(3)*sin(2*n*Pi/3))/9. - _Wesley Ivan Hurt_, Sep 30 2017
%p A063198 s0star := proc(n)
%p A063198     local pf,a,p,e ;
%p A063198     if n = 1 then
%p A063198         1;
%p A063198     else
%p A063198         a :=1 ;
%p A063198         for pf in ifactors(n)[2] do
%p A063198             p := op(1,pf) ;
%p A063198             e := op(2,pf) ;
%p A063198             if e =1 then
%p A063198                 a := a*(1-1/p) ;
%p A063198             elif e = 2 then
%p A063198                 a := a*(1-1/p-1/p^2) ;
%p A063198             else
%p A063198                 a := a*(1-1/p)*(1-1/p^2) ;
%p A063198             end if;
%p A063198         end do:
%p A063198         a ;
%p A063198     end if;
%p A063198 end proc:
%p A063198 nuInfstar := proc(n)
%p A063198     local pf,a,p,e ;
%p A063198     if n = 1 then
%p A063198         1;
%p A063198     else
%p A063198         a :=1 ;
%p A063198         for pf in ifactors(n)[2] do
%p A063198             p := op(1,pf) ;
%p A063198             e := op(2,pf) ;
%p A063198             if type(e,'odd') then
%p A063198                 return 0;
%p A063198             elif e = 2 then
%p A063198                 a := a*(p-2) ;
%p A063198             else
%p A063198                 a := a*(p-1)^2*p^(e/2-2) ;
%p A063198             end if;
%p A063198         end do:
%p A063198         a ;
%p A063198     end if;
%p A063198 end proc:
%p A063198 nu2star := proc(n)
%p A063198     local pf,a,p,e ;
%p A063198     if n = 1 then
%p A063198         1;
%p A063198     else
%p A063198         a :=1 ;
%p A063198         for pf in ifactors(n)[2] do
%p A063198             p := op(1,pf) ;
%p A063198             e := op(2,pf) ;
%p A063198             if p = 2 then
%p A063198                 if e =1 or e =2  then
%p A063198                     a := -a ;
%p A063198                 elif e =3 then
%p A063198                     ;
%p A063198                 else
%p A063198                     return 0 ;
%p A063198                 end if;
%p A063198             elif modp(p,4) = 1 then
%p A063198                 if e = 2 then
%p A063198                     a := -a ;
%p A063198                 else
%p A063198                     return 0;
%p A063198                 end if;
%p A063198             else
%p A063198                 if e = 1 then
%p A063198                     a := -2*a ;
%p A063198                 elif e = 2 then
%p A063198                     ;
%p A063198                 else
%p A063198                     return 0;
%p A063198                 end if;
%p A063198             end if;
%p A063198         end do:
%p A063198         a ;
%p A063198     end if;
%p A063198 end proc:
%p A063198 nu3star := proc(n)
%p A063198     local pf,a ;
%p A063198     if n = 1 then
%p A063198         1;
%p A063198     else
%p A063198         a :=1 ;
%p A063198         for pf in ifactors(n)[2] do
%p A063198             p := op(1,pf) ;
%p A063198             e := op(2,pf) ;
%p A063198             if p = 3 then
%p A063198                 if e =1 or e =2  then
%p A063198                     a := -a ;
%p A063198                 elif e =3 then
%p A063198                     ;
%p A063198                 else
%p A063198                     return 0 ;
%p A063198                 end if;
%p A063198             elif modp(p,3) = 1 then
%p A063198                 if e = 2 then
%p A063198                     a := -a ;
%p A063198                 else
%p A063198                     return 0;
%p A063198                 end if;
%p A063198             else
%p A063198                 if e = 1 then
%p A063198                     a := -2*a ;
%p A063198                 elif e = 2 then
%p A063198                     ;
%p A063198                 else
%p A063198                     return 0;
%p A063198                 end if;
%p A063198             end if;
%p A063198         end do:
%p A063198         a ;
%p A063198     end if;
%p A063198 end proc:
%p A063198 c2 := proc(k)
%p A063198     1/4+floor(k/4)-k/4 ;
%p A063198 end proc:
%p A063198 c3 := proc(k)
%p A063198     1/3+floor(k/3)-k/3 ;
%p A063198 end proc:
%p A063198 g0star := proc(k,N)
%p A063198     local a;
%p A063198     a := (k-1)/12*N*s0star(N) -nuInfstar(N)/2 +c2(k)*nu2star(N)+c3(k)*nu3star(N) ;
%p A063198     if k/2 = 1 then
%p A063198         a := a+numtheory[mobius](N) ;
%p A063198     end if;
%p A063198     a;
%p A063198 end proc:
%p A063198 A063198 := proc(n)
%p A063198     g0star(2*n,10) ;
%p A063198 end proc:
%p A063198 A063199 := proc(n)
%p A063198     g0star(2*n,11) ;
%p A063198 end proc:
%p A063198 A063200 := proc(n)
%p A063198     g0star(2*n,15) ;
%p A063198 end proc:
%p A063198 A063201 := proc(n)
%p A063198     g0star(2*n,18) ;
%p A063198 end proc:
%p A063198 A063205 := proc(n)
%p A063198     g0star(2*n,29) ;
%p A063198 end proc: # _R. J. Mathar_, Jul 19 2024
%Y A063198 Cf. A063942.
%K A063198 nonn,easy
%O A063198 1,3
%A A063198 _N. J. A. Sloane_, Jul 10 2001