/* package whatever; // don't place package name! */
/*
    
    This source code accompanies the article, "Using The Golay Error Detection And
    Correction Code", by Hank Wallace. This program demonstrates use of the Golay
    code.
      Usage: G DATA Encode/Correct/Verify/Test
        where DATA is the data to be encoded, codeword to be corrected,
        or codeword to be checked for errors. DATA is hexadecimal.
      Examples:
        G 555 E      encodes information value 555 and prints a codeword
        G ABC123 C   corrects codeword ABC123
        G ABC123 V   checks codeword ABC123 for errors
        G ABC123 T   tests routines, ABC123 is a dummy parameter
    This program may be freely incorporated into your programs as needed. It
    compiles under Borland's Turbo C 2.0. No warranty of any kind is granted.
    */
    #include "stdio.h"
    #include "conio.h"
    #define POLY  0xAE3  /* or use the other polynomial, 0xC75 */
 
    /* ====================================================== */
 
    unsigned long golay(unsigned long cw)
    /* This function calculates [23,12] Golay codewords.
      The format of the returned longint is
      [checkbits(11),data(12)]. */
    {
      int i;
      unsigned long c;
      cw&=0xfffl;
      c=cw; /* save original codeword */
      for (i=1; i<=12; i++)  /* examine each data bit */
        {
          if (cw & 1)        /* test data bit */
            cw^=POLY;        /* XOR polynomial */
          cw>>=1;            /* shift intermediate result */
        }
      return((cw<<12)|c);    /* assemble codeword */
    }
 
    /* ====================================================== */
 
    int parity(unsigned long cw)
    /* This function checks the overall parity of codeword cw.
      If parity is even, 0 is returned, else 1. */
    {
      unsigned char p;
 
      /* XOR the bytes of the codeword */
      p=*(unsigned char*)&cw;
      p^=*((unsigned char*)&cw+1);
      p^=*((unsigned char*)&cw+2);
 
      /* XOR the halves of the intermediate result */
      p=p ^ (p>>4);
      p=p ^ (p>>2);
      p=p ^ (p>>1);
 
      /* return the parity result */
      return(p & 1);
    }
 
    /* ====================================================== */
 
    unsigned long syndrome(unsigned long cw)
    /* This function calculates and returns the syndrome
      of a [23,12] Golay codeword. */
    {
      int i;
      cw&=0x7fffffl;
      for (i=1; i<=12; i++)  /* examine each data bit */
        {
          if (cw & 1)        /* test data bit */
            cw^=POLY;        /* XOR polynomial */
          cw>>=1;            /* shift intermediate result */
        }
      return(cw<<12);        /* value pairs with upper bits of cw */
    }
 
    /* ====================================================== */
 
    int weight(unsigned long cw)
    /* This function calculates the weight of
      23 bit codeword cw. */
    {
      int bits,k;
 
      /* nibble weight table */
      const char wgt[16] = {0,1,1,2,1,2,2,3,1,2,2,3,2,3,3,4};
 
      bits=0; /* bit counter */
      k=0;
      /* do all bits, six nibbles max */
      while ((k<6) && (cw))
        {
          bits=bits+wgt[cw & 0xf];
          cw>>=4;
          k++;
        }
 
      return(bits);
    }
 
    /* ====================================================== */
 
    unsigned long rotate_left(unsigned long cw, int n)
    /* This function rotates 23 bit codeword cw left by n bits. */
    {
      int i;
 
      if (n != 0)
        {
          for (i=1; i<=n; i++)
            {
              if ((cw & 0x400000l) != 0)
                cw=(cw << 1) | 1;
              else
                cw<<=1;
            }
        }
 
      return(cw & 0x7fffffl);
    }
 
    /* ====================================================== */
 
    unsigned long rotate_right(unsigned long cw, int n)
    /* This function rotates 23 bit codeword cw right by n bits. */
    {
      int i;
 
      if (n != 0)
        {
          for (i=1; i<=n; i++)
            {
              if ((cw & 1) != 0)
                cw=(cw >> 1) | 0x400000l;
              else
                cw>>=1;
            }
        }
 
      return(cw & 0x7fffffl);
    }
 
    /* ====================================================== */
 
    unsigned long correct(unsigned long cw, int *errs)
    /* This function corrects Golay [23,12] codeword cw, returning the
      corrected codeword. This function will produce the corrected codeword
      for three or fewer errors. It will produce some other valid Golay
      codeword for four or more errors, possibly not the intended
      one. *errs is set to the number of bit errors corrected. */
    {
      unsigned char
        w;                /* current syndrome limit weight, 2 or 3 */
      unsigned long
        mask;             /* mask for bit flipping */
      int
        i,j;              /* index */
      unsigned long
        s,                /* calculated syndrome */
        cwsaver;          /* saves initial value of cw */
 
      cwsaver=cw;         /* save */
      *errs=0;
      w=3;                /* initial syndrome weight threshold */
      j=-1;               /* -1 = no trial bit flipping on first pass */
      mask=1;
      while (j<23) /* flip each trial bit */
        {
          if (j != -1) /* toggle a trial bit */
            {
              if (j>0) /* restore last trial bit */
                {
                  cw=cwsaver ^ mask;
                  mask+=mask; /* point to next bit */
                }
              cw=cwsaver ^ mask; /* flip next trial bit */
              w=2; /* lower the threshold while bit diddling */
            }
 
          s=syndrome(cw); /* look for errors */
          if (s) /* errors exist */
            {
              for (i=0; i<23; i++) /* check syndrome of each cyclic shift */
                {
                  if ((*errs=weight(s)) <= w) /* syndrome matches error pattern */
                    {
                      cw=cw ^ s;              /* remove errors */
                      cw=rotate_right(cw,i);  /* unrotate data */
                      return(s=cw);
                    }
                  else
                    {
                      cw=rotate_left(cw,1);   /* rotate to next pattern */
                      s=syndrome(cw);         /* calc new syndrome */
                    }
                }
              j++; /* toggle next trial bit */
            }
          else
            return(cw); /* return corrected codeword */
        }
 
      return(cwsaver); /* return original if no corrections */
    } /* correct */
 
    /* ====================================================== */
 
    int decode(int correct_mode, int *errs, unsigned long *cw)
    /* This function decodes codeword *cw in one of two modes. If correct_mode
      is nonzero, error correction is attempted, with *errs set to the number of
      bits corrected, and returning 0 if no errors exist, or 1 if parity errors
      exist. If correct_mode is zero, error detection is performed on *cw,
      returning 0 if no errors exist, 1 if an overall parity error exists, and
      2 if a codeword error exists. */
    {
      unsigned long parity_bit;
 
      if (correct_mode)               /* correct errors */
        {
          parity_bit=*cw & 0x800000l; /* save parity bit */
          *cw&=~0x800000l;            /* remove parity bit for correction */
 
          *cw=correct(*cw, errs);     /* correct up to three bits */
          *cw|=parity_bit;            /* restore parity bit */
 
          /* check for 4 bit errors */
          if (parity(*cw))            /* odd parity is an error */
            return(1);
          return(0); /* no errors */
        }
      else /* detect errors only */
        {
          *errs=0;
          if (parity(*cw)) /* odd parity is an error */
            {
              *errs=1;
              return(1);
            }
          if (syndrome(*cw))
            {
              *errs=1;
              return(2);
            }
          else
            return(0); /* no errors */
        }
    } /* decode */
 
    /* ====================================================== */
 
    void golay_test(void)
    /* This function tests the Golay routines for detection and correction
      of various patterns of error_limit bit errors. The error_mask cycles
      over all possible values, and error_limit selects the maximum number
      of induced errors. */
    {
      unsigned long
        error_mask,         /* bitwise mask for inducing errors */
        trashed_codeword,   /* the codeword for trial correction */
        virgin_codeword;    /* the original codeword without errors */
      unsigned char
        pass=1,             /* assume test passes */
        error_limit=3;      /* select number of induced bit errors here */
      int
        error_count;        /* receives number of errors corrected */
 
      virgin_codeword=golay(0x555); /* make a test codeword */
      if (parity(virgin_codeword))
        virgin_codeword^=0x800000l;
      for (error_mask=0; error_mask<0x800000l; error_mask++)
        {
          /* filter the mask for the selected number of bit errors */
          if (weight(error_mask) <= error_limit) /* you can make this faster! */
            {
              trashed_codeword=virgin_codeword ^ error_mask; /* induce bit errors */
 
              decode(1,&error_count,&trashed_codeword); /* try to correct bit errors */
 
              if (trashed_codeword ^ virgin_codeword)
                {
                  printf("Unable to correct %d errors induced with error mask = 0x%lX\n",
                    weight(error_mask),error_mask);
                  pass=0;
                }
 
              if (kbhit()) /* look for user input */
                {
                  if (getch() == 27) return; /* escape exits */
 
                  /* other key prints status */
                  printf("Current test count = %ld of %ld\n",error_mask,0x800000l);
                }
            }
        }
      printf("Golay test %s!\n",pass?"PASSED":"FAILED");
    }
 
    /* ====================================================== */
 
    void main(int argument_count, char *argument[])
    {
      int i,j;
      unsigned long l,g;
      const char *errmsg = "Usage: G DATA Encode/Correct/Verify/Test\n\n"
                 "  where DATA is the data to be encoded, codeword to be corrected,\n"
                 "  or codeword to be checked for errors. DATA is hexadecimal.\n\n"
                 "Examples:\n\n"
                 "  G 555 E      encodes information value 555 and prints a codeword\n"
                 "  G ABC123 C   corrects codeword ABC123\n"
                 "  G ABC123 V   checks codeword ABC123 for errors\n"
                 "  G ABC123 T   tests routines, ABC123 is a dummy parameter\n\n";
 
      if (argument_count != 3)
        {
          printf(errmsg);
          exit(0);
        }
 
      if (sscanf(argument[1],"%lx",&l) != 1)
        {
          printf(errmsg);
          exit(0);
        }
 
      switch (toupper(*argument[2]))
        {
          case 'E': /* encode */
            l&=0xfff;
            l=golay(l);
            if (parity(l)) l^=0x800000l;
            printf("Codeword = %lX\n",l);
            break;
 
          case 'V': /* verify */
            if (decode(0,&i,&l))
              printf("Codeword %lX is not a Golay codeword.\n",l);
            else
              printf("Codeword %lX is a Golay codeword.\n",l);
            break;
 
          case 'C': /* correct */
            g=l; /* save initial codeword */
            j=decode(1,&i,&l);
            if ((j) && (i))
              printf("Codeword %lX had %d bits corrected,\n"
                     "resulting in codeword %lX with a parity error.\n",g,i,l);
            else
              if ((j == 0) && (i))
                printf("Codeword %lX had %d bits corrected, resulting in codeword %lX.\n",g,i,l);
              else
                if ((j) && (i == 0))
                  printf("Codeword %lX has a parity error. No bits were corrected.\n",g);
                else
                  if ((j == 0) && (i == 0))
                    printf("Codeword %lX does not require correction.\n",g);
            break;
 
          case 'T': /* test */
            printf("Press SPACE for status, ESC to exit test...\n");
            golay_test();
            break;
 
          default:
            printf(errmsg);
            exit(0);
        }
    }
 
    /* end of G.C */
 
 
import java.util.*;
import java.lang.*;
import java.io.*;
 
/* Name of the class has to be "Main" only if the class is public. */
class Ideone
{
	public static void main (String[] args) throws java.lang.Exception
	{
		// your code goes here
	}
}

LyogcGFja2FnZSB3aGF0ZXZlcjsgLy8gZG9uJ3QgcGxhY2UgcGFja2FnZSBuYW1lISAqLwovKgogICAgCiAgICBUaGlzIHNvdXJjZSBjb2RlIGFjY29tcGFuaWVzIHRoZSBhcnRpY2xlLCAiVXNpbmcgVGhlIEdvbGF5IEVycm9yIERldGVjdGlvbiBBbmQKICAgIENvcnJlY3Rpb24gQ29kZSIsIGJ5IEhhbmsgV2FsbGFjZS4gVGhpcyBwcm9ncmFtIGRlbW9uc3RyYXRlcyB1c2Ugb2YgdGhlIEdvbGF5CiAgICBjb2RlLgogICAgICBVc2FnZTogRyBEQVRBIEVuY29kZS9Db3JyZWN0L1ZlcmlmeS9UZXN0CiAgICAgICAgd2hlcmUgREFUQSBpcyB0aGUgZGF0YSB0byBiZSBlbmNvZGVkLCBjb2Rld29yZCB0byBiZSBjb3JyZWN0ZWQsCiAgICAgICAgb3IgY29kZXdvcmQgdG8gYmUgY2hlY2tlZCBmb3IgZXJyb3JzLiBEQVRBIGlzIGhleGFkZWNpbWFsLgogICAgICBFeGFtcGxlczoKICAgICAgICBHIDU1NSBFICAgICAgZW5jb2RlcyBpbmZvcm1hdGlvbiB2YWx1ZSA1NTUgYW5kIHByaW50cyBhIGNvZGV3b3JkCiAgICAgICAgRyBBQkMxMjMgQyAgIGNvcnJlY3RzIGNvZGV3b3JkIEFCQzEyMwogICAgICAgIEcgQUJDMTIzIFYgICBjaGVja3MgY29kZXdvcmQgQUJDMTIzIGZvciBlcnJvcnMKICAgICAgICBHIEFCQzEyMyBUICAgdGVzdHMgcm91dGluZXMsIEFCQzEyMyBpcyBhIGR1bW15IHBhcmFtZXRlcgogICAgVGhpcyBwcm9ncmFtIG1heSBiZSBmcmVlbHkgaW5jb3Jwb3JhdGVkIGludG8geW91ciBwcm9ncmFtcyBhcyBuZWVkZWQuIEl0CiAgICBjb21waWxlcyB1bmRlciBCb3JsYW5kJ3MgVHVyYm8gQyAyLjAuIE5vIHdhcnJhbnR5IG9mIGFueSBraW5kIGlzIGdyYW50ZWQuCiAgICAqLwogICAgI2luY2x1ZGUgInN0ZGlvLmgiCiAgICAjaW5jbHVkZSAiY29uaW8uaCIKICAgICNkZWZpbmUgUE9MWSAgMHhBRTMgIC8qIG9yIHVzZSB0aGUgb3RoZXIgcG9seW5vbWlhbCwgMHhDNzUgKi8KICAgIAogICAgLyogPT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09ICovCiAgICAKICAgIHVuc2lnbmVkIGxvbmcgZ29sYXkodW5zaWduZWQgbG9uZyBjdykKICAgIC8qIFRoaXMgZnVuY3Rpb24gY2FsY3VsYXRlcyBbMjMsMTJdIEdvbGF5IGNvZGV3b3Jkcy4KICAgICAgVGhlIGZvcm1hdCBvZiB0aGUgcmV0dXJuZWQgbG9uZ2ludCBpcwogICAgICBbY2hlY2tiaXRzKDExKSxkYXRhKDEyKV0uICovCiAgICB7CiAgICAgIGludCBpOwogICAgICB1bnNpZ25lZCBsb25nIGM7CiAgICAgIGN3Jj0weGZmZmw7CiAgICAgIGM9Y3c7IC8qIHNhdmUgb3JpZ2luYWwgY29kZXdvcmQgKi8KICAgICAgZm9yIChpPTE7IGk8PTEyOyBpKyspICAvKiBleGFtaW5lIGVhY2ggZGF0YSBiaXQgKi8KICAgICAgICB7CiAgICAgICAgICBpZiAoY3cgJiAxKSAgICAgICAgLyogdGVzdCBkYXRhIGJpdCAqLwogICAgICAgICAgICBjd149UE9MWTsgICAgICAgIC8qIFhPUiBwb2x5bm9taWFsICovCiAgICAgICAgICBjdz4+PTE7ICAgICAgICAgICAgLyogc2hpZnQgaW50ZXJtZWRpYXRlIHJlc3VsdCAqLwogICAgICAgIH0KICAgICAgcmV0dXJuKChjdzw8MTIpfGMpOyAgICAvKiBhc3NlbWJsZSBjb2Rld29yZCAqLwogICAgfQogICAgCiAgICAvKiA9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT0gKi8KICAgIAogICAgaW50IHBhcml0eSh1bnNpZ25lZCBsb25nIGN3KQogICAgLyogVGhpcyBmdW5jdGlvbiBjaGVja3MgdGhlIG92ZXJhbGwgcGFyaXR5IG9mIGNvZGV3b3JkIGN3LgogICAgICBJZiBwYXJpdHkgaXMgZXZlbiwgMCBpcyByZXR1cm5lZCwgZWxzZSAxLiAqLwogICAgewogICAgICB1bnNpZ25lZCBjaGFyIHA7CiAgICAKICAgICAgLyogWE9SIHRoZSBieXRlcyBvZiB0aGUgY29kZXdvcmQgKi8KICAgICAgcD0qKHVuc2lnbmVkIGNoYXIqKSZjdzsKICAgICAgcF49KigodW5zaWduZWQgY2hhciopJmN3KzEpOwogICAgICBwXj0qKCh1bnNpZ25lZCBjaGFyKikmY3crMik7CiAgICAKICAgICAgLyogWE9SIHRoZSBoYWx2ZXMgb2YgdGhlIGludGVybWVkaWF0ZSByZXN1bHQgKi8KICAgICAgcD1wIF4gKHA+PjQpOwogICAgICBwPXAgXiAocD4+Mik7CiAgICAgIHA9cCBeIChwPj4xKTsKICAgIAogICAgICAvKiByZXR1cm4gdGhlIHBhcml0eSByZXN1bHQgKi8KICAgICAgcmV0dXJuKHAgJiAxKTsKICAgIH0KICAgIAogICAgLyogPT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09ICovCiAgICAKICAgIHVuc2lnbmVkIGxvbmcgc3luZHJvbWUodW5zaWduZWQgbG9uZyBjdykKICAgIC8qIFRoaXMgZnVuY3Rpb24gY2FsY3VsYXRlcyBhbmQgcmV0dXJucyB0aGUgc3luZHJvbWUKICAgICAgb2YgYSBbMjMsMTJdIEdvbGF5IGNvZGV3b3JkLiAqLwogICAgewogICAgICBpbnQgaTsKICAgICAgY3cmPTB4N2ZmZmZmbDsKICAgICAgZm9yIChpPTE7IGk8PTEyOyBpKyspICAvKiBleGFtaW5lIGVhY2ggZGF0YSBiaXQgKi8KICAgICAgICB7CiAgICAgICAgICBpZiAoY3cgJiAxKSAgICAgICAgLyogdGVzdCBkYXRhIGJpdCAqLwogICAgICAgICAgICBjd149UE9MWTsgICAgICAgIC8qIFhPUiBwb2x5bm9taWFsICovCiAgICAgICAgICBjdz4+PTE7ICAgICAgICAgICAgLyogc2hpZnQgaW50ZXJtZWRpYXRlIHJlc3VsdCAqLwogICAgICAgIH0KICAgICAgcmV0dXJuKGN3PDwxMik7ICAgICAgICAvKiB2YWx1ZSBwYWlycyB3aXRoIHVwcGVyIGJpdHMgb2YgY3cgKi8KICAgIH0KICAgIAogICAgLyogPT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09ICovCiAgICAKICAgIGludCB3ZWlnaHQodW5zaWduZWQgbG9uZyBjdykKICAgIC8qIFRoaXMgZnVuY3Rpb24gY2FsY3VsYXRlcyB0aGUgd2VpZ2h0IG9mCiAgICAgIDIzIGJpdCBjb2Rld29yZCBjdy4gKi8KICAgIHsKICAgICAgaW50IGJpdHMsazsKICAgIAogICAgICAvKiBuaWJibGUgd2VpZ2h0IHRhYmxlICovCiAgICAgIGNvbnN0IGNoYXIgd2d0WzE2XSA9IHswLDEsMSwyLDEsMiwyLDMsMSwyLDIsMywyLDMsMyw0fTsKICAgIAogICAgICBiaXRzPTA7IC8qIGJpdCBjb3VudGVyICovCiAgICAgIGs9MDsKICAgICAgLyogZG8gYWxsIGJpdHMsIHNpeCBuaWJibGVzIG1heCAqLwogICAgICB3aGlsZSAoKGs8NikgJiYgKGN3KSkKICAgICAgICB7CiAgICAgICAgICBiaXRzPWJpdHMrd2d0W2N3ICYgMHhmXTsKICAgICAgICAgIGN3Pj49NDsKICAgICAgICAgIGsrKzsKICAgICAgICB9CiAgICAKICAgICAgcmV0dXJuKGJpdHMpOwogICAgfQogICAgCiAgICAvKiA9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT0gKi8KICAgIAogICAgdW5zaWduZWQgbG9uZyByb3RhdGVfbGVmdCh1bnNpZ25lZCBsb25nIGN3LCBpbnQgbikKICAgIC8qIFRoaXMgZnVuY3Rpb24gcm90YXRlcyAyMyBiaXQgY29kZXdvcmQgY3cgbGVmdCBieSBuIGJpdHMuICovCiAgICB7CiAgICAgIGludCBpOwogICAgCiAgICAgIGlmIChuICE9IDApCiAgICAgICAgewogICAgICAgICAgZm9yIChpPTE7IGk8PW47IGkrKykKICAgICAgICAgICAgewogICAgICAgICAgICAgIGlmICgoY3cgJiAweDQwMDAwMGwpICE9IDApCiAgICAgICAgICAgICAgICBjdz0oY3cgPDwgMSkgfCAxOwogICAgICAgICAgICAgIGVsc2UKICAgICAgICAgICAgICAgIGN3PDw9MTsKICAgICAgICAgICAgfQogICAgICAgIH0KICAgIAogICAgICByZXR1cm4oY3cgJiAweDdmZmZmZmwpOwogICAgfQogICAgCiAgICAvKiA9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT0gKi8KICAgIAogICAgdW5zaWduZWQgbG9uZyByb3RhdGVfcmlnaHQodW5zaWduZWQgbG9uZyBjdywgaW50IG4pCiAgICAvKiBUaGlzIGZ1bmN0aW9uIHJvdGF0ZXMgMjMgYml0IGNvZGV3b3JkIGN3IHJpZ2h0IGJ5IG4gYml0cy4gKi8KICAgIHsKICAgICAgaW50IGk7CiAgICAKICAgICAgaWYgKG4gIT0gMCkKICAgICAgICB7CiAgICAgICAgICBmb3IgKGk9MTsgaTw9bjsgaSsrKQogICAgICAgICAgICB7CiAgICAgICAgICAgICAgaWYgKChjdyAmIDEpICE9IDApCiAgICAgICAgICAgICAgICBjdz0oY3cgPj4gMSkgfCAweDQwMDAwMGw7CiAgICAgICAgICAgICAgZWxzZQogICAgICAgICAgICAgICAgY3c+Pj0xOwogICAgICAgICAgICB9CiAgICAgICAgfQogICAgCiAgICAgIHJldHVybihjdyAmIDB4N2ZmZmZmbCk7CiAgICB9CiAgICAKICAgIC8qID09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PSAqLwogICAgCiAgICB1bnNpZ25lZCBsb25nIGNvcnJlY3QodW5zaWduZWQgbG9uZyBjdywgaW50ICplcnJzKQogICAgLyogVGhpcyBmdW5jdGlvbiBjb3JyZWN0cyBHb2xheSBbMjMsMTJdIGNvZGV3b3JkIGN3LCByZXR1cm5pbmcgdGhlCiAgICAgIGNvcnJlY3RlZCBjb2Rld29yZC4gVGhpcyBmdW5jdGlvbiB3aWxsIHByb2R1Y2UgdGhlIGNvcnJlY3RlZCBjb2Rld29yZAogICAgICBmb3IgdGhyZWUgb3IgZmV3ZXIgZXJyb3JzLiBJdCB3aWxsIHByb2R1Y2Ugc29tZSBvdGhlciB2YWxpZCBHb2xheQogICAgICBjb2Rld29yZCBmb3IgZm91ciBvciBtb3JlIGVycm9ycywgcG9zc2libHkgbm90IHRoZSBpbnRlbmRlZAogICAgICBvbmUuICplcnJzIGlzIHNldCB0byB0aGUgbnVtYmVyIG9mIGJpdCBlcnJvcnMgY29ycmVjdGVkLiAqLwogICAgewogICAgICB1bnNpZ25lZCBjaGFyCiAgICAgICAgdzsgICAgICAgICAgICAgICAgLyogY3VycmVudCBzeW5kcm9tZSBsaW1pdCB3ZWlnaHQsIDIgb3IgMyAqLwogICAgICB1bnNpZ25lZCBsb25nCiAgICAgICAgbWFzazsgICAgICAgICAgICAgLyogbWFzayBmb3IgYml0IGZsaXBwaW5nICovCiAgICAgIGludAogICAgICAgIGksajsgICAgICAgICAgICAgIC8qIGluZGV4ICovCiAgICAgIHVuc2lnbmVkIGxvbmcKICAgICAgICBzLCAgICAgICAgICAgICAgICAvKiBjYWxjdWxhdGVkIHN5bmRyb21lICovCiAgICAgICAgY3dzYXZlcjsgICAgICAgICAgLyogc2F2ZXMgaW5pdGlhbCB2YWx1ZSBvZiBjdyAqLwogICAgCiAgICAgIGN3c2F2ZXI9Y3c7ICAgICAgICAgLyogc2F2ZSAqLwogICAgICAqZXJycz0wOwogICAgICB3PTM7ICAgICAgICAgICAgICAgIC8qIGluaXRpYWwgc3luZHJvbWUgd2VpZ2h0IHRocmVzaG9sZCAqLwogICAgICBqPS0xOyAgICAgICAgICAgICAgIC8qIC0xID0gbm8gdHJpYWwgYml0IGZsaXBwaW5nIG9uIGZpcnN0IHBhc3MgKi8KICAgICAgbWFzaz0xOwogICAgICB3aGlsZSAoajwyMykgLyogZmxpcCBlYWNoIHRyaWFsIGJpdCAqLwogICAgICAgIHsKICAgICAgICAgIGlmIChqICE9IC0xKSAvKiB0b2dnbGUgYSB0cmlhbCBiaXQgKi8KICAgICAgICAgICAgewogICAgICAgICAgICAgIGlmIChqPjApIC8qIHJlc3RvcmUgbGFzdCB0cmlhbCBiaXQgKi8KICAgICAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgICAgY3c9Y3dzYXZlciBeIG1hc2s7CiAgICAgICAgICAgICAgICAgIG1hc2srPW1hc2s7IC8qIHBvaW50IHRvIG5leHQgYml0ICovCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICAgY3c9Y3dzYXZlciBeIG1hc2s7IC8qIGZsaXAgbmV4dCB0cmlhbCBiaXQgKi8KICAgICAgICAgICAgICB3PTI7IC8qIGxvd2VyIHRoZSB0aHJlc2hvbGQgd2hpbGUgYml0IGRpZGRsaW5nICovCiAgICAgICAgICAgIH0KICAgIAogICAgICAgICAgcz1zeW5kcm9tZShjdyk7IC8qIGxvb2sgZm9yIGVycm9ycyAqLwogICAgICAgICAgaWYgKHMpIC8qIGVycm9ycyBleGlzdCAqLwogICAgICAgICAgICB7CiAgICAgICAgICAgICAgZm9yIChpPTA7IGk8MjM7IGkrKykgLyogY2hlY2sgc3luZHJvbWUgb2YgZWFjaCBjeWNsaWMgc2hpZnQgKi8KICAgICAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgICAgaWYgKCgqZXJycz13ZWlnaHQocykpIDw9IHcpIC8qIHN5bmRyb21lIG1hdGNoZXMgZXJyb3IgcGF0dGVybiAqLwogICAgICAgICAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgICAgICAgIGN3PWN3IF4gczsgICAgICAgICAgICAgIC8qIHJlbW92ZSBlcnJvcnMgKi8KICAgICAgICAgICAgICAgICAgICAgIGN3PXJvdGF0ZV9yaWdodChjdyxpKTsgIC8qIHVucm90YXRlIGRhdGEgKi8KICAgICAgICAgICAgICAgICAgICAgIHJldHVybihzPWN3KTsKICAgICAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICAgICAgIGVsc2UKICAgICAgICAgICAgICAgICAgICB7CiAgICAgICAgICAgICAgICAgICAgICBjdz1yb3RhdGVfbGVmdChjdywxKTsgICAvKiByb3RhdGUgdG8gbmV4dCBwYXR0ZXJuICovCiAgICAgICAgICAgICAgICAgICAgICBzPXN5bmRyb21lKGN3KTsgICAgICAgICAvKiBjYWxjIG5ldyBzeW5kcm9tZSAqLwogICAgICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgICBqKys7IC8qIHRvZ2dsZSBuZXh0IHRyaWFsIGJpdCAqLwogICAgICAgICAgICB9CiAgICAgICAgICBlbHNlCiAgICAgICAgICAgIHJldHVybihjdyk7IC8qIHJldHVybiBjb3JyZWN0ZWQgY29kZXdvcmQgKi8KICAgICAgICB9CiAgICAKICAgICAgcmV0dXJuKGN3c2F2ZXIpOyAvKiByZXR1cm4gb3JpZ2luYWwgaWYgbm8gY29ycmVjdGlvbnMgKi8KICAgIH0gLyogY29ycmVjdCAqLwogICAgCiAgICAvKiA9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT0gKi8KICAgIAogICAgaW50IGRlY29kZShpbnQgY29ycmVjdF9tb2RlLCBpbnQgKmVycnMsIHVuc2lnbmVkIGxvbmcgKmN3KQogICAgLyogVGhpcyBmdW5jdGlvbiBkZWNvZGVzIGNvZGV3b3JkICpjdyBpbiBvbmUgb2YgdHdvIG1vZGVzLiBJZiBjb3JyZWN0X21vZGUKICAgICAgaXMgbm9uemVybywgZXJyb3IgY29ycmVjdGlvbiBpcyBhdHRlbXB0ZWQsIHdpdGggKmVycnMgc2V0IHRvIHRoZSBudW1iZXIgb2YKICAgICAgYml0cyBjb3JyZWN0ZWQsIGFuZCByZXR1cm5pbmcgMCBpZiBubyBlcnJvcnMgZXhpc3QsIG9yIDEgaWYgcGFyaXR5IGVycm9ycwogICAgICBleGlzdC4gSWYgY29ycmVjdF9tb2RlIGlzIHplcm8sIGVycm9yIGRldGVjdGlvbiBpcyBwZXJmb3JtZWQgb24gKmN3LAogICAgICByZXR1cm5pbmcgMCBpZiBubyBlcnJvcnMgZXhpc3QsIDEgaWYgYW4gb3ZlcmFsbCBwYXJpdHkgZXJyb3IgZXhpc3RzLCBhbmQKICAgICAgMiBpZiBhIGNvZGV3b3JkIGVycm9yIGV4aXN0cy4gKi8KICAgIHsKICAgICAgdW5zaWduZWQgbG9uZyBwYXJpdHlfYml0OwogICAgCiAgICAgIGlmIChjb3JyZWN0X21vZGUpICAgICAgICAgICAgICAgLyogY29ycmVjdCBlcnJvcnMgKi8KICAgICAgICB7CiAgICAgICAgICBwYXJpdHlfYml0PSpjdyAmIDB4ODAwMDAwbDsgLyogc2F2ZSBwYXJpdHkgYml0ICovCiAgICAgICAgICAqY3cmPX4weDgwMDAwMGw7ICAgICAgICAgICAgLyogcmVtb3ZlIHBhcml0eSBiaXQgZm9yIGNvcnJlY3Rpb24gKi8KICAgIAogICAgICAgICAgKmN3PWNvcnJlY3QoKmN3LCBlcnJzKTsgICAgIC8qIGNvcnJlY3QgdXAgdG8gdGhyZWUgYml0cyAqLwogICAgICAgICAgKmN3fD1wYXJpdHlfYml0OyAgICAgICAgICAgIC8qIHJlc3RvcmUgcGFyaXR5IGJpdCAqLwogICAgCiAgICAgICAgICAvKiBjaGVjayBmb3IgNCBiaXQgZXJyb3JzICovCiAgICAgICAgICBpZiAocGFyaXR5KCpjdykpICAgICAgICAgICAgLyogb2RkIHBhcml0eSBpcyBhbiBlcnJvciAqLwogICAgICAgICAgICByZXR1cm4oMSk7CiAgICAgICAgICByZXR1cm4oMCk7IC8qIG5vIGVycm9ycyAqLwogICAgICAgIH0KICAgICAgZWxzZSAvKiBkZXRlY3QgZXJyb3JzIG9ubHkgKi8KICAgICAgICB7CiAgICAgICAgICAqZXJycz0wOwogICAgICAgICAgaWYgKHBhcml0eSgqY3cpKSAvKiBvZGQgcGFyaXR5IGlzIGFuIGVycm9yICovCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAqZXJycz0xOwogICAgICAgICAgICAgIHJldHVybigxKTsKICAgICAgICAgICAgfQogICAgICAgICAgaWYgKHN5bmRyb21lKCpjdykpCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAqZXJycz0xOwogICAgICAgICAgICAgIHJldHVybigyKTsKICAgICAgICAgICAgfQogICAgICAgICAgZWxzZQogICAgICAgICAgICByZXR1cm4oMCk7IC8qIG5vIGVycm9ycyAqLwogICAgICAgIH0KICAgIH0gLyogZGVjb2RlICovCiAgICAKICAgIC8qID09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PSAqLwogICAgCiAgICB2b2lkIGdvbGF5X3Rlc3Qodm9pZCkKICAgIC8qIFRoaXMgZnVuY3Rpb24gdGVzdHMgdGhlIEdvbGF5IHJvdXRpbmVzIGZvciBkZXRlY3Rpb24gYW5kIGNvcnJlY3Rpb24KICAgICAgb2YgdmFyaW91cyBwYXR0ZXJucyBvZiBlcnJvcl9saW1pdCBiaXQgZXJyb3JzLiBUaGUgZXJyb3JfbWFzayBjeWNsZXMKICAgICAgb3ZlciBhbGwgcG9zc2libGUgdmFsdWVzLCBhbmQgZXJyb3JfbGltaXQgc2VsZWN0cyB0aGUgbWF4aW11bSBudW1iZXIKICAgICAgb2YgaW5kdWNlZCBlcnJvcnMuICovCiAgICB7CiAgICAgIHVuc2lnbmVkIGxvbmcKICAgICAgICBlcnJvcl9tYXNrLCAgICAgICAgIC8qIGJpdHdpc2UgbWFzayBmb3IgaW5kdWNpbmcgZXJyb3JzICovCiAgICAgICAgdHJhc2hlZF9jb2Rld29yZCwgICAvKiB0aGUgY29kZXdvcmQgZm9yIHRyaWFsIGNvcnJlY3Rpb24gKi8KICAgICAgICB2aXJnaW5fY29kZXdvcmQ7ICAgIC8qIHRoZSBvcmlnaW5hbCBjb2Rld29yZCB3aXRob3V0IGVycm9ycyAqLwogICAgICB1bnNpZ25lZCBjaGFyCiAgICAgICAgcGFzcz0xLCAgICAgICAgICAgICAvKiBhc3N1bWUgdGVzdCBwYXNzZXMgKi8KICAgICAgICBlcnJvcl9saW1pdD0zOyAgICAgIC8qIHNlbGVjdCBudW1iZXIgb2YgaW5kdWNlZCBiaXQgZXJyb3JzIGhlcmUgKi8KICAgICAgaW50CiAgICAgICAgZXJyb3JfY291bnQ7ICAgICAgICAvKiByZWNlaXZlcyBudW1iZXIgb2YgZXJyb3JzIGNvcnJlY3RlZCAqLwogICAgCiAgICAgIHZpcmdpbl9jb2Rld29yZD1nb2xheSgweDU1NSk7IC8qIG1ha2UgYSB0ZXN0IGNvZGV3b3JkICovCiAgICAgIGlmIChwYXJpdHkodmlyZ2luX2NvZGV3b3JkKSkKICAgICAgICB2aXJnaW5fY29kZXdvcmRePTB4ODAwMDAwbDsKICAgICAgZm9yIChlcnJvcl9tYXNrPTA7IGVycm9yX21hc2s8MHg4MDAwMDBsOyBlcnJvcl9tYXNrKyspCiAgICAgICAgewogICAgICAgICAgLyogZmlsdGVyIHRoZSBtYXNrIGZvciB0aGUgc2VsZWN0ZWQgbnVtYmVyIG9mIGJpdCBlcnJvcnMgKi8KICAgICAgICAgIGlmICh3ZWlnaHQoZXJyb3JfbWFzaykgPD0gZXJyb3JfbGltaXQpIC8qIHlvdSBjYW4gbWFrZSB0aGlzIGZhc3RlciEgKi8KICAgICAgICAgICAgewogICAgICAgICAgICAgIHRyYXNoZWRfY29kZXdvcmQ9dmlyZ2luX2NvZGV3b3JkIF4gZXJyb3JfbWFzazsgLyogaW5kdWNlIGJpdCBlcnJvcnMgKi8KICAgIAogICAgICAgICAgICAgIGRlY29kZSgxLCZlcnJvcl9jb3VudCwmdHJhc2hlZF9jb2Rld29yZCk7IC8qIHRyeSB0byBjb3JyZWN0IGJpdCBlcnJvcnMgKi8KICAgIAogICAgICAgICAgICAgIGlmICh0cmFzaGVkX2NvZGV3b3JkIF4gdmlyZ2luX2NvZGV3b3JkKQogICAgICAgICAgICAgICAgewogICAgICAgICAgICAgICAgICBwcmludGYoIlVuYWJsZSB0byBjb3JyZWN0ICVkIGVycm9ycyBpbmR1Y2VkIHdpdGggZXJyb3IgbWFzayA9IDB4JWxYXG4iLAogICAgICAgICAgICAgICAgICAgIHdlaWdodChlcnJvcl9tYXNrKSxlcnJvcl9tYXNrKTsKICAgICAgICAgICAgICAgICAgcGFzcz0wOwogICAgICAgICAgICAgICAgfQogICAgCiAgICAgICAgICAgICAgaWYgKGtiaGl0KCkpIC8qIGxvb2sgZm9yIHVzZXIgaW5wdXQgKi8KICAgICAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgICAgaWYgKGdldGNoKCkgPT0gMjcpIHJldHVybjsgLyogZXNjYXBlIGV4aXRzICovCiAgICAKICAgICAgICAgICAgICAgICAgLyogb3RoZXIga2V5IHByaW50cyBzdGF0dXMgKi8KICAgICAgICAgICAgICAgICAgcHJpbnRmKCJDdXJyZW50IHRlc3QgY291bnQgPSAlbGQgb2YgJWxkXG4iLGVycm9yX21hc2ssMHg4MDAwMDBsKTsKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgfQogICAgICAgIH0KICAgICAgcHJpbnRmKCJHb2xheSB0ZXN0ICVzIVxuIixwYXNzPyJQQVNTRUQiOiJGQUlMRUQiKTsKICAgIH0KICAgIAogICAgLyogPT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09ICovCiAgICAKICAgIHZvaWQgbWFpbihpbnQgYXJndW1lbnRfY291bnQsIGNoYXIgKmFyZ3VtZW50W10pCiAgICB7CiAgICAgIGludCBpLGo7CiAgICAgIHVuc2lnbmVkIGxvbmcgbCxnOwogICAgICBjb25zdCBjaGFyICplcnJtc2cgPSAiVXNhZ2U6IEcgREFUQSBFbmNvZGUvQ29ycmVjdC9WZXJpZnkvVGVzdFxuXG4iCiAgICAgICAgICAgICAgICAgIiAgd2hlcmUgREFUQSBpcyB0aGUgZGF0YSB0byBiZSBlbmNvZGVkLCBjb2Rld29yZCB0byBiZSBjb3JyZWN0ZWQsXG4iCiAgICAgICAgICAgICAgICAgIiAgb3IgY29kZXdvcmQgdG8gYmUgY2hlY2tlZCBmb3IgZXJyb3JzLiBEQVRBIGlzIGhleGFkZWNpbWFsLlxuXG4iCiAgICAgICAgICAgICAgICAgIkV4YW1wbGVzOlxuXG4iCiAgICAgICAgICAgICAgICAgIiAgRyA1NTUgRSAgICAgIGVuY29kZXMgaW5mb3JtYXRpb24gdmFsdWUgNTU1IGFuZCBwcmludHMgYSBjb2Rld29yZFxuIgogICAgICAgICAgICAgICAgICIgIEcgQUJDMTIzIEMgICBjb3JyZWN0cyBjb2Rld29yZCBBQkMxMjNcbiIKICAgICAgICAgICAgICAgICAiICBHIEFCQzEyMyBWICAgY2hlY2tzIGNvZGV3b3JkIEFCQzEyMyBmb3IgZXJyb3JzXG4iCiAgICAgICAgICAgICAgICAgIiAgRyBBQkMxMjMgVCAgIHRlc3RzIHJvdXRpbmVzLCBBQkMxMjMgaXMgYSBkdW1teSBwYXJhbWV0ZXJcblxuIjsKICAgIAogICAgICBpZiAoYXJndW1lbnRfY291bnQgIT0gMykKICAgICAgICB7CiAgICAgICAgICBwcmludGYoZXJybXNnKTsKICAgICAgICAgIGV4aXQoMCk7CiAgICAgICAgfQogICAgCiAgICAgIGlmIChzc2NhbmYoYXJndW1lbnRbMV0sIiVseCIsJmwpICE9IDEpCiAgICAgICAgewogICAgICAgICAgcHJpbnRmKGVycm1zZyk7CiAgICAgICAgICBleGl0KDApOwogICAgICAgIH0KICAgIAogICAgICBzd2l0Y2ggKHRvdXBwZXIoKmFyZ3VtZW50WzJdKSkKICAgICAgICB7CiAgICAgICAgICBjYXNlICdFJzogLyogZW5jb2RlICovCiAgICAgICAgICAgIGwmPTB4ZmZmOwogICAgICAgICAgICBsPWdvbGF5KGwpOwogICAgICAgICAgICBpZiAocGFyaXR5KGwpKSBsXj0weDgwMDAwMGw7CiAgICAgICAgICAgIHByaW50ZigiQ29kZXdvcmQgPSAlbFhcbiIsbCk7CiAgICAgICAgICAgIGJyZWFrOwogICAgCiAgICAgICAgICBjYXNlICdWJzogLyogdmVyaWZ5ICovCiAgICAgICAgICAgIGlmIChkZWNvZGUoMCwmaSwmbCkpCiAgICAgICAgICAgICAgcHJpbnRmKCJDb2Rld29yZCAlbFggaXMgbm90IGEgR29sYXkgY29kZXdvcmQuXG4iLGwpOwogICAgICAgICAgICBlbHNlCiAgICAgICAgICAgICAgcHJpbnRmKCJDb2Rld29yZCAlbFggaXMgYSBHb2xheSBjb2Rld29yZC5cbiIsbCk7CiAgICAgICAgICAgIGJyZWFrOwogICAgCiAgICAgICAgICBjYXNlICdDJzogLyogY29ycmVjdCAqLwogICAgICAgICAgICBnPWw7IC8qIHNhdmUgaW5pdGlhbCBjb2Rld29yZCAqLwogICAgICAgICAgICBqPWRlY29kZSgxLCZpLCZsKTsKICAgICAgICAgICAgaWYgKChqKSAmJiAoaSkpCiAgICAgICAgICAgICAgcHJpbnRmKCJDb2Rld29yZCAlbFggaGFkICVkIGJpdHMgY29ycmVjdGVkLFxuIgogICAgICAgICAgICAgICAgICAgICAicmVzdWx0aW5nIGluIGNvZGV3b3JkICVsWCB3aXRoIGEgcGFyaXR5IGVycm9yLlxuIixnLGksbCk7CiAgICAgICAgICAgIGVsc2UKICAgICAgICAgICAgICBpZiAoKGogPT0gMCkgJiYgKGkpKQogICAgICAgICAgICAgICAgcHJpbnRmKCJDb2Rld29yZCAlbFggaGFkICVkIGJpdHMgY29ycmVjdGVkLCByZXN1bHRpbmcgaW4gY29kZXdvcmQgJWxYLlxuIixnLGksbCk7CiAgICAgICAgICAgICAgZWxzZQogICAgICAgICAgICAgICAgaWYgKChqKSAmJiAoaSA9PSAwKSkKICAgICAgICAgICAgICAgICAgcHJpbnRmKCJDb2Rld29yZCAlbFggaGFzIGEgcGFyaXR5IGVycm9yLiBObyBiaXRzIHdlcmUgY29ycmVjdGVkLlxuIixnKTsKICAgICAgICAgICAgICAgIGVsc2UKICAgICAgICAgICAgICAgICAgaWYgKChqID09IDApICYmIChpID09IDApKQogICAgICAgICAgICAgICAgICAgIHByaW50ZigiQ29kZXdvcmQgJWxYIGRvZXMgbm90IHJlcXVpcmUgY29ycmVjdGlvbi5cbiIsZyk7CiAgICAgICAgICAgIGJyZWFrOwogICAgCiAgICAgICAgICBjYXNlICdUJzogLyogdGVzdCAqLwogICAgICAgICAgICBwcmludGYoIlByZXNzIFNQQUNFIGZvciBzdGF0dXMsIEVTQyB0byBleGl0IHRlc3QuLi5cbiIpOwogICAgICAgICAgICBnb2xheV90ZXN0KCk7CiAgICAgICAgICAgIGJyZWFrOwogICAgCiAgICAgICAgICBkZWZhdWx0OgogICAgICAgICAgICBwcmludGYoZXJybXNnKTsKICAgICAgICAgICAgZXhpdCgwKTsKICAgICAgICB9CiAgICB9CiAgICAKICAgIC8qIGVuZCBvZiBHLkMgKi8KCgppbXBvcnQgamF2YS51dGlsLio7CmltcG9ydCBqYXZhLmxhbmcuKjsKaW1wb3J0IGphdmEuaW8uKjsKCi8qIE5hbWUgb2YgdGhlIGNsYXNzIGhhcyB0byBiZSAiTWFpbiIgb25seSBpZiB0aGUgY2xhc3MgaXMgcHVibGljLiAqLwpjbGFzcyBJZGVvbmUKewoJcHVibGljIHN0YXRpYyB2b2lkIG1haW4gKFN0cmluZ1tdIGFyZ3MpIHRocm93cyBqYXZhLmxhbmcuRXhjZXB0aW9uCgl7CgkJLy8geW91ciBjb2RlIGdvZXMgaGVyZQoJfQp9

I2luY2x1ZGUgJmx0O0xpcXVpZENyeXN0YWwuaCZndDsKCi8vIFBpbiBhc3NpZ25tZW50cwpjb25zdCBpbnQgZmxvd1NlbnNvclBpbiA9IDc7IC8vIENoYW5nZSBwaW4gY29ubmVjdGVkIHRvIHRoZSBmbG93IHNlbnNvciB0byBwaW4gNwpMaXF1aWRDcnlzdGFsIGxjZCgxMiwgMTEsIDUsIDQsIDMsIDIpOyAvLyBLZWVwIExDRCBwaW5zIGFzIHRoZXkgYXJlCgovLyBWYXJpYWJsZXMKdm9sYXRpbGUgaW50IGZsb3dDb3VudCA9IDA7IC8vIE51bWJlciBvZiBwdWxzZXMgY291bnRlZApmbG9hdCBmbG93UmF0ZSA9IDAuMDsgLy8gRmxvdyByYXRlIGluIGxpdGVycy9taW51dGUKZmxvYXQgdG90YWxWb2x1bWUgPSAwLjA7IC8vIFRvdGFsIHZvbHVtZSBpbiBsaXRlcnMKZmxvYXQgY29zdFBlckxpdGVyID0gMS41OyAvLyBDaGFuZ2UgdG8geW91ciBmdWVsIHByaWNlCmZsb2F0IHRvdGFsQ29zdCA9IDAuMDsgLy8gVG90YWwgY29zdAp1bnNpZ25lZCBsb25nIGxhc3RUaW1lID0gMDsgLy8gTGFzdCB0aW1lIHRoZSBmbG93IHJhdGUgd2FzIGNhbGN1bGF0ZWQKCnZvaWQgc2V0dXAoKSB7CiAgU2VyaWFsLmJlZ2luKDk2MDApOyAvLyBTdGFydCBzZXJpYWwgY29tbXVuaWNhdGlvbgogIGxjZC5iZWdpbigxNiwgMik7IC8vIEluaXRpYWxpemUgTENECiAgcGluTW9kZShmbG93U2Vuc29yUGluLCBJTlBVVCk7IC8vIFNldCBmbG93IHNlbnNvciBwaW4gYXMgaW5wdXQgIAogIGF0dGFjaEludGVycnVwdChkaWdpdGFsUGluVG9JbnRlcnJ1cHQoZmxvd1NlbnNvclBpbiksIGNvdW50RmxvdywgUklTSU5HKTsgLy8gQ291bnQgcHVsc2VzCn0KCnZvaWQgbG9vcCgpIHsKICAvLyBDYWxjdWxhdGUgZmxvdyByYXRlIGV2ZXJ5IHNlY29uZAogIGlmIChtaWxsaXMoKSAtIGxhc3RUaW1lICZndDs9IDEwMDApIHsKICAgIC8vIENhbGN1bGF0ZSBmbG93IHJhdGUgKGFkanVzdCBkaXZpc29yIGJhc2VkIG9uIHlvdXIgc2Vuc29yJ3Mgc3BlY3MpCiAgICBmbG93UmF0ZSA9IChmbG93Q291bnQgLyA3LjUpOyAvLyBFeGFtcGxlIGRpdmlzb3IgZm9yIFlGLVMyMDEKICAgIHRvdGFsVm9sdW1lICs9IChmbG93UmF0ZSAvIDYwKTsgLy8gQ29udmVydCBmbG93IHJhdGUgdG8gbGl0ZXJzIGFuZCBhZGQgdG8gdG90YWwKICAgIHRvdGFsQ29zdCA9IHRvdGFsVm9sdW1lICogY29zdFBlckxpdGVyOyAvLyBDYWxjdWxhdGUgdG90YWwgY29zdCBiYXNlZCBvbiBjb3N0IHBlciBsaXRlcgoKICAgIC8vIE91dHB1dCB0byBTZXJpYWwgTW9uaXRvcgogICAgU2VyaWFsLnByaW50KCZxdW90O0Zsb3cgUmF0ZTogJnF1b3Q7KTsKICAgIFNlcmlhbC5wcmludChmbG93UmF0ZSk7CiAgICBTZXJpYWwucHJpbnRsbigmcXVvdDsgTC9taW4mcXVvdDspOwogICAgU2VyaWFsLnByaW50KCZxdW90O1RvdGFsIFZvbHVtZTogJnF1b3Q7KTsKICAgIFNlcmlhbC5wcmludCh0b3RhbFZvbHVtZSk7CiAgICBTZXJpYWwucHJpbnRsbigmcXVvdDsgTCZxdW90Oyk7CiAgICBTZXJpYWwucHJpbnQoJnF1b3Q7VG90YWwgQ29zdDogJCZxdW90Oyk7CiAgICBTZXJpYWwucHJpbnRsbih0b3RhbENvc3QsIDIpOyAvLyBQcmludCBjb3N0IHdpdGggMiBkZWNpbWFsIHBsYWNlcwoKICAgIC8vIERpc3BsYXkgdmFsdWVzIG9uIExDRAogICAgbGNkLmNsZWFyKCk7CiAgICBsY2Quc2V0Q3Vyc29yKDAsIDApOwogICAgbGNkLnByaW50KCZxdW90O1ZvbDogJnF1b3Q7KTsKICAgIGxjZC5wcmludCh0b3RhbFZvbHVtZSwgMik7IC8vIFRvdGFsIHZvbHVtZQogICAgbGNkLnByaW50KCZxdW90OyBMJnF1b3Q7KTsKICAgIAogICAgbGNkLnNldEN1cnNvcigwLCAxKTsKICAgIGxjZC5wcmludCgmcXVvdDtDb3N0OiAkJnF1b3Q7KTsKICAgIGxjZC5wcmludCh0b3RhbENvc3QsIDIpOyAvLyBUb3RhbCBjb3N0CgogICAgLy8gUmVzZXQgZmxvdyBjb3VudAogICAgZmxvd0NvdW50ID0gMDsKICAgIGxhc3RUaW1lID0gbWlsbGlzKCk7IC8vIFVwZGF0ZSBsYXN0IHRpbWUKICB9Cn0KCnZvaWQgY291bnRGbG93KCkgewogIGZsb3dDb3VudCsrOyAvLyBJbmNyZW1lbnQgcHVsc2UgY291bnQKfQ==

#include <LiquidCrystal.h>

// Pin assignments
const int flowSensorPin = 7; // Change pin connected to the flow sensor to pin 7
LiquidCrystal lcd(12, 11, 5, 4, 3, 2); // Keep LCD pins as they are

// Variables
volatile int flowCount = 0; // Number of pulses counted
float flowRate = 0.0; // Flow rate in liters/minute
float totalVolume = 0.0; // Total volume in liters
float costPerLiter = 1.5; // Change to your fuel price
float totalCost = 0.0; // Total cost
unsigned long lastTime = 0; // Last time the flow rate was calculated

void setup() {
  Serial.begin(9600); // Start serial communication
  lcd.begin(16, 2); // Initialize LCD
  pinMode(flowSensorPin, INPUT); // Set flow sensor pin as input  
  attachInterrupt(digitalPinToInterrupt(flowSensorPin), countFlow, RISING); // Count pulses
}

void loop() {
  // Calculate flow rate every second
  if (millis() - lastTime >= 1000) {
    // Calculate flow rate (adjust divisor based on your sensor's specs)
    flowRate = (flowCount / 7.5); // Example divisor for YF-S201
    totalVolume += (flowRate / 60); // Convert flow rate to liters and add to total
    totalCost = totalVolume * costPerLiter; // Calculate total cost based on cost per liter

    // Output to Serial Monitor
    Serial.print("Flow Rate: ");
    Serial.print(flowRate);
    Serial.println(" L/min");
    Serial.print("Total Volume: ");
    Serial.print(totalVolume);
    Serial.println(" L");
    Serial.print("Total Cost: $");
    Serial.println(totalCost, 2); // Print cost with 2 decimal places

    // Display values on LCD
    lcd.clear();
    lcd.setCursor(0, 0);
    lcd.print("Vol: ");
    lcd.print(totalVolume, 2); // Total volume
    lcd.print(" L");
    
    lcd.setCursor(0, 1);
    lcd.print("Cost: $");
    lcd.print(totalCost, 2); // Total cost

    // Reset flow count
    flowCount = 0;
    lastTime = millis(); // Update last time
  }
}

void countFlow() {
  flowCount++; // Increment pulse count
}