#include <stdio.h>
#include <stdlib.h>

struct Element{
       char Name[900];
       int Nukleonenzahl[10];
       int Isoanzahl;
       int Zerfall[10]; //0=kein Zerfall, 1=alpha, 2= beta
       };

int main(int argc, char *argv[])
{
  struct Element Atom[200];
  int i;
  for(i=0;i<200;i++) Atom[i].Isoanzahl=0;
  
  strcpy(Atom[0].Name,"n");
  Atom[0].Isoanzahl=1;
  Atom[0].Nukleonenzahl[0]=1;
  Atom[0].Zerfall[0]=2;
  
  strcpy(Atom[1].Name,"H");
  Atom[1].Isoanzahl=3;
  Atom[1].Nukleonenzahl[0]=1;
  Atom[1].Nukleonenzahl[1]=2;
  Atom[1].Nukleonenzahl[2]=3;
  Atom[1].Zerfall[0]=0;
  Atom[1].Zerfall[1]=0;
  Atom[1].Zerfall[2]=2;
  
  strcpy(Atom[2].Name,"He");
  Atom[2].Isoanzahl=2;
  Atom[2].Nukleonenzahl[0]=3;
  Atom[2].Nukleonenzahl[1]=4;
  Atom[2].Zerfall[0]=0;
  Atom[2].Zerfall[1]=0;
  
  strcpy(Atom[3].Name,"Li");
  Atom[3].Isoanzahl=2;
  Atom[3].Nukleonenzahl[0]=6;
  Atom[3].Nukleonenzahl[1]=7;
  Atom[3].Zerfall[0]=0;
  Atom[3].Zerfall[1]=0;
  
  strcpy(Atom[4].Name,"Be");
  Atom[4].Isoanzahl=1;
  Atom[4].Nukleonenzahl[0]=9;
  Atom[4].Zerfall[0]=0;
  
  strcpy(Atom[5].Name,"B");
  Atom[5].Isoanzahl=2;
  Atom[5].Nukleonenzahl[0]=10;
  Atom[5].Nukleonenzahl[1]=11;
  Atom[5].Zerfall[0]=0;
  Atom[5].Zerfall[1]=0;
  
  strcpy(Atom[6].Name,"C");
  Atom[6].Isoanzahl=3;
  Atom[6].Nukleonenzahl[0]=12;
  Atom[6].Nukleonenzahl[1]=13;
  Atom[6].Nukleonenzahl[2]=14;
  Atom[6].Zerfall[0]=0;
  Atom[6].Zerfall[1]=0;
  Atom[6].Zerfall[2]=2;
  
  strcpy(Atom[7].Name,"N");
  Atom[7].Isoanzahl=2;
  Atom[7].Nukleonenzahl[0]=14;
  Atom[7].Nukleonenzahl[1]=15;
  Atom[7].Zerfall[0]=0;
  Atom[7].Zerfall[1]=0;
  
  strcpy(Atom[81].Name,"Tl");
  Atom[81].Isoanzahl=4;
  Atom[81].Nukleonenzahl[0]=203;
  Atom[81].Nukleonenzahl[1]=205;
  Atom[81].Nukleonenzahl[2]=207;
  Atom[81].Nukleonenzahl[3]=208;
  Atom[81].Zerfall[0]=0;
  Atom[81].Zerfall[1]=0;
  Atom[81].Zerfall[2]=2;
  Atom[81].Zerfall[3]=2;
  
  strcpy(Atom[82].Name,"Pb");
  Atom[82].Isoanzahl=8;
  Atom[82].Nukleonenzahl[0]=204;
  Atom[82].Nukleonenzahl[1]=206;
  Atom[82].Nukleonenzahl[2]=207;
  Atom[82].Nukleonenzahl[3]=208;
  Atom[82].Nukleonenzahl[4]=210;
  Atom[82].Nukleonenzahl[5]=211;
  Atom[82].Nukleonenzahl[6]=212;
  Atom[82].Nukleonenzahl[7]=214;
  Atom[82].Zerfall[0]=1;
  Atom[82].Zerfall[1]=0;
  Atom[82].Zerfall[2]=0;
  Atom[82].Zerfall[3]=0;
  Atom[82].Zerfall[4]=2;
  Atom[82].Zerfall[5]=2;
  Atom[82].Zerfall[6]=2;
  Atom[82].Zerfall[7]=2;
  
  strcpy(Atom[83].Name,"Bi");
  Atom[83].Isoanzahl=6;
  Atom[83].Nukleonenzahl[0]=209;
  Atom[83].Nukleonenzahl[1]=210;
  Atom[83].Nukleonenzahl[2]=211;
  Atom[83].Nukleonenzahl[3]=212;
  Atom[83].Nukleonenzahl[4]=214;
  Atom[83].Nukleonenzahl[5]=215;
  Atom[83].Zerfall[0]=0;
  Atom[83].Zerfall[1]=2;
  Atom[83].Zerfall[2]=1;
  Atom[83].Zerfall[3]=1;
  Atom[83].Zerfall[4]=2;
  Atom[83].Zerfall[5]=2;
  
  strcpy(Atom[84].Name,"Po");
  Atom[84].Isoanzahl=7;
  Atom[84].Nukleonenzahl[0]=210;
  Atom[84].Nukleonenzahl[1]=211;
  Atom[84].Nukleonenzahl[2]=212;
  Atom[84].Nukleonenzahl[3]=214;
  Atom[84].Nukleonenzahl[4]=215;
  Atom[84].Nukleonenzahl[5]=216;
  Atom[84].Nukleonenzahl[6]=218;
  Atom[84].Zerfall[0]=1;
  Atom[84].Zerfall[1]=1;
  Atom[84].Zerfall[2]=1;
  Atom[84].Zerfall[3]=1;
  Atom[84].Zerfall[4]=1;
  Atom[84].Zerfall[5]=1;
  Atom[84].Zerfall[6]=1;
  
  strcpy(Atom[85].Name,"At");
  Atom[85].Isoanzahl=1;
  Atom[85].Nukleonenzahl[0]=219;
  Atom[85].Zerfall[0]=1;
  
  strcpy(Atom[86].Name,"Rn");
  Atom[86].Isoanzahl=3;
  Atom[86].Nukleonenzahl[0]=219;
  Atom[86].Nukleonenzahl[1]=220;
  Atom[86].Nukleonenzahl[2]=222;
  Atom[86].Zerfall[0]=1;
  Atom[86].Zerfall[1]=1;
  Atom[86].Zerfall[2]=1;
  
  strcpy(Atom[87].Name,"Fr");
  Atom[87].Isoanzahl=1;
  Atom[87].Nukleonenzahl[0]=223;
  Atom[87].Zerfall[0]=2;
  
  strcpy(Atom[88].Name,"Ra");
  Atom[88].Isoanzahl=4;
  Atom[88].Nukleonenzahl[0]=223;
  Atom[88].Nukleonenzahl[1]=224;
  Atom[88].Nukleonenzahl[2]=226;
  Atom[88].Nukleonenzahl[3]=228;
  Atom[88].Zerfall[0]=1;
  Atom[88].Zerfall[1]=1;
  Atom[88].Zerfall[2]=1;
  Atom[88].Zerfall[3]=2;
  
  strcpy(Atom[89].Name,"Ac");
  Atom[89].Isoanzahl=2;
  Atom[89].Nukleonenzahl[0]=227;
  Atom[89].Nukleonenzahl[1]=228;
  Atom[89].Zerfall[0]=2;
  Atom[89].Zerfall[1]=2;
  
  strcpy(Atom[90].Name,"Th");
  Atom[90].Isoanzahl=6;
  Atom[90].Nukleonenzahl[0]=227;
  Atom[90].Nukleonenzahl[1]=228;
  Atom[90].Nukleonenzahl[2]=230;
  Atom[90].Nukleonenzahl[3]=231;
  Atom[90].Nukleonenzahl[4]=232;
  Atom[90].Nukleonenzahl[5]=234;
  Atom[90].Zerfall[0]=1;
  Atom[90].Zerfall[1]=1;
  Atom[90].Zerfall[2]=1;
  Atom[90].Zerfall[3]=2;
  Atom[90].Zerfall[4]=1;
  Atom[90].Zerfall[5]=2;
  
  strcpy(Atom[91].Name,"Pa");
  Atom[91].Isoanzahl=2;
  Atom[91].Nukleonenzahl[0]=231;
  Atom[91].Nukleonenzahl[1]=234;
  Atom[91].Zerfall[0]=1;
  Atom[91].Zerfall[1]=2;
  
  strcpy(Atom[92].Name,"U");
  Atom[92].Isoanzahl=3;
  Atom[92].Nukleonenzahl[0]=234;
  Atom[92].Nukleonenzahl[1]=235;
  Atom[92].Nukleonenzahl[2]=238;
  Atom[92].Zerfall[0]=1;
  Atom[92].Zerfall[1]=1;
  Atom[92].Zerfall[2]=1;
  
  int Ordnungszahl;
  int Nzahl;
  int a,b;
  char Wahl;
  
  do{
          printf("Ermittlung der Zerfallsreihen:\n\n");
          printf("Bitte geben Sie die Ordnungzahl des Elements an: ");
          scanf("%d",&Ordnungszahl);
          getchar();
          printf("Bitte geben Sie die Anzahl der Nukleonen des Elements an: ");
          scanf("%d",&Nzahl);
          getchar();
          
          printf("\n");
          
          if(Atom[Ordnungszahl].Isoanzahl==0) printf("Dieses Element kann nicht zerfallen.\n");
          else
          {
              a=0;
              b=0;
              do
              {
                        for(i=0;i<Atom[Ordnungszahl].Isoanzahl;i++) if(Atom[Ordnungszahl].Nukleonenzahl[i]==Nzahl)break; 
                        if(Atom[Ordnungszahl].Zerfall[i]==1)
                        {
                                   printf("Alpha-Zerfall:  %-2s%d -> He4 + ",Atom[Ordnungszahl].Name,Nzahl);
                                   Nzahl-=4;
                                   Ordnungszahl-=2;
                                   for(i=0;i<Atom[Ordnungszahl].Isoanzahl;i++) if(Atom[Ordnungszahl].Nukleonenzahl[i]==Nzahl)break; 
                                   printf("%-2s%d\n",Atom[Ordnungszahl].Name,Nzahl);  
                                   a++;                   
                        }  
                        if(Atom[Ordnungszahl].Zerfall[i]==2)
                        {
                                   printf("Beta-Zerfall:   %-2s%d -> e   + ",Atom[Ordnungszahl].Name,Nzahl);
                                   Ordnungszahl++;
                                   for(i=0;i<Atom[Ordnungszahl].Isoanzahl;i++) if(Atom[Ordnungszahl].Nukleonenzahl[i]==Nzahl)break; 
                                   printf("%-2s%d\n",Atom[Ordnungszahl].Name,Nzahl);   
                                   b++;                      
                        }
                        if(Atom[Ordnungszahl].Zerfall[i]==0) break;                             
              }while(1);
              
              printf("\nEs kam zu %d alpha Zerf\x84llen und zu %d beta Zerf\x84llen.\n\n",a,b);
          }
              
          printf("Wollen Sie erneut einen Zerfall berechnen(j/n): ");
          scanf("%c",&Wahl);
          getchar();
          system("cls");              
  }while(Wahl=='j');
  	
  return 0;
}