Eq_list_diff_B=1 (Numerical experiment up to 200 rounds).ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "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\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",
      "Test 0\n",
      "Test 1\n",
      "Test 2\n",
      "Test 3\n",
      "Test 4\n",
      "Test 5\n",
      "Test 6\n",
      "Test 7\n",
      "Test 8\n",
      "Test 9\n",
      "Test 10\n",
      "Test 11\n",
      "Test 12\n",
      "Test 13\n",
      "Test 14\n",
      "Test 15\n",
      "Test 16\n",
      "Test 17\n",
      "Test 18\n",
      "Test 19\n",
      "Test 20\n",
      "Test 21\n",
      "Test 22\n",
      "Test 23\n",
      "Test 24\n",
      "Test 25\n",
      "Test 26\n",
      "Test 27\n",
      "Test 28\n",
      "Test 29\n",
      "Test 30\n",
      "Test 31\n",
      "Test 32\n",
      "Test 33\n",
      "Test 34\n",
      "Test 35\n",
      "Test 36\n",
      "Test 37\n",
      "Test 38\n",
      "Test 39\n",
      "Test 40\n",
      "Test 41\n",
      "Test 42\n",
      "Test 43\n",
      "Test 44\n",
      "Test 45\n",
      "Test 46\n",
      "Test 47\n",
      "Test 48\n",
      "Test 49\n",
      "Test 50\n",
      "Test 51\n",
      "Test 52\n",
      "Test 53\n",
      "Test 54\n",
      "Test 55\n",
      "Test 56\n",
      "Test 57\n",
      "Test 58\n",
      "Test 59\n",
      "Test 60\n",
      "Test 61\n",
      "Test 62\n",
      "Test 63\n",
      "Test 64\n",
      "Test 65\n",
      "Test 66\n",
      "Test 67\n",
      "Test 68\n",
      "Test 69\n",
      "Test 70\n",
      "Test 71\n",
      "Test 72\n",
      "Test 73\n",
      "Test 74\n",
      "Test 75\n",
      "Test 76\n",
      "Test 77\n",
      "Test 78\n",
      "Test 79\n",
      "Test 80\n",
      "Test 81\n",
      "Test 82\n",
      "Test 83\n",
      "Test 84\n",
      "Test 85\n",
      "Test 86\n",
      "Test 87\n",
      "Test 88\n",
      "Test 89\n",
      "Test 90\n",
      "Test 91\n",
      "Test 92\n",
      "Test 93\n",
      "Test 94\n",
      "Test 95\n",
      "Test 96\n",
      "Test 97\n",
      "Test 98\n",
      "Test 99\n",
      "Fault target 0\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 1\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 2\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 3\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 4\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 5\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 6\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 7\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 8\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 9\n",
      "Output Indices: 6\n",
      "Output Indices: 72\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 10\n",
      "Output Indices: 48\n",
      "Output Indices: 64\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 11\n",
      "Output Indices: 8\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 12\n",
      "Output Indices: 50\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 13\n",
      "Output Indices: 10\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 14\n",
      "Output Indices: 52\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 15\n",
      "Output Indices: 12\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 16\n",
      "Output Indices: 0\n",
      "Output Indices: 54\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 17\n",
      "Output Indices: 14\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 18\n",
      "Output Indices: 2\n",
      "Output Indices: 56\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 19\n",
      "Output Indices: 16\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 20\n",
      "Output Indices: 4\n",
      "Output Indices: 58\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 21\n",
      "Output Indices: 18\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 22\n",
      "Output Indices: 6\n",
      "Output Indices: 60\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 23\n",
      "Output Indices: 20\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 24\n",
      "Output Indices: 8\n",
      "Output Indices: 62\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 25\n",
      "Output Indices: 22\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 26\n",
      "Output Indices: 10\n",
      "Output Indices: 38\n",
      "Output Indices: 54\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 27\n",
      "Output Indices: 24\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 28\n",
      "Output Indices: 12\n",
      "Output Indices: 40\n",
      "Output Indices: 56\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 29\n",
      "Output Indices: 26\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 30\n",
      "Output Indices: 14\n",
      "Output Indices: 42\n",
      "Output Indices: 58\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 31\n",
      "Output Indices: 28\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 32\n",
      "Output Indices: 16\n",
      "Output Indices: 44\n",
      "Output Indices: 60\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 33\n",
      "Output Indices: 30\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 34\n",
      "Output Indices: 18\n",
      "Output Indices: 46\n",
      "Output Indices: 62\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 35\n",
      "Output Indices: 32\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 36\n",
      "Output Indices: 20\n",
      "Output Indices: 48\n",
      "Output Indices: 64\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 37\n",
      "Output Indices: 0\n",
      "Output Indices: 34\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 38\n",
      "Output Indices: 22\n",
      "Output Indices: 50\n",
      "Output Indices: 66\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 39\n",
      "Output Indices: 2\n",
      "Output Indices: 36\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 40\n",
      "Output Indices: 24\n",
      "Output Indices: 52\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 41\n",
      "Output Indices: 4\n",
      "Output Indices: 38\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 42\n",
      "Output Indices: 26\n",
      "Output Indices: 54\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 43\n",
      "Output Indices: 6\n",
      "Output Indices: 40\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 44\n",
      "Output Indices: 28\n",
      "Output Indices: 56\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 45\n",
      "Output Indices: 8\n",
      "Output Indices: 42\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 46\n",
      "Output Indices: 0\n",
      "Output Indices: 30\n",
      "Output Indices: 58\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 47\n",
      "Output Indices: 10\n",
      "Output Indices: 44\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 48\n",
      "Output Indices: 2\n",
      "Output Indices: 60\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 49\n",
      "Output Indices: 12\n",
      "Output Indices: 46\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 50\n",
      "Output Indices: 4\n",
      "Output Indices: 62\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 51\n",
      "Output Indices: 14\n",
      "Output Indices: 48\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 52\n",
      "Output Indices: 6\n",
      "Output Indices: 64\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 53\n",
      "Output Indices: 16\n",
      "Output Indices: 50\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 54\n",
      "Output Indices: 8\n",
      "Output Indices: 66\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 55\n",
      "Output Indices: 18\n",
      "Output Indices: 52\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 56\n",
      "Output Indices: 10\n",
      "Output Indices: 38\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 57\n",
      "Output Indices: 20\n",
      "Output Indices: 54\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 58\n",
      "Output Indices: 12\n",
      "Output Indices: 40\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 59\n",
      "Output Indices: 22\n",
      "Output Indices: 56\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 60\n",
      "Output Indices: 14\n",
      "Output Indices: 42\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 61\n",
      "Output Indices: 24\n",
      "Output Indices: 58\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 62\n",
      "Output Indices: 16\n",
      "Output Indices: 44\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 63\n",
      "Output Indices: 26\n",
      "Output Indices: 60\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 64\n",
      "Output Indices: 18\n",
      "Output Indices: 46\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 65\n",
      "Output Indices: 0\n",
      "Output Indices: 28\n",
      "Output Indices: 62\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 66\n",
      "Output Indices: 20\n",
      "Output Indices: 48\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 67\n",
      "Output Indices: 2\n",
      "Output Indices: 30\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 68\n",
      "Output Indices: 22\n",
      "Output Indices: 50\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 69\n",
      "Output Indices: 4\n",
      "Output Indices: 32\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 70\n",
      "Output Indices: 24\n",
      "Output Indices: 52\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 71\n",
      "Output Indices: 6\n",
      "Output Indices: 34\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 72\n",
      "Output Indices: 26\n",
      "Output Indices: 54\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 73\n",
      "Output Indices: 8\n",
      "Output Indices: 36\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 74\n",
      "Output Indices: 0\n",
      "Output Indices: 28\n",
      "Output Indices: 56\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 75\n",
      "Output Indices: 10\n",
      "Output Indices: 38\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 76\n",
      "Output Indices: 2\n",
      "Output Indices: 30\n",
      "Output Indices: 58\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 77\n",
      "Output Indices: 12\n",
      "Output Indices: 40\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 78\n",
      "Output Indices: 4\n",
      "Output Indices: 60\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 79\n",
      "Output Indices: 14\n",
      "Output Indices: 42\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 80\n",
      "Output Indices: 6\n",
      "Output Indices: 62\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 81\n",
      "Output Indices: 16\n",
      "Output Indices: 44\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 82\n",
      "Output Indices: 8\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 83\n",
      "Output Indices: 18\n",
      "Output Indices: 46\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 84\n",
      "Output Indices: 10\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 85\n",
      "Output Indices: 20\n",
      "Output Indices: 48\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 86\n",
      "Output Indices: 12\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 87\n",
      "Output Indices: 22\n",
      "Output Indices: 50\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 88\n",
      "Output Indices: 14\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 89\n",
      "Output Indices: 24\n",
      "Output Indices: 52\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 90\n",
      "Output Indices: 0\n",
      "Output Indices: 16\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 91\n",
      "Output Indices: 26\n",
      "Output Indices: 38\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 92\n",
      "Output Indices: 2\n",
      "Output Indices: 18\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 93\n",
      "Output Indices: 28\n",
      "Output Indices: 40\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 94\n",
      "Output Indices: 4\n",
      "Output Indices: 20\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 95\n",
      "Output Indices: 30\n",
      "Output Indices: 42\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 96\n",
      "Output Indices: 6\n",
      "Output Indices: 22\n",
      "Output Indices: 38\n",
      "Output Indices: 54\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 97\n",
      "Output Indices: 32\n",
      "Output Indices: 44\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 98\n",
      "Output Indices: 8\n",
      "Output Indices: 24\n",
      "Output Indices: 40\n",
      "Output Indices: 56\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 99\n",
      "Output Indices: 34\n",
      "Output Indices: 46\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 100\n",
      "Output Indices: 10\n",
      "Output Indices: 26\n",
      "Output Indices: 42\n",
      "Output Indices: 58\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 101\n",
      "Output Indices: 36\n",
      "Output Indices: 48\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 102\n",
      "Output Indices: 12\n",
      "Output Indices: 28\n",
      "Output Indices: 44\n",
      "Output Indices: 60\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 103\n",
      "Output Indices: 38\n",
      "Output Indices: 50\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 104\n",
      "Output Indices: 14\n",
      "Output Indices: 30\n",
      "Output Indices: 46\n",
      "Output Indices: 62\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 105\n",
      "Output Indices: 40\n",
      "Output Indices: 52\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 106\n",
      "Output Indices: 16\n",
      "Output Indices: 32\n",
      "Output Indices: 48\n",
      "Output Indices: 64\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 107\n",
      "Output Indices: 42\n",
      "Output Indices: 54\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 108\n",
      "Output Indices: 18\n",
      "Output Indices: 34\n",
      "Output Indices: 50\n",
      "Output Indices: 66\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 109\n",
      "Output Indices: 44\n",
      "Output Indices: 56\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 110\n",
      "Output Indices: 20\n",
      "Output Indices: 36\n",
      "Output Indices: 52\n",
      "Output Indices: 68\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 111\n",
      "Output Indices: 46\n",
      "Output Indices: 58\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 112\n",
      "Output Indices: 22\n",
      "Output Indices: 38\n",
      "Output Indices: 54\n",
      "Output Indices: 70\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 113\n",
      "Output Indices: 48\n",
      "Output Indices: 60\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 114\n",
      "Output Indices: 24\n",
      "Output Indices: 40\n",
      "Output Indices: 56\n",
      "Output Indices: 72\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 115\n",
      "Output Indices: 50\n",
      "Output Indices: 62\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 116\n",
      "Output Indices: 26\n",
      "Output Indices: 42\n",
      "Output Indices: 58\n",
      "Output Indices: 74\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 117\n",
      "Output Indices: 52\n",
      "Output Indices: 64\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 118\n",
      "Output Indices: 28\n",
      "Output Indices: 44\n",
      "Output Indices: 60\n",
      "Output Indices: 76\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 119\n",
      "Output Indices: 54\n",
      "Output Indices: 66\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 120\n",
      "Output Indices: 30\n",
      "Output Indices: 46\n",
      "Output Indices: 62\n",
      "Output Indices: 78\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 121\n",
      "Output Indices: 56\n",
      "Output Indices: 68\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 122\n",
      "Output Indices: 32\n",
      "Output Indices: 48\n",
      "Output Indices: 64\n",
      "Output Indices: 80\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 123\n",
      "Output Indices: 58\n",
      "Output Indices: 70\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 124\n",
      "Output Indices: 34\n",
      "Output Indices: 50\n",
      "Output Indices: 66\n",
      "Output Indices: 82\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 125\n",
      "Output Indices: 60\n",
      "Output Indices: 72\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 126\n",
      "Output Indices: 36\n",
      "Output Indices: 52\n",
      "Output Indices: 68\n",
      "Output Indices: 84\n",
      "-----------------------------------------------------------------------------------------\n",
      "Fault target 127\n",
      "Output Indices: 62\n",
      "Output Indices: 74\n",
      "-----------------------------------------------------------------------------------------\n",
      "end\n"
     ]
    }
   ],
   "source": [
    "from IPython import get_ipython\n",
    "get_ipython().magic('reset -sf')\n",
    "from sage.crypto.boolean_function import BooleanFunction\n",
    "from random import randint\n",
    "import random\n",
    "import numpy\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> = 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",
    "#Generate pre-output Z\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",
    "    #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",
    "    return z\n",
    "\n",
    "num_test=100\n",
    "keystream_round = 200\n",
    "success = [[0 for i in range(128)] for j in range(keystream_round)]\n",
    "#print(success)\n",
    "for i in range(num_test): #Test for 100 random initial state\n",
    "    print('Test',i)\n",
    "    #initialize the internal states with random value (random Key and IV)\n",
    "    S_temp = [R(random.randint(0,1)) for u in range(128)]\n",
    "    B_temp = [R(random.randint(0,1)) for u in range(128,256,1)]\n",
    "    S = [C() for u in range(128)]\n",
    "    B = [C() for u in range(128)]\n",
    "    statef_S = [C() for u in range(128)]\n",
    "    statef_B = [C() for u in range(128,256,1)]\n",
    "    \n",
    "    #Fault target location\n",
    "    for faulty_state in range(9,128, 1):\n",
    "        #Restore the random state\n",
    "        for l in range(128):\n",
    "            S[l] = S_temp[l]\n",
    "            B[l] = B_temp[l]\n",
    "            statef_S[l] = S[l]\n",
    "            statef_B[l] = B[l]\n",
    "            \n",
    " \n",
    "        #Inject fault at the target location\n",
    "        #statef_S[faulty_state] = statef_S[faulty_state] + 1\n",
    "        statef_B[faulty_state] = statef_B[faulty_state] + 1\n",
    "        \n",
    "        #Generate the keystream for the desired number of rounds\n",
    "        for target_output in range(keystream_round):\n",
    "            zff = grain(S, B)\n",
    "            zf = grain(statef_S, statef_B)\n",
    "            if target_output%2==0:\n",
    "            #if keystream >=1 and keystream%2==0:\n",
    "                d = zf + zff\n",
    "                if d == 1:\n",
    "                    success[target_output][faulty_state] = success[target_output][faulty_state]+1\n",
    "\n",
    "#print('Unknown register location:', unknowns)\n",
    "#print('Success count:',success)\n",
    "for j in range(128):\n",
    "    print('Fault target', j)\n",
    "    for i in range(keystream_round):\n",
    "        if success[i][j] ==num_test:\n",
    "            print('Output Indices:',i)\n",
    "            #print('Fault target',j)\n",
    "    #print('Fault target', j)\n",
    "    #print('Keystream indices:', keystream_indices)\n",
    "    print('-----------------------------------------------------------------------------------------')\n",
    "       \n",
    "print('end')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "SageMath 9.0",
   "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
}