Get Number of Times of Output Differential = 1 for each Output Indices of Each LFSR Byte.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Defining s0, s1, s2, s3, s4, s5, s6, s7, s8, s9, s10, s11, s12, s13, s14, s15, s16, s17, s18, s19, s20, s21, s22, s23, s24, s25, s26, s27, s28, s29, s30, s31, s32, s33, s34, s35, s36, s37, s38, s39, s40, s41, s42, s43, s44, s45, s46, s47, s48, s49, s50, s51, s52, s53, s54, s55, s56, s57, s58, s59, s60, s61, s62, s63, s64, s65, s66, s67, s68, s69, s70, s71, s72, s73, s74, s75, s76, s77, s78, s79, s80, s81, s82, s83, s84, s85, s86, s87, s88, s89, s90, s91, s92, s93, s94, s95, s96, s97, s98, s99, s100, s101, s102, s103, s104, s105, s106, s107, s108, s109, s110, s111, s112, s113, s114, s115, s116, s117, s118, s119, s120, s121, s122, s123, s124, s125, s126, s127, b0, b1, b2, b3, b4, b5, b6, b7, b8, b9, b10, b11, b12, b13, b14, b15, b16, b17, b18, b19, b20, b21, b22, b23, b24, b25, b26, b27, b28, b29, b30, b31, b32, b33, b34, b35, b36, b37, b38, b39, b40, b41, b42, b43, b44, b45, b46, b47, b48, b49, b50, b51, b52, b53, b54, b55, b56, b57, b58, b59, b60, b61, b62, b63, b64, b65, b66, b67, b68, b69, b70, b71, b72, b73, b74, b75, b76, b77, b78, b79, b80, b81, b82, b83, b84, b85, b86, b87, b88, b89, b90, b91, b92, b93, b94, b95, b96, b97, b98, b99, b100, b101, b102, b103, b104, b105, b106, b107, b108, b109, b110, b111, b112, b113, b114, b115, b116, b117, b118, b119, b120, b121, b122, b123, b124, b125, b126, b127, k0, k1, k2, k3, k4, k5, k6, k7, k8, k9, k10, k11, k12, k13, k14, k15, k16, k17, k18, k19, k20, k21, k22, k23, k24, k25, k26, k27, k28, k29, k30, k31, k32, k33, k34, k35, k36, k37, k38, k39, k40, k41, k42, k43, k44, k45, k46, k47, k48, k49, k50, k51, k52, k53, k54, k55, k56, k57, k58, k59, k60, k61, k62, k63, k64, k65, k66, k67, k68, k69, k70, k71, k72, k73, k74, k75, k76, k77, k78, k79, k80, k81, k82, k83, k84, k85, k86, k87, k88, k89, k90, k91, k92, k93, k94, k95, k96, k97, k98, k99, k100, k101, k102, k103, k104, k105, k106, k107, k108, k109, k110, k111, k112, k113, k114, k115, k116, k117, k118, k119, k120, k121, k122, k123, k124, k125, k126, k127\n",
      "Defining s0, s1, s2, s3, s4, s5, s6, s7, s8, s9, s10, s11, s12, s13, s14, s15, s16, s17, s18, s19, s20, s21, s22, s23, s24, s25, s26, s27, s28, s29, s30, s31, s32, s33, s34, s35, s36, s37, s38, s39, s40, s41, s42, s43, s44, s45, s46, s47, s48, s49, s50, s51, s52, s53, s54, s55, s56, s57, s58, s59, s60, s61, s62, s63, s64, s65, s66, s67, s68, s69, s70, s71, s72, s73, s74, s75, s76, s77, s78, s79, s80, s81, s82, s83, s84, s85, s86, s87, s88, s89, s90, s91, s92, s93, s94, s95, s96, s97, s98, s99, s100, s101, s102, s103, s104, s105, s106, s107, s108, s109, s110, s111, s112, s113, s114, s115, s116, s117, s118, s119, s120, s121, s122, s123, s124, s125, s126, s127, b0, b1, b2, b3, b4, b5, b6, b7, b8, b9, b10, b11, b12, b13, b14, b15, b16, b17, b18, b19, b20, b21, b22, b23, b24, b25, b26, b27, b28, b29, b30, b31, b32, b33, b34, b35, b36, b37, b38, b39, b40, b41, b42, b43, b44, b45, b46, b47, b48, b49, b50, b51, b52, b53, b54, b55, b56, b57, b58, b59, b60, b61, b62, b63, b64, b65, b66, b67, b68, b69, b70, b71, b72, b73, b74, b75, b76, b77, b78, b79, b80, b81, b82, b83, b84, b85, b86, b87, b88, b89, b90, b91, b92, b93, b94, b95, b96, b97, b98, b99, b100, b101, b102, b103, b104, b105, b106, b107, b108, b109, b110, b111, b112, b113, b114, b115, b116, b117, b118, b119, b120, b121, b122, b123, b124, b125, b126, b127\n",
      "Byte 0\n",
      "0|0|0|100|0|0|0|49|12|24|43|32|20|68\n",
      "0|0|0|0|0|0|0|100|18|0|0|0|18|20\n",
      "100|0|0|0|100|0|0|0|8|27|25|52|25|20\n",
      "26|0|0|0|0|0|0|0|12|100|0|100|12|14\n",
      "41|100|0|0|0|100|0|0|100|7|28|61|100|25\n",
      "0|0|0|0|0|0|0|0|53|13|100|0|57|15\n",
      "0|0|100|48|0|0|100|0|6|3|10|20|27|100\n",
      "0|0|0|0|0|0|0|0|55|13|49|0|0|54\n",
      "Byte 1\n",
      "0|0|46|0|0|0|100|0|27|0|90|51|59|36|35|51\n",
      "100|0|0|0|47|0|0|0|0|0|46|0|0|0|12|77\n",
      "33|0|0|100|0|0|0|100|0|51|8|20|51|53|51|37\n",
      "0|100|0|14|0|58|0|0|0|100|14|100|0|100|14|14\n",
      "0|21|0|58|100|0|0|0|100|35|100|87|30|48|100|46\n",
      "0|0|100|0|0|0|0|0|0|0|52|50|100|0|54|14\n",
      "0|0|28|0|0|100|45|0|0|0|2|4|95|56|33|100\n",
      "0|0|0|0|0|0|0|0|0|0|46|11|46|0|0|58\n",
      "Byte 2\n",
      "0|0|52|0|0|0|100|0|29|0|83|49|40|40|40|43\n",
      "100|0|0|0|53|0|0|0|0|0|45|0|0|0|7|74\n",
      "25|0|0|100|0|0|0|100|0|58|6|31|46|43|48|37\n",
      "0|100|0|31|0|41|0|0|0|100|11|100|0|100|11|16\n",
      "0|25|0|52|100|0|0|0|100|24|100|87|32|52|100|52\n",
      "0|0|100|0|0|0|0|0|0|0|46|47|100|0|50|13\n",
      "0|0|24|0|0|100|55|0|0|0|6|6|90|53|38|100\n",
      "0|0|0|0|0|0|0|0|0|0|44|20|38|0|0|52\n",
      "Byte 3\n",
      "0|0|47|0|0|0|100|0|27|0|92|58|57|33|31|53\n",
      "100|0|0|0|49|0|0|0|0|0|52|0|0|0|11|75\n",
      "21|0|0|100|0|0|0|100|0|41|11|26|51|48|51|43\n",
      "0|100|0|17|0|46|0|0|0|100|13|100|0|100|13|19\n",
      "0|22|0|58|100|0|0|0|100|32|100|85|28|55|100|56\n",
      "0|0|100|0|0|0|0|0|0|0|46|53|100|0|53|10\n",
      "0|0|20|0|0|100|56|0|0|0|4|6|90|48|30|100\n",
      "0|0|0|0|0|0|0|0|0|0|49|12|50|0|0|61\n",
      "Byte 4\n",
      "0|0|0|49|0|0|0|0|100|0|25|46|0|89|51|54|38|39|35|47\n",
      "0|100|0|0|0|51|0|0|0|0|0|100|0|48|0|0|0|11|9|28\n",
      "0|26|0|0|100|0|0|47|0|100|0|26|43|5|25|56|56|100|54|45\n",
      "0|0|100|0|29|0|55|0|0|0|0|0|100|9|100|0|100|53|9|50\n",
      "0|0|31|0|46|100|0|0|0|0|100|0|27|100|85|16|55|44|100|8\n",
      "0|0|0|100|0|0|0|0|0|0|0|0|0|51|54|100|0|0|49|28\n",
      "100|0|0|25|0|0|100|0|49|0|0|0|0|55|10|57|45|42|46|48\n",
      "19|0|0|0|53|0|0|100|53|0|0|25|0|41|13|61|0|50|0|100\n",
      "Byte 5\n",
      "100|0|0|0|100|0|49|0|44|0|44|0|42|45|48|59|51|42|56|40|46\n",
      "28|0|0|0|70|51|0|0|50|0|45|0|100|49|0|51|47|100|54|58|24\n",
      "0|100|0|0|24|100|0|100|0|0|0|46|0|8|60|39|45|64|47|44|52\n",
      "0|17|0|0|0|73|40|25|0|100|0|48|53|10|100|0|54|38|100|57|46\n",
      "0|0|100|0|0|34|100|55|100|51|0|0|54|100|49|48|51|47|66|50|52\n",
      "45|0|20|0|0|0|76|0|0|49|100|0|77|48|39|100|0|42|35|100|55\n",
      "0|0|0|100|0|0|24|0|0|0|49|0|24|41|8|42|42|0|49|64|61\n",
      "0|44|0|23|0|0|0|48|0|0|48|100|56|53|22|52|100|70|53|31|100\n",
      "Byte 6\n",
      "100|0|0|0|100|0|41|0|62|0|50|0|51|56|53|46|49|43|53|56|42\n",
      "17|0|0|0|67|52|0|0|51|0|52|0|100|46|0|48|51|100|49|61|32\n",
      "0|100|0|0|24|100|0|100|0|0|0|52|0|9|53|47|45|63|58|44|53\n",
      "0|28|0|0|0|87|57|22|0|100|0|48|62|14|100|0|49|40|100|47|48\n",
      "0|0|100|0|0|18|100|51|100|43|0|0|54|100|58|54|43|41|68|42|51\n",
      "44|0|24|0|0|0|74|0|0|47|100|0|82|53|50|100|0|51|30|100|53\n",
      "0|0|0|100|0|0|26|0|0|0|49|0|29|52|12|52|55|0|48|60|54\n",
      "0|55|0|18|0|0|0|45|0|0|50|100|47|48|14|46|100|72|49|29|100\n",
      "Byte 7\n",
      "100|0|0|0|100|0|57|0|44|0|48|0|50|41|47|48|55|52|47|44|47\n",
      "25|0|0|0|80|42|0|0|53|0|49|0|100|49|0|57|44|100|52|51|23\n",
      "0|100|0|0|19|100|0|100|0|0|0|47|0|7|49|47|51|59|47|50|54\n",
      "0|20|0|0|0|78|44|19|0|100|0|51|45|17|100|0|49|34|100|48|52\n",
      "0|0|100|0|0|26|100|54|100|49|0|0|52|100|47|45|46|51|69|53|44\n",
      "58|0|25|0|0|0|81|0|0|45|100|0|70|50|37|100|0|56|26|100|48\n",
      "0|0|0|100|0|0|24|0|0|0|50|0|29|58|6|50|46|0|51|64|60\n",
      "0|50|0|26|0|0|0|48|0|0|45|100|54|45|9|57|100|71|50|27|100\n",
      "Byte 8\n",
      "0|100|0|0|0|0|44|0|42|0|46|0|50|50|49|56|56|42|53|49\n",
      "0|30|0|0|0|52|0|0|60|0|53|0|75|100|46|0|42|41|55|46\n",
      "0|0|100|0|0|100|0|100|0|0|0|42|22|0|13|46|51|0|47|48\n",
      "0|0|24|0|0|70|45|20|0|100|0|43|50|59|10|100|0|69|100|48\n",
      "0|0|0|100|0|23|100|43|100|52|0|0|0|53|100|51|46|32|74|50\n",
      "0|46|0|25|0|0|74|0|0|50|100|0|0|79|45|47|100|0|33|100\n",
      "100|0|0|0|0|0|18|0|51|0|24|57|100|72|5|55|67|54|55|49\n",
      "18|0|50|0|100|0|44|0|0|0|57|100|0|50|50|12|49|100|12|57\n",
      "Byte 9\n",
      "100|0|0|0|0|0|0|0|100|0|17|0|42|51|34|0|37|0|100|0|79|44|72|50|50|49|37|52\n",
      "22|0|0|0|100|0|53|0|29|0|51|0|0|0|46|0|0|49|0|46|58|0|48|0|100|0|58|48\n",
      "0|100|0|0|24|0|0|0|0|100|0|0|0|100|0|54|0|49|0|100|0|56|45|67|56|44|56|50\n",
      "0|29|0|0|0|100|0|51|0|22|0|100|0|36|0|100|0|0|0|0|49|100|12|46|24|100|12|48\n",
      "0|0|100|0|0|24|0|0|0|0|100|25|0|50|100|24|53|0|0|0|100|73|100|44|57|53|100|42\n",
      "0|0|24|0|0|0|100|0|0|0|31|56|100|0|58|52|100|0|52|0|0|47|45|8|0|19|48|8\n",
      "0|0|0|100|0|0|26|0|0|0|0|0|27|0|54|0|21|55|39|0|0|0|5|100|45|52|49|100\n",
      "41|0|0|20|0|0|0|100|54|0|0|0|49|0|0|0|49|100|0|48|0|59|46|56|0|0|0|57\n",
      "Byte 10\n",
      "0|0|0|0|100|0|0|0|45|0|0|0|0|0|0|0|0|18|0|55|50|0|47|0|0|0|70|57|49|1|49|54|61|44|52|41|50\n",
      "100|0|52|0|24|0|51|0|100|0|0|40|0|0|0|100|0|0|0|49|0|0|100|47|44|51|28|49|49|41|8|100|0|47|53|9|41\n",
      "26|0|0|0|0|100|0|0|30|49|0|0|100|0|58|27|45|0|0|0|100|57|20|53|42|100|0|42|0|100|46|9|29|45|54|21|43\n",
      "0|100|0|50|0|27|0|47|0|100|0|0|51|0|0|48|100|0|100|0|43|56|0|100|50|0|48|100|61|49|59|52|100|0|100|11|9\n",
      "0|21|0|0|0|0|100|0|0|23|46|0|0|100|0|54|26|100|56|0|50|45|41|25|49|0|100|48|50|52|32|42|12|19|53|100|20\n",
      "0|0|100|0|0|0|29|0|0|0|100|0|46|54|0|0|60|32|0|100|52|0|0|0|100|52|0|24|100|56|45|13|53|100|54|45|13\n",
      "0|0|30|0|0|0|0|100|0|0|18|55|0|0|100|0|38|0|0|60|0|100|0|50|29|36|0|0|50|0|100|52|46|9|48|28|100\n",
      "0|0|0|100|58|0|0|30|50|0|0|100|0|49|59|0|0|0|0|0|0|43|40|0|0|37|50|52|25|53|45|57|11|45|0|54|45\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Byte 11\n",
      "0|0|0|0|0|100|0|0|0|52|0|0|0|100|0|47|0|42|0|64|0|54|50|45|0|51|54|0|0|52|0|47|24|46|47|20|45|47|27|44|55\n",
      "0|100|0|50|0|25|0|52|0|100|0|0|53|42|0|0|0|100|0|0|0|45|0|49|0|100|54|47|51|30|48|52|46|47|15|100|0|54|43|14|47\n",
      "0|34|0|0|0|0|100|0|0|18|42|0|0|0|100|0|47|29|49|0|0|0|100|52|48|24|48|50|100|0|47|0|100|29|47|11|26|48|50|27|49\n",
      "0|0|100|0|51|0|27|0|53|0|100|0|0|49|46|0|0|53|100|0|100|0|42|0|49|0|100|47|0|57|100|50|52|50|54|44|100|0|100|16|13\n",
      "0|0|23|0|0|0|0|100|0|0|26|46|0|0|0|100|0|52|19|100|48|0|48|100|49|48|25|48|0|100|44|41|56|100|30|47|11|23|48|100|26\n",
      "0|0|0|100|0|0|0|21|0|0|0|100|0|0|51|56|0|0|54|26|0|100|40|44|0|0|0|100|50|0|24|100|48|58|54|13|47|100|45|57|13\n",
      "100|0|0|27|0|0|0|0|100|0|0|25|53|0|0|0|100|0|49|0|0|57|0|57|100|0|63|23|45|0|0|48|0|48|100|45|51|10|50|22|100\n",
      "24|0|0|0|100|49|0|0|25|49|0|0|100|0|0|52|60|0|0|0|0|0|0|62|45|50|0|0|64|52|48|27|49|56|56|52|20|55|0|48|46\n",
      "Byte 12\n",
      "100|0|0|0|100|0|0|0|0|0|0|0|100|0|0|0|0|0|100|0|0|39|0|42|0|58|0|100|43|0|57|100|20|58|81\n",
      "33|0|0|0|28|47|0|0|100|0|48|0|24|0|50|49|0|0|23|54|0|100|0|0|62|51|40|23|21|41|9|72|58|33|39\n",
      "0|100|0|0|0|100|0|0|51|0|0|0|0|100|0|100|0|0|0|100|0|46|51|0|100|0|48|0|100|40|0|26|100|22|58\n",
      "0|28|0|0|0|31|52|0|0|100|0|49|0|26|0|20|48|0|0|23|0|49|100|38|75|46|57|46|27|25|54|44|73|56|29\n",
      "0|0|100|0|0|0|100|0|0|51|0|0|0|0|100|0|100|0|0|0|47|48|57|100|60|100|0|0|0|100|43|48|31|100|28\n",
      "0|0|25|0|0|0|19|58|0|0|100|0|0|0|29|45|23|55|52|0|51|0|54|24|0|77|55|16|56|19|4|58|49|79|52\n",
      "0|0|0|100|0|0|0|100|0|0|51|0|0|0|0|0|0|100|0|0|100|0|52|0|0|42|100|0|0|0|100|50|45|29|100\n",
      "0|0|0|27|0|0|0|21|0|0|0|100|47|0|0|0|56|22|53|55|28|50|0|0|44|0|78|53|25|44|45|61|50|50|77\n",
      "Byte 13\n",
      "100|0|0|0|100|0|0|0|0|0|0|0|100|0|0|0|0|0|100|0|0|55|0|50|0|58|0|100|50|0|46|100|24|49|71\n",
      "28|0|0|0|26|53|0|0|100|0|59|0|23|0|48|46|0|0|19|53|0|100|0|0|53|53|51|30|27|46|13|72|48|24|55\n",
      "0|100|0|0|0|100|0|0|55|0|0|0|0|100|0|100|0|0|0|100|0|54|53|0|100|0|52|0|100|54|0|20|100|21|53\n",
      "0|22|0|0|0|29|54|0|0|100|0|54|0|25|0|23|56|0|0|33|0|55|100|54|80|45|50|55|24|16|41|46|76|55|31\n",
      "0|0|100|0|0|0|100|0|0|50|0|0|0|0|100|0|100|0|0|0|50|61|49|100|45|100|0|0|0|100|52|50|23|100|27\n",
      "0|0|26|0|0|0|27|50|0|0|100|0|0|0|35|46|22|45|53|0|54|0|44|27|0|78|53|30|49|23|13|45|44|77|52\n",
      "0|0|0|100|0|0|0|100|0|0|45|0|0|0|0|0|0|100|0|0|100|0|40|0|0|53|100|0|0|0|100|53|49|22|100\n",
      "0|0|0|27|0|0|0|29|0|0|0|100|42|0|0|0|45|30|48|42|17|58|0|0|45|0|68|42|27|52|48|54|51|61|79\n",
      "Byte 14\n",
      "100|0|0|0|100|0|0|0|0|0|0|0|100|0|0|0|0|0|100|0|0|51|0|47|0|41|0|100|46|0|53|100|31|51|75\n",
      "28|0|0|0|28|53|0|0|100|0|51|0|28|0|50|51|0|0|26|48|0|100|0|0|50|48|46|22|23|44|13|79|50|29|50\n",
      "0|100|0|0|0|100|0|0|47|0|0|0|0|100|0|100|0|0|0|100|0|61|44|0|100|0|44|0|100|53|0|17|100|23|53\n",
      "0|27|0|0|0|29|47|0|0|100|0|50|0|29|0|29|48|0|0|23|0|54|100|52|75|52|58|47|24|29|42|53|76|56|27\n",
      "0|0|100|0|0|0|100|0|0|50|0|0|0|0|100|0|100|0|0|0|52|55|46|100|56|100|0|0|0|100|47|54|22|100|24\n",
      "0|0|20|0|0|0|27|45|0|0|100|0|0|0|23|53|26|46|46|0|55|0|48|24|0|77|61|21|45|22|14|46|59|80|51\n",
      "0|0|0|100|0|0|0|100|0|0|53|0|0|0|0|0|0|100|0|0|100|0|48|0|0|56|100|0|0|0|100|43|53|23|100\n",
      "0|0|0|31|0|0|0|31|0|0|0|100|50|0|0|0|59|19|59|49|27|44|0|0|50|0|78|47|31|50|47|49|48|42|84\n"
     ]
    }
   ],
   "source": [
    "# get number of time output differential = 1 grouped by single output index\n",
    "\n",
    "from sage.crypto.boolean_function import BooleanFunction\n",
    "from random import randint\n",
    "import random\n",
    "import numpy\n",
    "import csv\n",
    "import pandas as pd\n",
    "\n",
    "R.<s0, s1, s2, s3, s4, s5, s6, s7, s8, s9, s10, s11, s12, s13, s14, s15, s16, s17, s18, s19, s20, s21, s22, s23, s24, s25, s26, s27, s28, s29, s30, s31, s32, s33, s34, s35, s36, s37, s38, s39, s40, s41, s42, s43, s44, s45, s46, s47, s48, s49, s50, s51, s52, s53, s54, s55, s56, s57, s58, s59, s60, s61, s62, s63, s64, s65, s66, s67, s68, s69, s70, s71, s72, s73, s74, s75, s76, s77, s78, s79, s80, s81, s82, s83, s84, s85, s86, s87, s88, s89, s90, s91, s92, s93, s94, s95, s96, s97, s98, s99, s100, s101, s102, s103, s104, s105, s106, s107, s108, s109, s110, s111, s112, s113, s114, s115, s116, s117, s118, s119, s120, s121, s122, s123, s124, s125, s126, s127, b0, b1, b2, b3, b4, b5, b6, b7, b8, b9, b10, b11, b12, b13, b14, b15, b16, b17, b18, b19, b20, b21, b22, b23, b24, b25, b26, b27, b28, b29, b30, b31, b32, b33, b34, b35, b36, b37, b38, b39, b40, b41, b42, b43, b44, b45, b46, b47, b48, b49, b50, b51, b52, b53, b54, b55, b56, b57, b58, b59, b60, b61, b62, b63, b64, b65, b66, b67, b68, b69, b70, b71, b72, b73, b74, b75, b76, b77, b78, b79, b80, b81, b82, b83, b84, b85, b86, b87, b88, b89, b90, b91, b92, b93, b94, b95, b96, b97, b98, b99, b100, b101, b102, b103, b104, b105, b106, b107, b108, b109, b110, b111, b112, b113, b114, b115, b116, b117, b118, b119, b120, b121, b122, b123, b124, b125, b126, b127, k0, k1, k2, k3, k4, k5, k6, k7, k8, k9, k10, k11, k12, k13, k14, k15, k16, k17, k18, k19, k20, k21, k22, k23, k24, k25, k26, k27, k28, k29, k30, k31, k32, k33, k34, k35, k36, k37, k38, k39, k40, k41, k42, k43, k44, k45, k46, k47, k48, k49, k50, k51, k52, k53, k54, k55, k56, k57, k58, k59, k60, k61, k62, k63, k64, k65, k66, k67, k68, k69, k70, k71, k72, k73, k74, k75, k76, k77, k78, k79, k80, k81, k82, k83, k84, k85, k86, k87, k88, k89, k90, k91, k92, k93, k94, k95, k96, k97, k98, k99, k100, k101, k102, k103, k104, k105, k106, k107, k108, k109, k110, k111, k112, k113, k114, k115, k116, k117, k118, k119, k120, k121, k122, k123, k124, k125, k126, k127> = BooleanPolynomialRing()\n",
    "R.inject_variables()\n",
    "\n",
    "C.<s0, s1, s2, s3, s4, s5, s6, s7, s8, s9, s10, s11, s12, s13, s14, s15, s16, s17, s18, s19, s20, s21, s22, s23, s24, s25, s26, s27, s28, s29, s30, s31, s32, s33, s34, s35, s36, s37, s38, s39, s40, s41, s42, s43, s44, s45, s46, s47, s48, s49, s50, s51, s52, s53, s54, s55, s56, s57, s58, s59, s60, s61, s62, s63, s64, s65, s66, s67, s68, s69, s70, s71, s72, s73, s74, s75, s76, s77, s78, s79, s80, s81, s82, s83, s84, s85, s86, s87, s88, s89, s90, s91, s92, s93, s94, s95, s96, s97, s98, s99, s100, s101, s102, s103, s104, s105, s106, s107, s108, s109, s110, s111, s112, s113, s114, s115, s116, s117, s118, s119, s120, s121, s122, s123, s124, s125, s126, s127, b0, b1, b2, b3, b4, b5, b6, b7, b8, b9, b10, b11, b12, b13, b14, b15, b16, b17, b18, b19, b20, b21, b22, b23, b24, b25, b26, b27, b28, b29, b30, b31, b32, b33, b34, b35, b36, b37, b38, b39, b40, b41, b42, b43, b44, b45, b46, b47, b48, b49, b50, b51, b52, b53, b54, b55, b56, b57, b58, b59, b60, b61, b62, b63, b64, b65, b66, b67, b68, b69, b70, b71, b72, b73, b74, b75, b76, b77, b78, b79, b80, b81, b82, b83, b84, b85, b86, b87, b88, b89, b90, b91, b92, b93, b94, b95, b96, b97, b98, b99, b100, b101, b102, b103, b104, b105, b106, b107, b108, b109, b110, b111, b112, b113, b114, b115, b116, b117, b118, b119, b120, b121, b122, b123, b124, b125, b126, b127> = BooleanPolynomialRing()\n",
    "C.inject_variables()\n",
    "\n",
    "S_256 = [R.gens()[u] for u in range(128)] #Original State at S^[256]\n",
    "S = [R.gens()[u] for u in range(128)] #State at S^[256]\n",
    "K = [R.gens()[u] for u in range(256,384,1)]\n",
    "B = [R.gens()[u] for u in range(128,256,1)]\n",
    "\n",
    "def grain(State_S, State_B):\n",
    "    #update function of LFSR and NFSR\n",
    "    feedback_L = State_S[0] + State_S[7] + State_S[38] + State_S[70] + State_S[81] + State_S[96]\n",
    "    feedback_N = State_S[0] + State_B[0] + State_B[26] + State_B[56] + State_B[91] + State_B[96] + State_B[3] * State_B[67] + State_B[11] * State_B[13] + State_B[17] * State_B[18] + State_B[27] * State_B[59] + State_B[40] * State_B[48] + State_B[61] * State_B[65] + State_B[68] * State_B[84] + State_B[22] * State_B[24] * State_B[25] + State_B[70] * State_B[78] * State_B[82] + State_B[88] * State_B[92] * State_B[93] * State_B[95]\n",
    "\n",
    "    #feedback update\n",
    "    for j in range(128):\n",
    "        if j<=126:\n",
    "            State_S[j] = State_S[j+1]\n",
    "            State_B[j] = State_B[j+1]\n",
    "        else:\n",
    "            State_S[127] = feedback_L\n",
    "            State_B[127] = feedback_N\n",
    "        \n",
    "    #h function\n",
    "    h = State_B[12]*State_S[8] + State_S[13]*State_S[20] + State_B[95]*State_S[42] + State_S[60]*State_S[79] + State_B[12]*State_B[95]*State_S[94]\n",
    "    \n",
    "    #pre-output function\n",
    "    z = h + State_S[93] + State_B[2]+ State_B[15] + State_B[36] + State_B[45] + State_B[64] + State_B[73] + State_B[89]\n",
    "    \n",
    "    return z\n",
    "\n",
    "# byte = 1\n",
    "# val_list = [36, 38, 40, 42, 44, 46, 48, 50, 52, 62, 64, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96, 98, 100]\n",
    "# print ('|'.join(map(str, val_list)))\n",
    "\n",
    "# 2D array that holds 16 groups of unique output indices of each byte\n",
    "unique_output_indices = []\n",
    "\n",
    "with open(\"./csv/unique_output_indices.csv\", \"r\", newline='') as f:\n",
    "    reader = csv.reader(f, delimiter=',')\n",
    "    for l in reader:\n",
    "        int_list = [int(i) for i in l]\n",
    "        unique_output_indices.append(int_list[:])\n",
    "\n",
    "\n",
    "# total 16 bytes in 128 registers\n",
    "for byte_index, row in enumerate(unique_output_indices, start=0):\n",
    "# for row in unique_output_indices:\n",
    "    byte = byte_index\n",
    "    val_list = row  \n",
    "    \n",
    "    array_2d = []\n",
    "\n",
    "    for reg in range(8):\n",
    "        for output_index in range(len(val_list)): \n",
    "            faulty_state = (8 * byte) + reg\n",
    "            result = [0]*(len(val_list))\n",
    "\n",
    "            for i in range(100): \n",
    "                #Test for 100 random initial state\n",
    "                #initialize the internal states with random value (random Key and IV)\n",
    "\n",
    "                S = [R(random.randint(0,1)) for u in range(128)]\n",
    "                B = [R(random.randint(0,1)) for u in range(128,256,1)]\n",
    "                statef_S = [C() for u in range(128)]\n",
    "                statef_B = [C() for u in range(128,256,1)]\n",
    "                for l in range(128):\n",
    "                    statef_S[l] = S[l]\n",
    "                    statef_B[l] = B[l]\n",
    "\n",
    "                statef_S[faulty_state] = statef_S[faulty_state] + 1\n",
    "\n",
    "                #Generate the keystream for the desired number of rounds\n",
    "                j=0\n",
    "                for keystream in range(max(val_list)+1):\n",
    "                    zff = grain(S, B)\n",
    "                    zf = grain(statef_S, statef_B)\n",
    "                    d = zf + zff\n",
    "                    if keystream == val_list[j]:\n",
    "                        if d == 1:\n",
    "                            result[j] += 1\n",
    "                        if j < (len(val_list)-1):\n",
    "                            j+=1\n",
    "\n",
    "        array_2d.append(result)\n",
    "\n",
    "        \n",
    "    # # array printed separated with | \n",
    "    # # can be pasted in excel with delimeter\n",
    "    print(\"Byte %s\"  % (byte_index))\n",
    "    for k in range(len(array_2d)):\n",
    "        print ('|'.join(map(str, array_2d[k])))\n",
    "\n",
    "    # store number of times of each outout indices = 1 in 100 tests\n",
    "    csv_path = \"./csv/byte_%s_result.csv\" % (byte_index)\n",
    "    with open(csv_path, mode='w') as csv_file:\n",
    "        reg_writer = csv.writer(\n",
    "            csv_file, delimiter=',', quotechar='\"', quoting=csv.QUOTE_MINIMAL)\n",
    "\n",
    "        # write the header\n",
    "        # the header consists of the output indices\n",
    "        reg_writer.writerow(val_list) \n",
    "        for row in array_2d: \n",
    "            reg_writer.writerow(row)\n",
    "    \n",
    "print(\"End\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "SageMath 9.1",
   "language": "sage",
   "name": "sagemath"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}