program irrep
c  Calculates the number of repetitions of L for a given tau.
c  I know (now, dec '03) it might have been much shorter...have a look at
c  irhar.for, for a shorter example. 
c  Written by L.Fortunato, Feb. '03 and Jul. '03.
        implicit integer(a-z)
        parameter(MAXD=1000,MAXTAU=50,maxn=500)
        dimension vec(0:MAXTAU,MAXD),conta(0:maxn)


        do tau=0,MAXTAU
           c=0
           do i=1,MAXD
              vec(tau,i)=-10
           enddo
           do n=0,maxn
              conta(n)=0
           enddo 
           numax=int(tau/3)
           do nu=0,numax
              K=tau-3*nu
              do J=K,2*K
                 if (J.eq.2*K-1) goto 100
                 c=c+1 
                 vec(tau,c)=J
 100          enddo
           enddo
c           print *,c
           cmax=c
           do n=0,maxn
              do c=1,cmax  
                 if (vec(tau,c).eq.n) then
                           conta(n)=conta(n)+1                   
                 endif  
              enddo
           enddo
           sum=0
           do n=0,maxn
              if (conta(n).gt.0) sum=sum+conta(n)             
           enddo
           print *,tau,' ',sum
        enddo



        end