#!/usr/bin/python
import sys,os

infile=sys.stdin
outfile=sys.stdout
if len(sys.argv)>1:
    infile=open(sys.argv[1],"r")
if len(sys.argv)>2:
    outfile=open(sys.argv[2],"w")

anames=[ "H","He", "Li","Be","B","C","N","O","F","Ne",
         "Na","Mg","Al","Si","P","S","Cl","Ar",
         "K","Ca","Sc","Ti","V","Cr","Mn","Fe","Co","Ni","Cu","Zn","Ga","Ge","As","Se","Br","Kr",
         "Rb","Sr","Y","Zr","Nb","Mo","Tc","Ru","Rh","Pd","Ag","Cd","In","Sn","Sb","Te","I","Xe",
         "Cs","Ba",
         "La","Ce","Pr","Nd","Pm","Sm","Eu","Gd","Tb","Dy","Ho","Er","Tm","Yb","Lu",
         "Hf","Ta","W","Re","Os","Ir","Pt","Au","Hg","Tl","Pb","Bi","Po","At","Rn",
         "Fr","Ra",
         "Ac","Th","Pa","U","Np","Pu","Am","Cm","Bk","Cf","Es","Fm","Md","No","Lr",
         "Rf","Db","Sg","Bh","Hs","Mt" ]

btoa=0.529177

line=infile.readline()
atoms=[]
if line[:6]=="$coord":
    while 1:
        line=infile.readline()
        if not line:
            break
        if line[0]=="$":
            break
        (x,y,z,a)=line.split()
        A=a.capitalize()
        anr=0
        for i in range(len(anames)):
            if A==anames[i]:
                anr=i+1
                break
        if not anr:
            sys.stderr.write("Cannot determine atom number for atom "+A+"\n")
            sys.exit(1)
        atoms.append((float(x)*btoa,float(y)*btoa,float(z)*btoa,anr))
    outfile.write("1\n"+`len(atoms)`+"\n")
    for x in atoms:
        outfile.write(`x[3]`+" "+`x[0]`+" "+`x[1]`+" "+`x[2]`+"\n")
else:
    sys.stderr.write("Input is not a turbomole coord file!\n")
    sys.exit(1)
