hillclimb_implement_generalized_collection.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "f01dd075",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "from sympy import Matrix, init_printing\n",
    "import numpy as np\n",
    "import random"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "972212a0",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pickle\n",
    "with open(\"trigram.pickle\",\"rb\") as f:\n",
    "    trigram = pickle.load(f)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "79375f36",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pickle\n",
    "with open(\"twogram.pickle\",\"rb\") as f:\n",
    "    twogram = pickle.load(f)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "dd07d53d",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pickle\n",
    "with open(\"onegram.pickle\",\"rb\") as f:\n",
    "    onegram = pickle.load(f)/100"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "id": "ba438e8c",
   "metadata": {},
   "outputs": [],
   "source": [
    "def decrypt_cipher(ciphertext,key):\n",
    "    alphabet = string.ascii_lowercase\n",
    "    plaintext = ciphertext.translate(str.maketrans(alphabet,key))\n",
    "    return plaintext\n",
    "\n",
    "\n",
    "def generate_di_freq_matrix(ciphertext, key):\n",
    "    \n",
    "    alphabet = string.ascii_lowercase\n",
    "    plaintext = decrypt_cipher(ciphertext,key)\n",
    "\n",
    "    freq_mat = np.zeros(shape=(26,26))\n",
    "\n",
    "    for i in range(len(plaintext)-1):\n",
    "\n",
    "        if (plaintext[i] in alphabet) and (plaintext[i+1] in alphabet):\n",
    "            freq_mat[ord(plaintext[i+1])-97][ord(plaintext[i])-97]+=1\n",
    "            \n",
    "    freq_mat = freq_mat/freq_mat.sum()\n",
    "    return freq_mat\n",
    "\n",
    "\n",
    "def generate_trig_freq_matrix(ciphertext, key):\n",
    "    \n",
    "    alphabet = string.ascii_lowercase\n",
    "    plaintext = decrypt_cipher(ciphertext,key)\n",
    "\n",
    "    freq_mat = np.zeros(shape=(26,26,26))\n",
    "\n",
    "    for i in range(len(plaintext)-2):\n",
    "\n",
    "        if (plaintext[i] in alphabet) and (plaintext[i+1] in alphabet) and (plaintext[i+2] in alphabet):\n",
    "            freq_mat[ord(plaintext[i])-97][ord(plaintext[i+1])-97][ord(plaintext[i+2])-97]+=1\n",
    "            \n",
    "    freq_mat = freq_mat/freq_mat.sum()\n",
    "    return freq_mat\n",
    "\n",
    "def generate_quad_freq_matrix(ciphertext, key):\n",
    "    \n",
    "    alphabet = string.ascii_lowercase\n",
    "    plaintext = decrypt_cipher(ciphertext,key)\n",
    "\n",
    "    freq_mat = np.zeros(shape=(26,26,26,26))\n",
    "\n",
    "    for i in range(len(plaintext)-3):\n",
    "\n",
    "        if (plaintext[i] in alphabet) and (plaintext[i+1] in alphabet) and (plaintext[i+2] in alphabet)and (plaintext[i+3] in alphabet):\n",
    "            freq_mat[ord(plaintext[i])-97][ord(plaintext[i+1])-97][ord(plaintext[i+2])-97][ord(plaintext[i+3])-97]+=1\n",
    "            \n",
    "    freq_mat = freq_mat/freq_mat.sum()\n",
    "    return freq_mat\n",
    "\n",
    "\n",
    "def swap_chars(st,a,b):\n",
    "    result = \"\"\n",
    "    for i in st:\n",
    "        if i == a:\n",
    "            result += b\n",
    "        elif i == b:\n",
    "            result += a\n",
    "        else:\n",
    "            result +=i\n",
    "    return result\n",
    "\n",
    "def fswap_dig(freq_mat,a,b):\n",
    "    freq_mat = freq_mat.copy()\n",
    "    a_cord = ord(a)-97\n",
    "    b_cord = ord(b)-97\n",
    "    freq_mat[[a_cord,b_cord]]=freq_mat[[b_cord,a_cord]]\n",
    "    freq_mat[:,[a_cord,b_cord]]=freq_mat[:,[b_cord,a_cord]]\n",
    "    return freq_mat\n",
    "\n",
    "def fswap_trig(freq_mat,a,b):\n",
    "    freq_mat = freq_mat.copy()\n",
    "    a_cord = ord(a)-97\n",
    "    b_cord = ord(b)-97\n",
    "    freq_mat[[a_cord,b_cord]]=freq_mat[[b_cord,a_cord]]\n",
    "    freq_mat[:,[a_cord,b_cord]]=freq_mat[:,[b_cord,a_cord]]\n",
    "    freq_mat[:,:,[a_cord,b_cord]]=freq_mat[:,:,[b_cord,a_cord]]\n",
    "    return freq_mat\n",
    "\n",
    "def get_score(freq_mat,ngram):\n",
    "    return np.abs(freq_mat-ngram).sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "id": "6d444549",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "nopqrstuvwxydabcjezghifklm a very long sentence that quid quid solver is unable to jecrydt ha ha because it is quite jumb 99.37424843040051\n",
      "spfurmtawqxvihdcbezygljkon h werv opsy mestesce tght uali uali mpower lm ashdoe tp becrvit gh gh dechame lt lm ualte band 99.35407472049451\n",
      "spgfrmuiwqxavhdtbezycnjkol h wera opsy mesueste uchu finv finv mpower nm ishdoe up betravu ch ch dethime nu nm finue bild 99.35188907649125\n",
      "spvfrmuiwqxaghdtbezycnjkol h wera opsy mesueste uchu fing fing mpower nm ishdoe up betragu ch ch dethime nu nm finue bild 99.28686106833656\n",
      "spvkrmuiwqxaghdtbezycnjfol h wera opsy mesueste uchu king king mpower nm ishdoe up betragu ch ch dethime nu nm kinue bild 99.2779887908495\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[1;32m~\\AppData\\Local\\Temp/ipykernel_2484/3825719912.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m\u001b[0m\n\u001b[0;32m     18\u001b[0m     \u001b[0mputative_key\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mswap_chars\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     19\u001b[0m     \u001b[0mputative_freq_mat\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mfswap_trig\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfreq_mat\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 20\u001b[1;33m     \u001b[0mputative_score\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mget_score\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mputative_freq_mat\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mtrigram\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m     21\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     22\u001b[0m     \u001b[1;32mif\u001b[0m \u001b[0mk\u001b[0m\u001b[1;33m>\u001b[0m\u001b[1;36m10000\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m~\\AppData\\Local\\Temp/ipykernel_2484/2680380359.py\u001b[0m in \u001b[0;36mget_score\u001b[1;34m(freq_mat, ngram)\u001b[0m\n\u001b[0;32m     66\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     67\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mget_score\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfreq_mat\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mngram\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 68\u001b[1;33m     \u001b[1;32mreturn\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mabs\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfreq_mat\u001b[0m\u001b[1;33m-\u001b[0m\u001b[0mngram\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msum\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "##slow version\n",
    "import string\n",
    "import numpy as np\n",
    "from sympy import Matrix, init_printing\n",
    "init_printing()\n",
    "ciphertext=\"n irel ybat fragrapr gung dhvm dhvm fbyire vf hanoyr gb qrpelmg un un orpnhfr vg vf dhvgr qhzo\"\n",
    "# key = string.ascii_lowercase\n",
    "key = \"nopqrstuvwxyzabcdefghijklm\"\n",
    "freq_mat = generate_trig_freq_matrix(ciphertext,key)\n",
    "c_score=get_score(freq_mat,trigram)\n",
    "g_score=c_score\n",
    "k = 0\n",
    "tot=0\n",
    "while True:\n",
    "    k+=1\n",
    "    tot+=1\n",
    "    a,b=np.random.choice(list(string.ascii_lowercase),replace=False,size=2,p=onegram)\n",
    "    putative_key = swap_chars(key,a,b)\n",
    "    putative_freq_mat = fswap_trig(freq_mat,a,b)\n",
    "    putative_score = get_score(putative_freq_mat,trigram)\n",
    "    \n",
    "    if k>10000:\n",
    "        key_=list(key)\n",
    "        random.shuffle(key_)\n",
    "        key=''.join(key_)\n",
    "        freq_mat = generate_trig_freq_matrix(ciphertext,key)\n",
    "        c_score=get_score(freq_mat,trigram)    \n",
    "        k=0\n",
    "        \n",
    "    if putative_score < c_score:\n",
    "        key=putative_key\n",
    "        freq_mat = putative_freq_mat\n",
    "        c_score = putative_score\n",
    "        k=0\n",
    "        #print(key,decrypt_cipher(ciphertext,key))\n",
    "    if putative_score < g_score:\n",
    "        gkey=putative_key\n",
    "        g_score=putative_score\n",
    "        print(key,decrypt_cipher(ciphertext,gkey),g_score)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "id": "205aef1e",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "nopqrstuvwxyzabcjefghidklm a very long sentence that quiz quiz solver is unable to jecryzt ha ha because it is quite jumb 2599.501850541988\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[1;32m~\\AppData\\Local\\Temp/ipykernel_2484/2355962410.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m\u001b[0m\n\u001b[0;32m     17\u001b[0m     \u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrandom\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mchoice\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mlist\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mstring\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mascii_lowercase\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mreplace\u001b[0m\u001b[1;33m=\u001b[0m\u001b[1;32mFalse\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0msize\u001b[0m\u001b[1;33m=\u001b[0m\u001b[1;36m2\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0monegram\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     18\u001b[0m     \u001b[0mputative_key\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mswap_chars\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 19\u001b[1;33m     \u001b[0mputative_freq_mat\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mgenerate_quad_freq_matrix\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mciphertext\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mputative_key\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m     20\u001b[0m     \u001b[0mputative_score\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mget_score\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mputative_freq_mat\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mtrigram\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     21\u001b[0m     \u001b[1;32mif\u001b[0m \u001b[0mk\u001b[0m\u001b[1;33m>\u001b[0m\u001b[1;36m10000\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m~\\AppData\\Local\\Temp/ipykernel_2484/476348549.py\u001b[0m in \u001b[0;36mgenerate_quad_freq_matrix\u001b[1;34m(ciphertext, key)\u001b[0m\n\u001b[0;32m     48\u001b[0m             \u001b[0mfreq_mat\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mord\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mplaintext\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m-\u001b[0m\u001b[1;36m97\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mord\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mplaintext\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m+\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m-\u001b[0m\u001b[1;36m97\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mord\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mplaintext\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m+\u001b[0m\u001b[1;36m2\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m-\u001b[0m\u001b[1;36m97\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mord\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mplaintext\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m+\u001b[0m\u001b[1;36m3\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m-\u001b[0m\u001b[1;36m97\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m+=\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     49\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 50\u001b[1;33m     \u001b[0mfreq_mat\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mfreq_mat\u001b[0m\u001b[1;33m/\u001b[0m\u001b[0mfreq_mat\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msum\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m     51\u001b[0m     \u001b[1;32mreturn\u001b[0m \u001b[0mfreq_mat\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     52\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m~\\AppData\\Local\\Programs\\Python\\Python39\\lib\\site-packages\\numpy\\core\\_methods.py\u001b[0m in \u001b[0;36m_sum\u001b[1;34m(a, axis, dtype, out, keepdims, initial, where)\u001b[0m\n\u001b[0;32m     46\u001b[0m def _sum(a, axis=None, dtype=None, out=None, keepdims=False,\n\u001b[0;32m     47\u001b[0m          initial=_NoValue, where=True):\n\u001b[1;32m---> 48\u001b[1;33m     \u001b[1;32mreturn\u001b[0m \u001b[0mumr_sum\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0maxis\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mdtype\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mout\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mkeepdims\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0minitial\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mwhere\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m     49\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     50\u001b[0m def _prod(a, axis=None, dtype=None, out=None, keepdims=False,\n",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "##slow version\n",
    "import string\n",
    "import numpy as np\n",
    "from sympy import Matrix, init_printing\n",
    "init_printing()\n",
    "ciphertext=\"n irel ybat fragrapr gung dhvm dhvm fbyire vf hanoyr gb qrpelmg un un orpnhfr vg vf dhvgr qhzo\"\n",
    "# key = string.ascii_lowercase\n",
    "key = \"nopqrstuvwxyzabcdefghijklm\"\n",
    "freq_mat = generate_quad_freq_matrix(ciphertext,key)\n",
    "c_score=get_score(freq_mat,trigram)\n",
    "g_score=c_score\n",
    "k = 0\n",
    "tot=0\n",
    "while True:\n",
    "    k+=1\n",
    "    tot+=1\n",
    "    a,b=np.random.choice(list(string.ascii_lowercase),replace=False,size=2,p=onegram)\n",
    "    putative_key = swap_chars(key,a,b)\n",
    "    putative_freq_mat = generate_quad_freq_matrix(ciphertext,putative_key)\n",
    "    putative_score = get_score(putative_freq_mat,trigram)\n",
    "    if k>10000:\n",
    "        key_=list(key)\n",
    "        random.shuffle(key_)\n",
    "        key=''.join(key_)\n",
    "        freq_mat = generate_quad_freq_matrix(ciphertext,key)\n",
    "        c_score=get_score(freq_mat,trigram)    \n",
    "        k=0\n",
    "        \n",
    "    if putative_score < c_score:\n",
    "        key=putative_key\n",
    "        freq_mat = putative_freq_mat\n",
    "        c_score = putative_score\n",
    "        k=0\n",
    "        #print(key,decrypt_cipher(ciphertext,key))\n",
    "    if putative_score < g_score:\n",
    "        gkey=putative_key\n",
    "        g_score=putative_score\n",
    "        print(key,decrypt_cipher(ciphertext,gkey),g_score)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "efc958ca",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "pbhreavwjumqgcdnysxitkfolz nmivcprhkwihjihsijfnjauoeauoewrpmivowuhnypijrxisvcejfnfnyisnuwiojowauojixuby\n",
      "pbhrnavwjumqgcdeysxitkfolz emivcprhkwihjihsijfejauonauonwrpmivowuheypijrxisvcnjfefeyiseuwiojowauojixuby\n",
      "pbhrnavwjdmqgcueysxitkfolz emivcprhkwihjihsijfejadonadonwrpmivowdheypijrxisvcnjfefeyisedwiojowadojixdby\n",
      "pbhrnavwjdmqgtueysxickfolz emivtprhkwihjihsijfejadonadonwrpmivowdheypijrxisvtnjfefeyisedwiojowadojixdby\n",
      "pbhanrvwjdmqgtueysxickfolz emivtpahkwihjihsijfejrdonrdonwapmivowdheypijaxisvtnjfefeyisedwiojowrdojixdby\n",
      "pbhanrvwjdmqgsueytxickfolz emivspahkwihjihtijfejrdonrdonwapmivowdheypijaxitvsnjfefeyitedwiojowrdojixdby\n",
      "pbmanrvwjdhqgsueytxickfolz ehivspamkwimjimtijfejrdonrdonwaphivowdmeypijaxitvsnjfefeyitedwiojowrdojixdby\n",
      "pbmanrvwjdhqgsfeytxickuolz ehivspamkwimjimtijuejrdonrdonwaphivowdmeypijaxitvsnjueueyitedwiojowrdojixdby\n",
      "pbcanrvwjdhqgsfeytximkuolz ehivspackwicjictijuejrdonrdonwaphivowdceypijaxitvsnjueueyitedwiojowrdojixdby\n",
      "pbcinrvwjdhqgsfeytxamkuolz ehavspickwacjactajuejrdonrdonwiphavowdceypajixatvsnjueueyatedwaojowrdojaxdby\n",
      "pblinrvwjdhqgsfeytxamkuocz ehavspilkwaljaltajuejrdonrdonwiphavowdleypajixatvsnjueueyatedwaojowrdojaxdby\n",
      "pblinrvwjthqgsfeydxamkuocz ehavspilkwaljaldajuejrtonrtonwiphavowtleypajixadvsnjueueyadetwaojowrtojaxtby\n",
      "pblinsvwjthqgrfeydxamkuocz ehavrpilkwaljaldajuejstonstonwiphavowtleypajixadvrnjueueyadetwaojowstojaxtby\n",
      "pblinevwjthqgrfsydxamkuocz shavrpilkwaljaldajusjetonetonwiphavowtlsypajixadvrnjususyadstwaojowetojaxtby\n",
      "pblinrvwjthqgefsydxamkuocz shavepilkwaljaldajusjrtonrtonwiphavowtlsypajixadvenjususyadstwaojowrtojaxtby\n",
      "pblinsvwjthqgefrydxamkuocz rhavepilkwaljaldajurjstonstonwiphavowtlrypajixadvenjururyadrtwaojowstojaxtby\n",
      "pblinsvwdthqgefryjxamkuocz rhavepilkwaldaljadurdstonstonwiphavowtlrypadixajvendururyajrtwaodowstodaxtby\n",
      "pblirsvwdthqgefnyjxamkuocz nhavepilkwaldaljadundstorstorwiphavowtlnypadixajverdununyajntwaodowstodaxtby\n",
      "pblirsvwdthqgafnyjxemkuocz nhevapilkweldeljedundstorstorwiphevowtlnypedixejvardununyejntweodowstodextby\n",
      "pbnirsvwdthqgaflyjxemkuocz lhevapinkwendenjeduldstorstorwiphevowtnlypedixejvardululyejltweodowstodextby\n",
      "pbnirsvwdthqgafcyjxemkuolz chevapinkwendenjeducdstorstorwiphevowtncypedixejvarducucyejctweodowstodextby\n",
      "pbnirsvwmthqgafcyjxedkuolz chevapinkwenmenjemucmstorstorwiphevowtncypemixejvarmucucyejctweomowstomextby\n",
      "pbnirsfwmthqgavcyjxedkuolz chefapinkwenmenjemucmstorstorwiphefowtncypemixejfarmucucyejctweomowstomextby\n",
      "plnirsfwmthqgavcyjxedkuobz chefapinkwenmenjemucmstorstorwiphefowtncypemixejfarmucucyejctweomowstomextly\n",
      "plniasfwmthqgrvcyjxedkuobz chefrpinkwenmenjemucmstoastoawiphefowtncypemixejframucucyejctweomowstomextly\n",
      "pwniasflmthqgrvcyjxedkuobz chefrpinklenmenjemucmstoastoaliphefoltncypemixejframucucyejctleomolstomextwy\n",
      "pwniasfymthqgrvcljxedkuobz chefrpinkyenmenjemucmstoastoayiphefoytnclpemixejframucuclejctyeomoystomextwl\n",
      "pwniacfymthqgrvsljxedkuobz shefrpinkyenmenjemusmctoactoayiphefoytnslpemixejframususlejstyeomoyctomextwl\n",
      "pwniacdymthqgrvsljxefkuobz shedrpinkyenmenjemusmctoactoayiphedoytnslpemixejdramususlejstyeomoyctomextwl\n",
      "pwaincdymthqgrvsljxefkuobz shedrpiakyeameajemusmctonctonyiphedoytaslpemixejdrnmususlejstyeomoyctomextwl\n",
      "pwsincdymthqgrvaljxefkuobz ahedrpiskyesmesjemuamctonctonyiphedoytsalpemixejdrnmuaualejatyeomoyctomextwl\n",
      "pwniscdymthqgrvaljxefkuobz ahedrpinkyenmenjemuamctosctosyiphedoytnalpemixejdrsmuaualejatyeomoyctomextwl\n",
      "cwnispdymthqgrvaljxefkuobz ahedrcinkyenmenjemuamptosptosyichedoytnalcemixejdrsmuaualejatyeomoyptomextwl\n",
      "cwnilpdymthqgrvasjxefkuobz ahedrcinkyenmenjemuamptolptolyichedoytnascemixejdrlmuauasejatyeomoyptomextws\n",
      "cwniypdlmthqgrvasjxefkuobz ahedrcinklenmenjemuamptoyptoylichedoltnascemixejdrymuauasejatleomolptomextws\n",
      "cwniypdsmthqgrvaljxefkuobz ahedrcinksenmenjemuamptoyptoysichedostnalcemixejdrymuaualejatseomosptomextwl\n",
      "cwniypdsmtjqgrvalhxefkuobz ajedrcinksenmenhemuamptoyptoysicjedostnalcemixehdrymuaualehatseomosptomextwl\n",
      "cwniypdsmtjqgrvalfxehkuobz ajedrcinksenmenfemuamptoyptoysicjedostnalcemixefdrymuaualefatseomosptomextwl\n",
      "cwniypgsmtjqdrvalfxehkuobz ajegrcinksenmenfemuamptoyptoysicjegostnalcemixefgrymuaualefatseomosptomextwl\n",
      "cwriypgsmtjqdnvalfxehkuobz ajegncirksermerfemuamptoyptoysicjegostralcemixefgnymuaualefatseomosptomextwl\n",
      "dwriypgsmtjqcnvalfxehkuobz ajegndirksermerfemuamptoyptoysidjegostraldemixefgnymuaualefatseomosptomextwl\n",
      "fwriypgsmtjqcnvaldxehkuobz ajegnfirksermerdemuamptoyptoysifjegostralfemixedgnymuaualedatseomosptomextwl\n",
      "fwriypgsmtjqcnvaldxehkbouz ajegnfirksermerdembamptoyptoysifjegostralfemixedgnymbabaledatseomosptomextwl\n",
      "uwriypgsmtjqcnvaldxehkbofz ajegnuirksermerdembamptoyptoysiujegostraluemixedgnymbabaledatseomosptomextwl\n",
      "dwriypgsmtjqcnvaluxehkbofz ajegndirksermeruembamptoyptoysidjegostraldemixeugnymbabaleuatseomosptomextwl\n",
      "dwriypgsmthqcnvaluxejkbofz ahegndirksermeruembamptoyptoysidhegostraldemixeugnymbabaleuatseomosptomextwl\n",
      "cwriypgsmthqdnvaluxejkbofz ahegncirksermeruembamptoyptoysichegostralcemixeugnymbabaleuatseomosptomextwl\n",
      "cwriypgsmthqunvaldxejkbofz ahegncirksermerdembamptoyptoysichegostralcemixedgnymbabaledatseomosptomextwl\n",
      "cwliypgsmthqunvardxejkbofz ahegncilkselmeldembamptoyptoysichegostlarcemixedgnymbabaredatseomosptomextwr\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[1;32m~\\AppData\\Local\\Temp/ipykernel_2484/1420222530.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m\u001b[0m\n\u001b[0;32m     11\u001b[0m     \u001b[0mputative_key\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mswap_chars\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     12\u001b[0m     \u001b[0mputative_freq_mat\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mfswap_dig\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfreq_mat\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 13\u001b[1;33m     \u001b[0mputative_score\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mget_score\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mputative_freq_mat\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mtwogram\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m     14\u001b[0m     \u001b[1;32mif\u001b[0m \u001b[0mputative_score\u001b[0m \u001b[1;33m<\u001b[0m \u001b[0mc_score\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     15\u001b[0m         \u001b[0mkey\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mputative_key\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m~\\AppData\\Local\\Temp/ipykernel_2484/2680380359.py\u001b[0m in \u001b[0;36mget_score\u001b[1;34m(freq_mat, ngram)\u001b[0m\n\u001b[0;32m     66\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     67\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mget_score\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfreq_mat\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mngram\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 68\u001b[1;33m     \u001b[1;32mreturn\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mabs\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfreq_mat\u001b[0m\u001b[1;33m-\u001b[0m\u001b[0mngram\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msum\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "# ciphertext=\"nn va gur rirag gung gur chepunfre qrsnhygf va gur cnlzrag bs nal vafgnyzrag bs chepunfr cevpr, gnkrf, vafhenapr, vagrerfg, be gur naahny punetr qrfpevorq ryfrjurer urerva, be funyy qrsnhyg va gur cresbeznapr bs nal bgure boyvtngvbaf frg sbegu va guvf pbagenpg, gur fryyre znl: ng uvf bcgvba: (n) qrpyner vzzrqvngryl qhr naq cnlnoyr gur ragver hacnvq onynapr bs chepunfr cevpr, jvgu nppehrq vagrerfg, gnkrf, naq naahny punetr, naq qrznaq shyy cnlzrag gurerbs, naq rasbepr pbairlnapr bs gur ynaq ol grezvangvba bs gur pbagenpg be nppbeqvat gb gur grezf urerbs, va juvpu pnfr gur chepunfre funyy nyfb or yvnoyr gb gur fryyre sbe ernfbanoyr nggbearl'f srrf sbe freivprf eraqrerq ol nal nggbearl ba orunys bs gur fryyre, be (o) fryy fnvq ynaq naq cerzvfrf be nal cneg gurerbs ng choyvp nhpgvba, va fhpu znaare, ng fhpu gvzr naq cynpr, hcba fhpu grezf naq pbaqvgvbaf, naq hcba fhpu choyvp abgvpr nf gur fryyre znl qrrz orfg sbe gur vagrerfg bs nyy pbaprearq, pbafvfgvat bs nqiregvfrzrag va n arjfcncre bs trareny pvephyngvba va gur pbhagl be pvgl va juvpu gur frphevgl cebcregl vf ybpngrq ng yrnfg bapr n jrrx sbe guerr (3) fhpprffvir jrrxf be sbe fhpu crevbq nf nccyvpnoyr ynj znl erdhver naq, va pnfr bs qrsnhyg bs nal chepunfre, gb er-fryy jvgu fhpu cbfgcbarzrag bs fnyr be erfnyr naq hcba fhpu choyvp abgvpr gurerbs nf gur fryyre znl qrgrezvar, naq hcba pbzcyvnapr ol gur chepunfre jvgu gur grezf bs fnyr, naq hcba whqvpvny nccebiny nf znl or erdhverq ol ynj, pbairl fnvq ynaq naq cerzvfrf va srr fvzcyr gb naq ng gur pbfg bs gur chepunfre, jub funyy abg or yvnoyr gb frr gb gur nccyvpngvba bs gur chepunfr zbarl; naq sebz gur cebprrqf bs gur fnyr: svefg gb cnl nyy cebcre pbfgf naq punetrf, vapyhqvat ohg abg yvzvgrq gb pbheg pbfgf, nqiregvfvat rkcrafrf, nhpgvbarre'f nyybjnapr, gur rkcrafrf, vs nal erdhverq gb pbeerpg nal veerthynevgl va gur gvgyr, cerzvhz sbe fryyre'f obaq, nhqvgbe'f srr, nggbearl'f srr, naq nyy bgure rkcrafrf bs fnyr bppheerq va naq nobhg gur cebgrpgvba naq rkrphgvba bs guvf pbagenpg, naq nyy zbarlf nqinaprq sbe gnkrf, nffrffzragf, vafhenapr, naq jvgu vagrerfg gurerba nf cebivqrq urerva, naq nyy gnkrf qhr hcba fnvq ynaq naq cerzvfrf ng gvzr bs fnyr, naq gb ergnva nf pbzcrafngvba n pbzzvffvba bs svir creprag (5%) ba gur nzbhag bs fnvq fnyr be fnyrf; frpbaq, gb cnl gur jubyr nzbhag gura erznvavat hacnvq bs gur cevapvcny bs fnvq pbagenpg, naq vagrerfg gurerba gb qngr bs cnlzrag, jurgure gur fnzr funyy or qhr be abg, vg orvat haqrefgbbq naq nterrq gung hcba fhpu fnyr orsber znghevgl bs gur pbagenpg gur onynapr gurerbs funyy or vzzrqvngryl qhr naq cnlnoyr; guveq, gb cnl yvraf bs erpbeq ntnvafg gur frphevgl cebcregl nppbeqvat gb gurve cevbevgl bs yvra naq gb gur rkgrag gung shaqf erznvavat va gur unaqf bs gur fryyre ner ninvynoyr; naq ynfg, gb cnl gur erznvaqre bs fnvq cebprrqf, vs nal, gb gur iraqbe, uvf urvef, crefbanyf ercerfragngvirf, fhpprffbef be nffvtaf hcba gur qryvirel naq fheeraqre gb gur iraqrr bs cbffrffvba bs gur ynaq naq cerzvfrf, yrff pbfgf naq rkprff bs bognvavat cbffrffvba\"\n",
    "# ciphertext=\"regrt eaf yofr ashiaztmoeas lwzlmomwmogfl dart mg a mtvm zb afasbnofu stmmtk yktjwtfeotl.\"\n",
    "ciphertext=\"pktgnadcvhtcitcrtiwpifjxefjxehdaktgxhjcpqatidstrgneiwpwpqtrpjhtxixhfjxitsjbq\"\n",
    "key=\"pbhrtavwjumqgcdnysxiekfolz\"\n",
    "freq_mat = generate_di_freq_matrix(ciphertext,key)\n",
    "c_score=get_score(freq_mat,twogram)\n",
    "k = 0\n",
    "while True:\n",
    "    k+=1\n",
    "    a,b=np.random.choice(list(string.ascii_lowercase),replace=False,size=2,p=onegram)\n",
    "    putative_key = swap_chars(key,a,b)\n",
    "    putative_freq_mat = fswap_dig(freq_mat,a,b)\n",
    "    putative_score = get_score(putative_freq_mat,twogram)\n",
    "    if putative_score < c_score:\n",
    "        key=putative_key\n",
    "        freq_mat = putative_freq_mat\n",
    "        c_score = putative_score\n",
    "        #print(key,putative_score,k)\n",
    "        k=0\n",
    "        print(key,decrypt_cipher(ciphertext,key))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3bdbc55a",
   "metadata": {},
   "outputs": [],
   "source": [
    "# ##slow version\n",
    "# import string\n",
    "# import numpy as np\n",
    "# from sympy import Matrix, init_printing\n",
    "# init_printing()\n",
    "# # ciphertext=\"nn va gur rirag gung gur chepunfre qrsnhygf va gur cnlzrag bs nal vafgnyzrag bs chepunfr cevpr, gnkrf, vafhenapr, vagrerfg, be gur naahny punetr qrfpevorq ryfrjurer urerva, be funyy qrsnhyg va gur cresbeznapr bs nal bgure boyvtngvbaf frg sbegu va guvf pbagenpg, gur fryyre znl: ng uvf bcgvba: (n) qrpyner vzzrqvngryl qhr naq cnlnoyr gur ragver hacnvq onynapr bs chepunfr cevpr, jvgu nppehrq vagrerfg, gnkrf, naq naahny punetr, naq qrznaq shyy cnlzrag gurerbs, naq rasbepr pbairlnapr bs gur ynaq ol grezvangvba bs gur pbagenpg be nppbeqvat gb gur grezf urerbs, va juvpu pnfr gur chepunfre funyy nyfb or yvnoyr gb gur fryyre sbe ernfbanoyr nggbearl'f srrf sbe freivprf eraqrerq ol nal nggbearl ba orunys bs gur fryyre, be (o) fryy fnvq ynaq naq cerzvfrf be nal cneg gurerbs ng choyvp nhpgvba, va fhpu znaare, ng fhpu gvzr naq cynpr, hcba fhpu grezf naq pbaqvgvbaf, naq hcba fhpu choyvp abgvpr nf gur fryyre znl qrrz orfg sbe gur vagrerfg bs nyy pbaprearq, pbafvfgvat bs nqiregvfrzrag va n arjfcncre bs trareny pvephyngvba va gur pbhagl be pvgl va juvpu gur frphevgl cebcregl vf ybpngrq ng yrnfg bapr n jrrx sbe guerr (3) fhpprffvir jrrxf be sbe fhpu crevbq nf nccyvpnoyr ynj znl erdhver naq, va pnfr bs qrsnhyg bs nal chepunfre, gb er-fryy jvgu fhpu cbfgcbarzrag bs fnyr be erfnyr naq hcba fhpu choyvp abgvpr gurerbs nf gur fryyre znl qrgrezvar, naq hcba pbzcyvnapr ol gur chepunfre jvgu gur grezf bs fnyr, naq hcba whqvpvny nccebiny nf znl or erdhverq ol ynj, pbairl fnvq ynaq naq cerzvfrf va srr fvzcyr gb naq ng gur pbfg bs gur chepunfre, jub funyy abg or yvnoyr gb frr gb gur nccyvpngvba bs gur chepunfr zbarl; naq sebz gur cebprrqf bs gur fnyr: svefg gb cnl nyy cebcre pbfgf naq punetrf, vapyhqvat ohg abg yvzvgrq gb pbheg pbfgf, nqiregvfvat rkcrafrf, nhpgvbarre'f nyybjnapr, gur rkcrafrf, vs nal erdhverq gb pbeerpg nal veerthynevgl va gur gvgyr, cerzvhz sbe fryyre'f obaq, nhqvgbe'f srr, nggbearl'f srr, naq nyy bgure rkcrafrf bs fnyr bppheerq va naq nobhg gur cebgrpgvba naq rkrphgvba bs guvf pbagenpg, naq nyy zbarlf nqinaprq sbe gnkrf, nffrffzragf, vafhenapr, naq jvgu vagrerfg gurerba nf cebivqrq urerva, naq nyy gnkrf qhr hcba fnvq ynaq naq cerzvfrf ng gvzr bs fnyr, naq gb ergnva nf pbzcrafngvba n pbzzvffvba bs svir creprag (5%) ba gur nzbhag bs fnvq fnyr be fnyrf; frpbaq, gb cnl gur jubyr nzbhag gura erznvavat hacnvq bs gur cevapvcny bs fnvq pbagenpg, naq vagrerfg gurerba gb qngr bs cnlzrag, jurgure gur fnzr funyy or qhr be abg, vg orvat haqrefgbbq naq nterrq gung hcba fhpu fnyr orsber znghevgl bs gur pbagenpg gur onynapr gurerbs funyy or vzzrqvngryl qhr naq cnlnoyr; guveq, gb cnl yvraf bs erpbeq ntnvafg gur frphevgl cebcregl nppbeqvat gb gurve cevbevgl bs yvra naq gb gur rkgrag gung shaqf erznvavat va gur unaqf bs gur fryyre ner ninvynoyr; naq ynfg, gb cnl gur erznvaqre bs fnvq cebprrqf, vs nal, gb gur iraqbe, uvf urvef, crefbanyf ercerfragngvirf, fhpprffbef be nffvtaf hcba gur qryvirel naq fheeraqre gb gur iraqrr bs cbffrffvba bs gur ynaq naq cerzvfrf, yrff pbfgf naq rkprff bs bognvavat cbffrffvba\"\n",
    "# ciphertext=\"regrt eaf yofr ashiaztmoeas lwzlmomwmogfl dart mg a mtvm zb afasbnofu stmmtk yktjwtfeotl.\"\n",
    "# # key = string.ascii_lowercase\n",
    "# key=\"abjwcnophqrltvifysmekxuzgd\"\n",
    "# freq_mat = generate_trig_freq_matrix(ciphertext,key)\n",
    "# c_score=get_score(freq_mat,trigram)\n",
    "# k = 0\n",
    "# tot=0\n",
    "# while True:\n",
    "#     k+=1\n",
    "#     tot+=1\n",
    "#     a,b=np.random.choice(list(string.ascii_lowercase),size=2)\n",
    "#     putative_key = swap_chars(key,a,b)\n",
    "#     putative_freq_mat = fswap_trig(freq_mat,a,b)\n",
    "#     putative_score = get_score(putative_freq_mat,trigram)\n",
    "#     assert np.all(generate_trig_freq_matrix(ciphertext,putative_key)==putative_freq_mat), \"holly mother of god\"\n",
    "#     if putative_score < c_score:\n",
    "#         key=putative_key\n",
    "#         freq_mat = putative_freq_mat\n",
    "#         c_score = putative_score\n",
    "#         k=0\n",
    "#         print(key)\n",
    "#         #print(key,decrypt_cipher(ciphertext,key))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "id": "e2ef7a89",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'ahman her xtra esbleunithes wouwitioitmrw yean im e indi up erespftrv sniing xgnconrhtnw.'"
      ]
     },
     "execution_count": 79,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "decrypt_cipher(ciphertext,\"epjyhrmblcgwiftzqasnvdokxu\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "b83b447f",
   "metadata": {},
   "outputs": [],
   "source": [
    "# assert np.all(generate_freq_matrix(ciphertext,swap_chars(key,\"a\",\"d\"))==fswap(freq_mat,\"a\",\"d\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 77,
   "id": "26610e08",
   "metadata": {
    "collapsed": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "abcdefghijkwmnopqrstuvlxyz\n"
     ]
    },
    {
     "data": {
      "text/latex": [
       "$\\displaystyle \\left[\\begin{array}{cccccccccccccccccccccccccc}0.00156657963446475 & 0.0 & 0.000522193211488251 & 0.00417754569190601 & 0.002088772845953 & 0.0120104438642298 & 0.00156657963446475 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.00783289817232376 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0130548302872063 & 0.0 & 0.00365535248041775 & 0.0 & 0.0 & 0.0 & 0.00469973890339426 & 0.00939947780678851 & 0.000522193211488251\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00469973890339426 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.000522193211488251 & 0.0 & 0.00156657963446475 & 0.002088772845953 & 0.000522193211488251 & 0.002088772845953 & 0.0 & 0.00417754569190601 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.00156657963446475 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0\\\\0.0120104438642298 & 0.0 & 0.00626631853785901 & 0.000522193211488251 & 0.00365535248041775 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.00574412532637076 & 0.00835509138381201 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.002088772845953 & 0.00261096605744125 & 0.000522193211488251 & 0.00469973890339426 & 0.0 & 0.0 & 0.0 & 0.012532637075718 & 0.002088772845953 & 0.0010443864229765\\\\0.00887728459530026 & 0.002088772845953 & 0.00626631853785901 & 0.0 & 0.00731070496083551 & 0.0 & 0.014621409921671 & 0.00261096605744125 & 0.0010443864229765 & 0.0 & 0.002088772845953 & 0.0010443864229765 & 0.0 & 0.0318537859007833 & 0.0 & 0.00261096605744125 & 0.00626631853785901 & 0.00626631853785901 & 0.00783289817232376 & 0.025065274151436 & 0.00574412532637076 & 0.0 & 0.0 & 0.012532637075718 & 0.0198433420365535 & 0.00522193211488251\\\\0.00417754569190601 & 0.0010443864229765 & 0.00731070496083551 & 0.0 & 0.0010443864229765 & 0.000522193211488251 & 0.00783289817232376 & 0.0010443864229765 & 0.000522193211488251 & 0.0 & 0.00156657963446475 & 0.0031331592689295 & 0.00156657963446475 & 0.0109660574412533 & 0.0 & 0.0 & 0.00626631853785901 & 0.00156657963446475 & 0.000522193211488251 & 0.00417754569190601 & 0.00156657963446475 & 0.0 & 0.0 & 0.00261096605744125 & 0.00887728459530026 & 0.0010443864229765\\\\0.000522193211488251 & 0.0 & 0.00261096605744125 & 0.000522193211488251 & 0.00678851174934726 & 0.016710182767624 & 0.0109660574412533 & 0.0031331592689295 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.00574412532637076\\\\0.000522193211488251 & 0.000522193211488251 & 0.00574412532637076 & 0.00261096605744125 & 0.0 & 0.0031331592689295 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.000522193211488251 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0010443864229765 & 0.002088772845953 & 0.000522193211488251 & 0.0 & 0.0 & 0.000522193211488251 & 0.00261096605744125 & 0.00678851174934726 & 0.000522193211488251\\\\0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0010443864229765 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0198433420365535 & 0.002088772845953 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0109660574412533 & 0.0 & 0.0 & 0.0 & 0.00261096605744125 & 0.00783289817232376 & 0.00678851174934726 & 0.00156657963446475 & 0.002088772845953 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.00365535248041775 & 0.0 & 0.0 & 0.00469973890339426 & 0.00156657963446475 & 0.00417754569190601 & 0.00574412532637076 & 0.0010443864229765 & 0.0 & 0.0 & 0.002088772845953 & 0.00417754569190601 & 0.0\\\\0.00469973890339426 & 0.0 & 0.0 & 0.0 & 0.00365535248041775 & 0.00731070496083551 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00417754569190601 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.002088772845953\\\\0.0349869451697128 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0031331592689295 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0109660574412533 & 0.00261096605744125 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0031331592689295 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.002088772845953 & 0.00626631853785901 & 0.00156657963446475 & 0.000522193211488251 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0031331592689295 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0031331592689295 & 0.0 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0\\\\0.0 & 0.000522193211488251 & 0.0 & 0.022976501305483 & 0.0130548302872063 & 0.0308093994778068 & 0.0 & 0.00365535248041775 & 0.0 & 0.0 & 0.0 & 0.0224543080939948 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.00835509138381201 & 0.00156657963446475 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.00469973890339426 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0255874673629243 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.00261096605744125 & 0.000522193211488251 & 0.00939947780678851 & 0.0172323759791123 & 0.0287206266318538 & 0.00365535248041775 & 0.0 & 0.00835509138381201 & 0.0 & 0.0 & 0.000522193211488251 & 0.00522193211488251 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.000522193211488251 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.00156657963446475 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.002088772845953 & 0.00417754569190601 & 0.00469973890339426 & 0.0 & 0.00626631853785901 & 0.0 & 0.0 & 0.0 & 0.00678851174934726 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.00678851174934726 & 0.0 & 0.00574412532637076 & 0.0 & 0.0 & 0.0 & 0.0031331592689295 & 0.000522193211488251 & 0.0\\\\0.00261096605744125 & 0.0 & 0.0 & 0.00365535248041775 & 0.016710182767624 & 0.0109660574412533 & 0.00156657963446475 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.00678851174934726 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.00678851174934726 & 0.002088772845953 & 0.00261096605744125 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.00626631853785901 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.00417754569190601 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\end{array}\\right]$"
      ],
      "text/plain": [
       "⎡0.00156657963446475           0.0           0.000522193211488251  0.004177545\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢0.000522193211488251          0.0           0.00156657963446475    0.00208877\n",
       "⎢                                                                             \n",
       "⎢ 0.0120104438642298           0.0           0.00626631853785901   0.000522193\n",
       "⎢                                                                             \n",
       "⎢0.00887728459530026    0.002088772845953    0.00626631853785901           0.0\n",
       "⎢                                                                             \n",
       "⎢0.00417754569190601    0.0010443864229765   0.00731070496083551           0.0\n",
       "⎢                                                                             \n",
       "⎢0.000522193211488251          0.0           0.00261096605744125   0.000522193\n",
       "⎢                                                                             \n",
       "⎢0.000522193211488251  0.000522193211488251  0.00574412532637076   0.002610966\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0           0.000522193\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0            0.01984334\n",
       "⎢                                                                             \n",
       "⎢ 0.0109660574412533           0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢0.00469973890339426           0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢ 0.0349869451697128           0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0            0.00208877\n",
       "⎢                                                                             \n",
       "⎢        0.0           0.000522193211488251          0.0            0.02297650\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0            0.00104438\n",
       "⎢                                                                             \n",
       "⎢0.00261096605744125   0.000522193211488251  0.00939947780678851    0.01723237\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0            0.00104438\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0            0.00208877\n",
       "⎢                                                                             \n",
       "⎢0.00261096605744125           0.0                   0.0           0.003655352\n",
       "⎢                                                                             \n",
       "⎣        0.0                   0.0                   0.0            0.00104438\n",
       "\n",
       "69190601    0.002088772845953     0.0120104438642298   0.00156657963446475   0\n",
       "                                                                              \n",
       "           0.000522193211488251   0.0010443864229765           0.0            \n",
       "                                                                              \n",
       "2845953    0.000522193211488251   0.002088772845953            0.0           0\n",
       "                                                                              \n",
       "211488251  0.00365535248041775           0.0            0.0010443864229765    \n",
       "                                                                              \n",
       "           0.00731070496083551           0.0            0.014621409921671    0\n",
       "                                                                              \n",
       "            0.0010443864229765   0.000522193211488251  0.00783289817232376    \n",
       "                                                                              \n",
       "211488251  0.00678851174934726    0.016710182767624     0.0109660574412533    \n",
       "                                                                              \n",
       "05744125           0.0            0.0031331592689295   0.000522193211488251   \n",
       "                                                                              \n",
       "211488251   0.0010443864229765    0.0010443864229765           0.0            \n",
       "                                                                              \n",
       "                   0.0                   0.0                   0.0            \n",
       "                                                                              \n",
       "20365535    0.002088772845953            0.0           0.000522193211488251   \n",
       "                                                                              \n",
       "           0.00261096605744125   0.00783289817232376   0.00678851174934726   0\n",
       "                                                                              \n",
       "           0.00365535248041775   0.00731070496083551    0.0010443864229765    \n",
       "                                                                              \n",
       "           0.000522193211488251          0.0                   0.0            \n",
       "                                                                              \n",
       "           0.00156657963446475           0.0                   0.0            \n",
       "                                                                              \n",
       "            0.0031331592689295    0.002088772845953            0.0            \n",
       "                                                                              \n",
       "2845953    0.00626631853785901   0.00156657963446475   0.000522193211488251  0\n",
       "                                                                              \n",
       "1305483     0.0130548302872063    0.0308093994778068           0.0           0\n",
       "                                                                              \n",
       "64229765   0.00835509138381201   0.00156657963446475           0.0           0\n",
       "                                                                              \n",
       "59791123    0.0287206266318538   0.00365535248041775           0.0           0\n",
       "                                                                              \n",
       "64229765   0.000522193211488251  0.000522193211488251          0.0            \n",
       "                                                                              \n",
       "            0.0010443864229765           0.0                   0.0            \n",
       "                                                                              \n",
       "                   0.0                   0.0                   0.0            \n",
       "                                                                              \n",
       "2845953    0.00417754569190601   0.00469973890339426           0.0           0\n",
       "                                                                              \n",
       "48041775    0.016710182767624     0.0109660574412533   0.00156657963446475    \n",
       "                                                                              \n",
       "64229765           0.0           0.00417754569190601           0.0           0\n",
       "\n",
       ".00156657963446475           0.0           0.0          0.0           0.007832\n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0            0.00104\n",
       "                                                                              \n",
       ".00417754569190601           0.0           0.0          0.0           0.000522\n",
       "                                                                              \n",
       "       0.0                   0.0           0.0  0.00574412532637076   0.008355\n",
       "                                                                              \n",
       ".00261096605744125    0.0010443864229765   0.0   0.002088772845953     0.00104\n",
       "                                                                              \n",
       "0.0010443864229765   0.000522193211488251  0.0  0.00156657963446475    0.00313\n",
       "                                                                              \n",
       "0.0031331592689295           0.0           0.0          0.0           0.000522\n",
       "                                                                              \n",
       "       0.0                   0.0           0.0   0.0010443864229765   0.000522\n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0                   \n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0                   \n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0            0.00104\n",
       "                                                                              \n",
       ".00156657963446475    0.002088772845953    0.0   0.0010443864229765           \n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0                   \n",
       "                                                                              \n",
       "       0.0            0.0031331592689295   0.0          0.0                   \n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0                   \n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0                   \n",
       "                                                                              \n",
       ".000522193211488251          0.0           0.0          0.0            0.00313\n",
       "                                                                              \n",
       ".00365535248041775           0.0           0.0          0.0            0.02245\n",
       "                                                                              \n",
       ".00156657963446475           0.0           0.0          0.0           0.004699\n",
       "                                                                              \n",
       ".00835509138381201           0.0           0.0  0.000522193211488251  0.005221\n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0            0.00208\n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0                   \n",
       "                                                                              \n",
       "       0.0                   0.0           0.0          0.0                   \n",
       "                                                                              \n",
       ".00626631853785901           0.0           0.0          0.0           0.006788\n",
       "                                                                              \n",
       "       0.0           0.000522193211488251  0.0          0.0           0.006788\n",
       "                                                                              \n",
       ".00156657963446475           0.0           0.0          0.0           0.000522\n",
       "\n",
       "89817232376           0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "43864229765           0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "193211488251          0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "09138381201           0.0            0.0010443864229765           0.0         \n",
       "                                                                              \n",
       "43864229765           0.0            0.0318537859007833           0.0         \n",
       "                                                                              \n",
       "31592689295   0.00156657963446475    0.0109660574412533           0.0         \n",
       "                                                                              \n",
       "193211488251          0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "193211488251          0.0                   0.0           0.00156657963446475 \n",
       "                                                                              \n",
       "0.0                   0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "0.0                   0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "43864229765           0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "0.0                   0.0           0.00365535248041775           0.0         \n",
       "                                                                              \n",
       "0.0                   0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "0.0                   0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "0.0                   0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "0.0                   0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "31592689295   0.00156657963446475           0.0                   0.0         \n",
       "                                                                              \n",
       "43080939948           0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "73890339426           0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "93211488251           0.0           0.000522193211488251          0.0         \n",
       "                                                                              \n",
       "8772845953            0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "0.0                   0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "0.0                   0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "51174934726           0.0                   0.0                   0.0         \n",
       "                                                                              \n",
       "51174934726   0.000522193211488251          0.0                   0.0         \n",
       "                                                                              \n",
       "193211488251          0.0                   0.0                   0.0         \n",
       "\n",
       " 0.000522193211488251          0.0           0.0130548302872063           0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0          0.00469973890339426           0.0 \n",
       "                                                                              \n",
       " 0.00156657963446475   0.00156657963446475   0.0010443864229765           0.0 \n",
       "                                                                              \n",
       "         0.0            0.002088772845953   0.00261096605744125   0.0005221932\n",
       "                                                                              \n",
       " 0.00261096605744125   0.00626631853785901  0.00626631853785901   0.0078328981\n",
       "                                                                              \n",
       "         0.0           0.00626631853785901  0.00156657963446475   0.0005221932\n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0           0.0010443864229765    0.002088772\n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0          0.000522193211488251          0.0 \n",
       "                                                                              \n",
       "         0.0           0.00469973890339426  0.00156657963446475   0.0041775456\n",
       "                                                                              \n",
       "         0.0                   0.0          0.00417754569190601           0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "                                                                              \n",
       "         0.0           0.00156657963446475          0.0                   0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0          0.00156657963446475           0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0           0.0255874673629243           0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0           0.0010443864229765   0.0015665796\n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "                                                                              \n",
       " 0.000522193211488251          0.0          0.00678851174934726           0.0 \n",
       "                                                                              \n",
       "         0.0           0.00156657963446475  0.00678851174934726    0.002088772\n",
       "                                                                              \n",
       "         0.0                   0.0                  0.0                   0.0 \n",
       "\n",
       "          0.00365535248041775           0.0                  0.0              \n",
       "                                                                              \n",
       "          0.00156657963446475           0.0                  0.0              \n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "                                                                              \n",
       "11488251  0.00469973890339426           0.0                  0.0              \n",
       "                                                                              \n",
       "7232376    0.025065274151436    0.00574412532637076          0.0              \n",
       "                                                                              \n",
       "11488251  0.00417754569190601   0.00156657963446475          0.0              \n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "                                                                              \n",
       "845953    0.000522193211488251          0.0                  0.0           0.0\n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "                                                                              \n",
       "          0.000522193211488251          0.0                  0.0              \n",
       "                                                                              \n",
       "9190601   0.00574412532637076   0.0010443864229765           0.0              \n",
       "                                                                              \n",
       "          0.000522193211488251          0.0                  0.0              \n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "                                                                              \n",
       "           0.0031331592689295           0.0                  0.0              \n",
       "                                                                              \n",
       "           0.002088772845953            0.0                  0.0              \n",
       "                                                                              \n",
       "           0.002088772845953            0.0                  0.0              \n",
       "                                                                              \n",
       "           0.002088772845953            0.0                  0.0              \n",
       "                                                                              \n",
       "3446475   0.000522193211488251          0.0                  0.0              \n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "                                                                              \n",
       "          0.00574412532637076           0.0                  0.0              \n",
       "                                                                              \n",
       "845953    0.00261096605744125           0.0          0.000522193211488251     \n",
       "                                                                              \n",
       "                  0.0                   0.0                  0.0              \n",
       "\n",
       "     0.0           0.00469973890339426  0.00939947780678851   0.00052219321148\n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0          0.000522193211488251          0.0     \n",
       "                                                                              \n",
       "     0.0            0.012532637075718    0.002088772845953     0.0010443864229\n",
       "                                                                              \n",
       "     0.0            0.012532637075718    0.0198433420365535   0.00522193211488\n",
       "                                                                              \n",
       "     0.0           0.00261096605744125  0.00887728459530026    0.0010443864229\n",
       "                                                                              \n",
       "     0.0           0.0010443864229765           0.0           0.00574412532637\n",
       "                                                                              \n",
       "00522193211488251  0.00261096605744125  0.00678851174934726   0.00052219321148\n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0            0.002088772845953   0.00417754569190601           0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0            0.0020887728459\n",
       "                                                                              \n",
       "     0.0           0.0109660574412533   0.00261096605744125           0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0          0.000522193211488251          0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0           0.0010443864229765           0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "                                                                              \n",
       "     0.0           0.0031331592689295   0.000522193211488251          0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0          0.00626631853785901           0.0     \n",
       "                                                                              \n",
       "     0.0                   0.0                  0.0                   0.0     \n",
       "\n",
       "8251⎤\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "765 ⎥\n",
       "    ⎥\n",
       "251 ⎥\n",
       "    ⎥\n",
       "765 ⎥\n",
       "    ⎥\n",
       "076 ⎥\n",
       "    ⎥\n",
       "8251⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "53  ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎦"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import string\n",
    "import numpy as np\n",
    "from sympy import Matrix, init_printing\n",
    "init_printing()\n",
    "ciphertext=\"nn va gur rirag gung gur chepunfre qrsnhygf va gur cnlzrag bs nal vafgnyzrag bs chepunfr cevpr, gnkrf, vafhenapr, vagrerfg, be gur naahny punetr qrfpevorq ryfrjurer urerva, be funyy qrsnhyg va gur cresbeznapr bs nal bgure boyvtngvbaf frg sbegu va guvf pbagenpg, gur fryyre znl: ng uvf bcgvba: (n) qrpyner vzzrqvngryl qhr naq cnlnoyr gur ragver hacnvq onynapr bs chepunfr cevpr, jvgu nppehrq vagrerfg, gnkrf, naq naahny punetr, naq qrznaq shyy cnlzrag gurerbs, naq rasbepr pbairlnapr bs gur ynaq ol grezvangvba bs gur pbagenpg be nppbeqvat gb gur grezf urerbs, va juvpu pnfr gur chepunfre funyy nyfb or yvnoyr gb gur fryyre sbe ernfbanoyr nggbearl'f srrf sbe freivprf eraqrerq ol nal nggbearl ba orunys bs gur fryyre, be (o) fryy fnvq ynaq naq cerzvfrf be nal cneg gurerbs ng choyvp nhpgvba, va fhpu znaare, ng fhpu gvzr naq cynpr, hcba fhpu grezf naq pbaqvgvbaf, naq hcba fhpu choyvp abgvpr nf gur fryyre znl qrrz orfg sbe gur vagrerfg bs nyy pbaprearq, pbafvfgvat bs nqiregvfrzrag va n arjfcncre bs trareny pvephyngvba va gur pbhagl be pvgl va juvpu gur frphevgl cebcregl vf ybpngrq ng yrnfg bapr n jrrx sbe guerr (3) fhpprffvir jrrxf be sbe fhpu crevbq nf nccyvpnoyr ynj znl erdhver naq, va pnfr bs qrsnhyg bs nal chepunfre, gb er-fryy jvgu fhpu cbfgcbarzrag bs fnyr be erfnyr naq hcba fhpu choyvp abgvpr gurerbs nf gur fryyre znl qrgrezvar, naq hcba pbzcyvnapr ol gur chepunfre jvgu gur grezf bs fnyr, naq hcba whqvpvny nccebiny nf znl or erdhverq ol ynj, pbairl fnvq ynaq naq cerzvfrf va srr fvzcyr gb naq ng gur pbfg bs gur chepunfre, jub funyy abg or yvnoyr gb frr gb gur nccyvpngvba bs gur chepunfr zbarl; naq sebz gur cebprrqf bs gur fnyr: svefg gb cnl nyy cebcre pbfgf naq punetrf, vapyhqvat ohg abg yvzvgrq gb pbheg pbfgf, nqiregvfvat rkcrafrf, nhpgvbarre'f nyybjnapr, gur rkcrafrf, vs nal erdhverq gb pbeerpg nal veerthynevgl va gur gvgyr, cerzvhz sbe fryyre'f obaq, nhqvgbe'f srr, nggbearl'f srr, naq nyy bgure rkcrafrf bs fnyr bppheerq va naq nobhg gur cebgrpgvba naq rkrphgvba bs guvf pbagenpg, naq nyy zbarlf nqinaprq sbe gnkrf, nffrffzragf, vafhenapr, naq jvgu vagrerfg gurerba nf cebivqrq urerva, naq nyy gnkrf qhr hcba fnvq ynaq naq cerzvfrf ng gvzr bs fnyr, naq gb ergnva nf pbzcrafngvba n pbzzvffvba bs svir creprag (5%) ba gur nzbhag bs fnvq fnyr be fnyrf; frpbaq, gb cnl gur jubyr nzbhag gura erznvavat hacnvq bs gur cevapvcny bs fnvq pbagenpg, naq vagrerfg gurerba gb qngr bs cnlzrag, jurgure gur fnzr funyy or qhr be abg, vg orvat haqrefgbbq naq nterrq gung hcba fhpu fnyr orsber znghevgl bs gur pbagenpg gur onynapr gurerbs funyy or vzzrqvngryl qhr naq cnlnoyr; guveq, gb cnl yvraf bs erpbeq ntnvafg gur frphevgl cebcregl nppbeqvat gb gurve cevbevgl bs yvra naq gb gur rkgrag gung shaqf erznvavat va gur unaqf bs gur fryyre ner ninvynoyr; naq ynfg, gb cnl gur erznvaqre bs fnvq cebprrqf, vs nal, gb gur iraqbe, uvf urvef, crefbanyf ercerfragngvirf, fhpprffbef be nffvtaf hcba gur qryvirel naq fheeraqre gb gur iraqrr bs cbffrffvba bs gur ynaq naq cerzvfrf, yrff pbfgf naq rkprff bs bognvavat cbffrffvba\"\n",
    "# key=\"nopqrstuvwxyzabcdefghijklm\"\n",
    "key = string.ascii_lowercase\n",
    "freq_mat = generate_freq_matrix(ciphertext,key)\n",
    "c_score=get_score(freq_mat,twogram)\n",
    "while c_score>=91.4393:\n",
    "    a,b=np.random.choice(list(string.ascii_lowercase),size=2)\n",
    "    putative_key = swap_chars(key,a,b)\n",
    "    putative_freq_mat = generate_freq_matrix(ciphertext,putative_key)\n",
    "    putative_score = get_score(putative_freq_mat,twogram)\n",
    "    #print(putative_key)\n",
    "    if putative_score < c_score:\n",
    "        key=putative_key\n",
    "        c_score = putative_score\n",
    "        print(key)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "id": "1d75950d",
   "metadata": {},
   "outputs": [],
   "source": [
    "##slow version\n",
    "import string\n",
    "import numpy as np\n",
    "from sympy import Matrix, init_printing\n",
    "init_printing()\n",
    "ciphertext=\"nn va gur rirag gung gur chepunfre qrsnhygf va gur cnlzrag bs nal vafgnyzrag bs chepunfr cevpr, gnkrf, vafhenapr, vagrerfg, be gur naahny punetr qrfpevorq ryfrjurer urerva, be funyy qrsnhyg va gur cresbeznapr bs nal bgure boyvtngvbaf frg sbegu va guvf pbagenpg, gur fryyre znl: ng uvf bcgvba: (n) qrpyner vzzrqvngryl qhr naq cnlnoyr gur ragver hacnvq onynapr bs chepunfr cevpr, jvgu nppehrq vagrerfg, gnkrf, naq naahny punetr, naq qrznaq shyy cnlzrag gurerbs, naq rasbepr pbairlnapr bs gur ynaq ol grezvangvba bs gur pbagenpg be nppbeqvat gb gur grezf urerbs, va juvpu pnfr gur chepunfre funyy nyfb or yvnoyr gb gur fryyre sbe ernfbanoyr nggbearl'f srrf sbe freivprf eraqrerq ol nal nggbearl ba orunys bs gur fryyre, be (o) fryy fnvq ynaq naq cerzvfrf be nal cneg gurerbs ng choyvp nhpgvba, va fhpu znaare, ng fhpu gvzr naq cynpr, hcba fhpu grezf naq pbaqvgvbaf, naq hcba fhpu choyvp abgvpr nf gur fryyre znl qrrz orfg sbe gur vagrerfg bs nyy pbaprearq, pbafvfgvat bs nqiregvfrzrag va n arjfcncre bs trareny pvephyngvba va gur pbhagl be pvgl va juvpu gur frphevgl cebcregl vf ybpngrq ng yrnfg bapr n jrrx sbe guerr (3) fhpprffvir jrrxf be sbe fhpu crevbq nf nccyvpnoyr ynj znl erdhver naq, va pnfr bs qrsnhyg bs nal chepunfre, gb er-fryy jvgu fhpu cbfgcbarzrag bs fnyr be erfnyr naq hcba fhpu choyvp abgvpr gurerbs nf gur fryyre znl qrgrezvar, naq hcba pbzcyvnapr ol gur chepunfre jvgu gur grezf bs fnyr, naq hcba whqvpvny nccebiny nf znl or erdhverq ol ynj, pbairl fnvq ynaq naq cerzvfrf va srr fvzcyr gb naq ng gur pbfg bs gur chepunfre, jub funyy abg or yvnoyr gb frr gb gur nccyvpngvba bs gur chepunfr zbarl; naq sebz gur cebprrqf bs gur fnyr: svefg gb cnl nyy cebcre pbfgf naq punetrf, vapyhqvat ohg abg yvzvgrq gb pbheg pbfgf, nqiregvfvat rkcrafrf, nhpgvbarre'f nyybjnapr, gur rkcrafrf, vs nal erdhverq gb pbeerpg nal veerthynevgl va gur gvgyr, cerzvhz sbe fryyre'f obaq, nhqvgbe'f srr, nggbearl'f srr, naq nyy bgure rkcrafrf bs fnyr bppheerq va naq nobhg gur cebgrpgvba naq rkrphgvba bs guvf pbagenpg, naq nyy zbarlf nqinaprq sbe gnkrf, nffrffzragf, vafhenapr, naq jvgu vagrerfg gurerba nf cebivqrq urerva, naq nyy gnkrf qhr hcba fnvq ynaq naq cerzvfrf ng gvzr bs fnyr, naq gb ergnva nf pbzcrafngvba n pbzzvffvba bs svir creprag (5%) ba gur nzbhag bs fnvq fnyr be fnyrf; frpbaq, gb cnl gur jubyr nzbhag gura erznvavat hacnvq bs gur cevapvcny bs fnvq pbagenpg, naq vagrerfg gurerba gb qngr bs cnlzrag, jurgure gur fnzr funyy or qhr be abg, vg orvat haqrefgbbq naq nterrq gung hcba fhpu fnyr orsber znghevgl bs gur pbagenpg gur onynapr gurerbs funyy or vzzrqvngryl qhr naq cnlnoyr; guveq, gb cnl yvraf bs erpbeq ntnvafg gur frphevgl cebcregl nppbeqvat gb gurve cevbevgl bs yvra naq gb gur rkgrag gung shaqf erznvavat va gur unaqf bs gur fryyre ner ninvynoyr; naq ynfg, gb cnl gur erznvaqre bs fnvq cebprrqf, vs nal, gb gur iraqbe, uvf urvef, crefbanyf ercerfragngvirf, fhpprffbef be nffvtaf hcba gur qryvirel naq fheeraqre gb gur iraqrr bs cbffrffvba bs gur ynaq naq cerzvfrf, yrff pbfgf naq rkprff bs bognvavat cbffrffvba\"\n",
    "# key=\"nopqrstuvwxyzabcdefghijklm\"\n",
    "key = string.ascii_lowercase\n",
    "freq_mat = generate_freq_matrix(ciphertext,key)\n",
    "c_score=get_score(freq_mat,twogram)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "id": "9365c380",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$\\displaystyle \\left[\\begin{array}{cccccccccccccccccccccccccc}0.00156657963446475 & 0.022976501305483 & 0.0 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.00365535248041775 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0308093994778068 & 0.0 & 0.0 & 0.0 & 0.0130548302872063 & 0.0 & 0.000522193211488251 & 0.0 & 0.0224543080939948 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.00261096605744125 & 0.000522193211488251 & 0.00626631853785901 & 0.0 & 0.00469973890339426 & 0.002088772845953 & 0.0120104438642298 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.012532637075718 & 0.000522193211488251 & 0.00365535248041775 & 0.00574412532637076 & 0.0 & 0.0010443864229765 & 0.00835509138381201 & 0.0 & 0.0 & 0.0010443864229765 & 0.002088772845953\\\\0.0010443864229765 & 0.002088772845953 & 0.00156657963446475 & 0.0 & 0.0 & 0.000522193211488251 & 0.000522193211488251 & 0.00417754569190601 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.00156657963446475\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0172323759791123 & 0.00939947780678851 & 0.0 & 0.002088772845953 & 0.0 & 0.00261096605744125 & 0.00835509138381201 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00365535248041775 & 0.0 & 0.0010443864229765 & 0.0 & 0.0287206266318538 & 0.000522193211488251 & 0.000522193211488251 & 0.000522193211488251 & 0.00522193211488251 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.00678851174934726 & 0.00365535248041775 & 0.0 & 0.0 & 0.00261096605744125 & 0.00626631853785901 & 0.00261096605744125 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0109660574412533 & 0.0 & 0.0 & 0.002088772845953 & 0.016710182767624 & 0.0 & 0.0 & 0.0 & 0.00678851174934726 & 0.000522193211488251 & 0.000522193211488251 & 0.00156657963446475 & 0.00156657963446475\\\\0.0130548302872063 & 0.00417754569190601 & 0.000522193211488251 & 0.0 & 0.00365535248041775 & 0.00939947780678851 & 0.00156657963446475 & 0.00156657963446475 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0120104438642298 & 0.000522193211488251 & 0.00469973890339426 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.00783289817232376 & 0.0 & 0.0 & 0.00156657963446475 & 0.0\\\\0.0010443864229765 & 0.00261096605744125 & 0.00574412532637076 & 0.00156657963446475 & 0.000522193211488251 & 0.00678851174934726 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0031331592689295 & 0.000522193211488251 & 0.00261096605744125 & 0.002088772845953 & 0.0 & 0.0010443864229765 & 0.000522193211488251 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.000522193211488251 & 0.0\\\\0.0010443864229765 & 0.0010443864229765 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.00156657963446475 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.0031331592689295 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.00156657963446475 & 0.0 & 0.00731070496083551 & 0.0 & 0.00417754569190601 & 0.00887728459530026 & 0.00417754569190601 & 0.0010443864229765 & 0.00156657963446475 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0010443864229765 & 0.00261096605744125 & 0.000522193211488251 & 0.0010443864229765 & 0.00156657963446475 & 0.0010443864229765 & 0.0109660574412533 & 0.0031331592689295 & 0.00156657963446475 & 0.0 & 0.00783289817232376 & 0.00626631853785901\\\\0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00417754569190601 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.00678851174934726 & 0.002088772845953 & 0.0 & 0.0 & 0.00574412532637076 & 0.000522193211488251 & 0.0 & 0.00626631853785901 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.00469973890339426 & 0.0 & 0.0031331592689295 & 0.0 & 0.00417754569190601 & 0.0 & 0.0 & 0.0 & 0.00678851174934726 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0255874673629243 & 0.0010443864229765 & 0.0 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.00835509138381201 & 0.0 & 0.0 & 0.0 & 0.00469973890339426 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.00626631853785901 & 0.0 & 0.00626631853785901 & 0.0 & 0.025065274151436 & 0.0198433420365535 & 0.00887728459530026 & 0.00261096605744125 & 0.00574412532637076 & 0.0010443864229765 & 0.00261096605744125 & 0.0 & 0.0 & 0.0 & 0.00522193211488251 & 0.012532637075718 & 0.00783289817232376 & 0.00731070496083551 & 0.002088772845953 & 0.002088772845953 & 0.0318537859007833 & 0.0010443864229765 & 0.0 & 0.0 & 0.014621409921671 & 0.00626631853785901\\\\0.000522193211488251 & 0.0198433420365535 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.000522193211488251 & 0.0\\\\0.00469973890339426 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00261096605744125 & 0.0349869451697128 & 0.0 & 0.0 & 0.0031331592689295 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0109660574412533 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.00574412532637076 & 0.00417754569190601 & 0.0109660574412533 & 0.00156657963446475 & 0.0010443864229765 & 0.002088772845953 & 0.0 & 0.0 & 0.0 & 0.00783289817232376 & 0.0 & 0.002088772845953 & 0.00417754569190601 & 0.00261096605744125 & 0.0010443864229765 & 0.0 & 0.00365535248041775 & 0.0 & 0.0 & 0.0 & 0.00678851174934726 & 0.00469973890339426\\\\0.00417754569190601 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.00469973890339426 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00731070496083551 & 0.002088772845953 & 0.0 & 0.0 & 0.00365535248041775 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0010443864229765 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & 0.000522193211488251 & 0.00261096605744125 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0031331592689295 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.016710182767624 & 0.00574412532637076 & 0.0010443864229765 & 0.0 & 0.00678851174934726 & 0.0 & 0.0 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0109660574412533 & 0.0\\\\0.0 & 0.002088772845953 & 0.0 & 0.0 & 0.0031331592689295 & 0.000522193211488251 & 0.0 & 0.000522193211488251 & 0.0 & 0.0 & 0.0 & 0.0 & 0.0 & 0.00156657963446475 & 0.0 & 0.0 & 0.0 & 0.00626631853785901 & 0.0 & 0.0 & 0.0 & 0.0031331592689295 & 0.00156657963446475 & 0.0 & 0.000522193211488251 & 0.00156657963446475\\end{array}\\right]$"
      ],
      "text/plain": [
       "⎡0.00156657963446475    0.022976501305483            0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢0.00261096605744125   0.000522193211488251  0.00626631853785901           0.0\n",
       "⎢                                                                             \n",
       "⎢ 0.0010443864229765    0.002088772845953    0.00156657963446475           0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0            0.0172323759791123   0.00939947780678851           0.0\n",
       "⎢                                                                             \n",
       "⎢0.00678851174934726   0.00365535248041775           0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢ 0.0130548302872063   0.00417754569190601   0.000522193211488251          0.0\n",
       "⎢                                                                             \n",
       "⎢ 0.0010443864229765   0.00261096605744125   0.00574412532637076   0.001566579\n",
       "⎢                                                                             \n",
       "⎢ 0.0010443864229765    0.0010443864229765           0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0           0.000522193211488251          0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢0.00156657963446475           0.0           0.00731070496083551           0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0            0.0010443864229765           0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢0.00678851174934726    0.002088772845953            0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢ 0.0255874673629243    0.0010443864229765           0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢0.00626631853785901           0.0           0.00626631853785901           0.0\n",
       "⎢                                                                             \n",
       "⎢0.000522193211488251   0.0198433420365535           0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢0.00469973890339426           0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢0.00156657963446475           0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢0.00417754569190601           0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0                   0.0                   0.0                   0.0\n",
       "⎢                                                                             \n",
       "⎢        0.0           0.000522193211488251  0.00261096605744125           0.0\n",
       "⎢                                                                             \n",
       "⎣        0.0            0.002088772845953            0.0                   0.0\n",
       "\n",
       "           0.002088772845953            0.0                   0.0           0.\n",
       "                                                                              \n",
       "          0.00469973890339426    0.002088772845953     0.0120104438642298     \n",
       "                                                                              \n",
       "                  0.0           0.000522193211488251  0.000522193211488251  0.\n",
       "                                                                              \n",
       "                  0.0                   0.0                   0.0             \n",
       "                                                                              \n",
       "           0.002088772845953            0.0           0.00261096605744125   0.\n",
       "                                                                              \n",
       "          0.00261096605744125   0.00626631853785901   0.00261096605744125     \n",
       "                                                                              \n",
       "          0.00365535248041775   0.00939947780678851   0.00156657963446475   0.\n",
       "                                                                              \n",
       "63446475  0.000522193211488251  0.00678851174934726   0.000522193211488251    \n",
       "                                                                              \n",
       "          0.000522193211488251          0.0                   0.0             \n",
       "                                                                              \n",
       "                  0.0                   0.0                   0.0             \n",
       "                                                                              \n",
       "                  0.0                   0.0                   0.0             \n",
       "                                                                              \n",
       "                  0.0                   0.0                   0.0             \n",
       "                                                                              \n",
       "                  0.0                   0.0                   0.0             \n",
       "                                                                              \n",
       "          0.00417754569190601   0.00887728459530026   0.00417754569190601    0\n",
       "                                                                              \n",
       "                  0.0                   0.0                   0.0           0.\n",
       "                                                                              \n",
       "          0.00574412532637076   0.000522193211488251          0.0           0.\n",
       "                                                                              \n",
       "           0.002088772845953            0.0                   0.0           0.\n",
       "                                                                              \n",
       "           0.025065274151436     0.0198433420365535   0.00887728459530026   0.\n",
       "                                                                              \n",
       "          0.000522193211488251          0.0                   0.0             \n",
       "                                                                              \n",
       "          0.00156657963446475           0.0                   0.0             \n",
       "                                                                              \n",
       "                  0.0           0.00261096605744125    0.0349869451697128     \n",
       "                                                                              \n",
       "          0.00574412532637076   0.00417754569190601    0.0109660574412533   0.\n",
       "                                                                              \n",
       "          0.000522193211488251          0.0           0.00469973890339426     \n",
       "                                                                              \n",
       "                  0.0                   0.0                   0.0             \n",
       "                                                                              \n",
       "                  0.0                   0.0           0.000522193211488251   0\n",
       "                                                                              \n",
       "           0.0031331592689295   0.000522193211488251          0.0           0.\n",
       "\n",
       "00365535248041775           0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "00417754569190601           0.0                  0.0           0.0015665796344\n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "00835509138381201           0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0          0.000522193211488251          0.0    \n",
       "                                                                              \n",
       "00156657963446475           0.0                  0.0           0.0005221932114\n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       ".0010443864229765   0.00156657963446475  0.000522193211488251          0.0    \n",
       "                                                                              \n",
       "00156657963446475           0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "00626631853785901           0.0                  0.0           0.0005221932114\n",
       "                                                                              \n",
       "00156657963446475           0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "00261096605744125   0.00574412532637076   0.0010443864229765   0.0026109660574\n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0           0.0031331592689295           0.0    \n",
       "                                                                              \n",
       "00156657963446475   0.0010443864229765    0.002088772845953            0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "      0.0                   0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       ".0031331592689295           0.0                  0.0                   0.0    \n",
       "                                                                              \n",
       "000522193211488251          0.0                  0.0                   0.0    \n",
       "\n",
       "               0.0           0.0   0.0308093994778068           0.0           \n",
       "                                                                              \n",
       "               0.0           0.0          0.0            0.0010443864229765   \n",
       "                                                                              \n",
       "6475           0.0           0.0   0.002088772845953            0.0           \n",
       "                                                                              \n",
       "               0.0           0.0          0.0                   0.0           \n",
       "                                                                              \n",
       "               0.0           0.0  0.00365535248041775           0.0           \n",
       "                                                                              \n",
       "               0.0           0.0   0.0109660574412533           0.0           \n",
       "                                                                              \n",
       "88251          0.0           0.0   0.0120104438642298   0.000522193211488251  \n",
       "                                                                              \n",
       "       0.000522193211488251  0.0   0.0031331592689295   0.000522193211488251  \n",
       "                                                                              \n",
       "               0.0           0.0  0.000522193211488251          0.0           \n",
       "                                                                              \n",
       "               0.0           0.0   0.0010443864229765           0.0           \n",
       "                                                                              \n",
       "               0.0           0.0   0.002088772845953            0.0           \n",
       "                                                                              \n",
       "               0.0           0.0          0.0                   0.0           \n",
       "                                                                              \n",
       "               0.0           0.0          0.0                   0.0           \n",
       "                                                                              \n",
       "               0.0           0.0  0.000522193211488251   0.0010443864229765   \n",
       "                                                                              \n",
       "               0.0           0.0  0.00417754569190601           0.0           \n",
       "                                                                              \n",
       "88251          0.0           0.0  0.00469973890339426           0.0           \n",
       "                                                                              \n",
       "               0.0           0.0  0.00156657963446475           0.0           \n",
       "                                                                              \n",
       "4125           0.0           0.0          0.0           0.00522193211488251   \n",
       "                                                                              \n",
       "               0.0           0.0          0.0                   0.0           \n",
       "                                                                              \n",
       "               0.0           0.0   0.0010443864229765           0.0           \n",
       "                                                                              \n",
       "               0.0           0.0          0.0                   0.0           \n",
       "                                                                              \n",
       "               0.0           0.0  0.00783289817232376           0.0           \n",
       "                                                                              \n",
       "               0.0           0.0  0.00731070496083551    0.002088772845953    \n",
       "                                                                              \n",
       "               0.0           0.0          0.0                   0.0           \n",
       "                                                                              \n",
       "               0.0           0.0   0.016710182767624    0.00574412532637076   \n",
       "                                                                              \n",
       "               0.0           0.0  0.00156657963446475           0.0           \n",
       "\n",
       "        0.0                  0.0            0.0130548302872063           0.0  \n",
       "                                                                              \n",
       " 0.012532637075718   0.000522193211488251  0.00365535248041775   0.00574412532\n",
       "                                                                              \n",
       "        0.0                  0.0           0.000522193211488251          0.0  \n",
       "                                                                              \n",
       "        0.0                  0.0           0.00156657963446475           0.0  \n",
       "                                                                              \n",
       "0.0010443864229765           0.0            0.0287206266318538   0.00052219321\n",
       "                                                                              \n",
       "        0.0           0.002088772845953     0.016710182767624            0.0  \n",
       "                                                                              \n",
       "0.00469973890339426          0.0            0.002088772845953            0.0  \n",
       "                                                                              \n",
       "0.00261096605744125   0.002088772845953            0.0            0.0010443864\n",
       "                                                                              \n",
       "        0.0          0.00156657963446475   0.000522193211488251          0.0  \n",
       "                                                                              \n",
       "        0.0                  0.0            0.0010443864229765           0.0  \n",
       "                                                                              \n",
       "        0.0                  0.0            0.0031331592689295           0.0  \n",
       "                                                                              \n",
       "        0.0                  0.0                   0.0                   0.0  \n",
       "                                                                              \n",
       "        0.0                  0.0                   0.0                   0.0  \n",
       "                                                                              \n",
       "0.00261096605744125  0.000522193211488251   0.0010443864229765   0.00156657963\n",
       "                                                                              \n",
       "        0.0                  0.0                   0.0                   0.0  \n",
       "                                                                              \n",
       "0.0031331592689295           0.0           0.00417754569190601           0.0  \n",
       "                                                                              \n",
       "        0.0                  0.0           0.00835509138381201           0.0  \n",
       "                                                                              \n",
       " 0.012532637075718   0.00783289817232376   0.00731070496083551    0.0020887728\n",
       "                                                                              \n",
       "        0.0                  0.0            0.002088772845953            0.0  \n",
       "                                                                              \n",
       "        0.0                  0.0           0.000522193211488251          0.0  \n",
       "                                                                              \n",
       "0.0109660574412533           0.0           0.000522193211488251          0.0  \n",
       "                                                                              \n",
       " 0.002088772845953   0.00417754569190601   0.00261096605744125    0.0010443864\n",
       "                                                                              \n",
       "        0.0                  0.0           0.00365535248041775           0.0  \n",
       "                                                                              \n",
       "        0.0                  0.0            0.0010443864229765           0.0  \n",
       "                                                                              \n",
       "0.0010443864229765           0.0           0.00678851174934726           0.0  \n",
       "                                                                              \n",
       "        0.0                  0.0           0.00626631853785901           0.0  \n",
       "\n",
       "         0.000522193211488251          0.0            0.0224543080939948      \n",
       "                                                                              \n",
       "637076           0.0            0.0010443864229765   0.00835509138381201      \n",
       "                                                                              \n",
       "                 0.0                   0.0           0.000522193211488251     \n",
       "                                                                              \n",
       "                 0.0                   0.0                   0.0              \n",
       "                                                                              \n",
       "1488251  0.000522193211488251  0.000522193211488251  0.00522193211488251      \n",
       "                                                                              \n",
       "                 0.0                   0.0           0.00678851174934726   0.0\n",
       "                                                                              \n",
       "                 0.0                   0.0           0.00783289817232376      \n",
       "                                                                              \n",
       "229765   0.000522193211488251          0.0           0.000522193211488251     \n",
       "                                                                              \n",
       "                 0.0                   0.0            0.002088772845953       \n",
       "                                                                              \n",
       "                 0.0                   0.0                   0.0              \n",
       "                                                                              \n",
       "                 0.0                   0.0                   0.0              \n",
       "                                                                              \n",
       "                 0.0                   0.0                   0.0              \n",
       "                                                                              \n",
       "                 0.0                   0.0                   0.0              \n",
       "                                                                              \n",
       "446475    0.0010443864229765    0.0109660574412533    0.0031331592689295   0.0\n",
       "                                                                              \n",
       "                 0.0                   0.0           0.000522193211488251     \n",
       "                                                                              \n",
       "                 0.0                   0.0           0.00678851174934726      \n",
       "                                                                              \n",
       "                 0.0                   0.0           0.00469973890339426      \n",
       "                                                                              \n",
       "45953     0.002088772845953     0.0318537859007833    0.0010443864229765      \n",
       "                                                                              \n",
       "                 0.0                   0.0            0.0010443864229765      \n",
       "                                                                              \n",
       "                 0.0                   0.0            0.0010443864229765      \n",
       "                                                                              \n",
       "                 0.0                   0.0                   0.0              \n",
       "                                                                              \n",
       "229765           0.0           0.00365535248041775           0.0              \n",
       "                                                                              \n",
       "                 0.0                   0.0                   0.0              \n",
       "                                                                              \n",
       "                 0.0                   0.0                   0.0              \n",
       "                                                                              \n",
       "                 0.0                   0.0           0.000522193211488251     \n",
       "                                                                              \n",
       "                 0.0                   0.0            0.0031331592689295   0.0\n",
       "\n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0            0.0010443864229765    0.002088772845\n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0           0.0015665796344\n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "00522193211488251  0.000522193211488251  0.00156657963446475   0.0015665796344\n",
       "                                                                              \n",
       "     0.0                   0.0           0.00156657963446475           0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0           0.000522193211488251          0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "0156657963446475           0.0           0.00783289817232376   0.0062663185378\n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0            0.014621409921671    0.0062663185378\n",
       "                                                                              \n",
       "     0.0                   0.0           0.000522193211488251          0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0           0.00678851174934726   0.0046997389033\n",
       "                                                                              \n",
       "     0.0                   0.0            0.0010443864229765           0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0                   0.0                   0.0    \n",
       "                                                                              \n",
       "     0.0                   0.0            0.0109660574412533           0.0    \n",
       "                                                                              \n",
       "0156657963446475           0.0           0.000522193211488251  0.0015665796344\n",
       "\n",
       "    ⎤\n",
       "    ⎥\n",
       "953 ⎥\n",
       "    ⎥\n",
       "6475⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "6475⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "5901⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "5901⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "9426⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "    ⎥\n",
       "6475⎦"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(Matrix(generate_freq_matrix(ciphertext,\"abcdefghijkwmnopqrstuvlxyz\")))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 85,
   "id": "3e6fa244",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[0.00156658, 0.0229765 , 0.        , 0.        , 0.00208877,\n",
       "        0.        , 0.        , 0.00365535, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.0308094 , 0.        ,\n",
       "        0.        , 0.        , 0.01305483, 0.        , 0.00052219,\n",
       "        0.        , 0.02245431, 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.00261097, 0.00052219, 0.00626632, 0.        , 0.00469974,\n",
       "        0.00208877, 0.01201044, 0.        , 0.        , 0.        ,\n",
       "        0.00208877, 0.        , 0.        , 0.        , 0.00104439,\n",
       "        0.01253264, 0.00052219, 0.00365535, 0.00574413, 0.        ,\n",
       "        0.00104439, 0.00835509, 0.        , 0.        , 0.00104439,\n",
       "        0.        ],\n",
       "       [0.00104439, 0.00208877, 0.00156658, 0.        , 0.        ,\n",
       "        0.00052219, 0.00052219, 0.00417755, 0.        , 0.        ,\n",
       "        0.00156658, 0.        , 0.        , 0.00208877, 0.        ,\n",
       "        0.        , 0.        , 0.00052219, 0.        , 0.        ,\n",
       "        0.        , 0.00052219, 0.        , 0.        , 0.        ,\n",
       "        0.00156658],\n",
       "       [0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.00156658, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.        , 0.01723238, 0.00939948, 0.        , 0.00208877,\n",
       "        0.        , 0.00261097, 0.00835509, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.00365535, 0.        ,\n",
       "        0.00104439, 0.        , 0.02872063, 0.00052219, 0.00052219,\n",
       "        0.00052219, 0.00522193, 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.00678851, 0.00365535, 0.        , 0.        , 0.00261097,\n",
       "        0.00626632, 0.00261097, 0.        , 0.        , 0.00052219,\n",
       "        0.00156658, 0.00052219, 0.        , 0.01096606, 0.        ,\n",
       "        0.        , 0.00208877, 0.01671018, 0.        , 0.        ,\n",
       "        0.        , 0.00678851, 0.        , 0.00052219, 0.00156658,\n",
       "        0.        ],\n",
       "       [0.01305483, 0.00417755, 0.00052219, 0.        , 0.00365535,\n",
       "        0.00939948, 0.00156658, 0.00156658, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.01201044, 0.00052219,\n",
       "        0.00469974, 0.        , 0.00208877, 0.        , 0.        ,\n",
       "        0.        , 0.0078329 , 0.        , 0.        , 0.00156658,\n",
       "        0.00052219],\n",
       "       [0.00104439, 0.00261097, 0.00574413, 0.00156658, 0.00052219,\n",
       "        0.00678851, 0.00052219, 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.00313316, 0.00052219,\n",
       "        0.00261097, 0.00208877, 0.        , 0.00104439, 0.00052219,\n",
       "        0.        , 0.00052219, 0.00052219, 0.        , 0.00052219,\n",
       "        0.        ],\n",
       "       [0.00104439, 0.00104439, 0.        , 0.        , 0.00052219,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.00052219, 0.        ,\n",
       "        0.        , 0.00156658, 0.00052219, 0.        , 0.        ,\n",
       "        0.        , 0.00208877, 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.        , 0.00052219, 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.00104439, 0.        ,\n",
       "        0.        , 0.        , 0.00104439, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.        , 0.00208877, 0.        , 0.        , 0.00313316,\n",
       "        0.00052219, 0.        , 0.00052219, 0.        , 0.        ,\n",
       "        0.00156658, 0.00156658, 0.        , 0.00156658, 0.        ,\n",
       "        0.        , 0.        , 0.00626632, 0.        , 0.        ,\n",
       "        0.        , 0.00313316, 0.        , 0.        , 0.00052219,\n",
       "        0.        ],\n",
       "       [0.00417755, 0.        , 0.        , 0.        , 0.00052219,\n",
       "        0.        , 0.00469974, 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.0073107 , 0.00208877,\n",
       "        0.        , 0.        , 0.00365535, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.00104439,\n",
       "        0.        ],\n",
       "       [0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.00156658, 0.        , 0.0073107 , 0.        , 0.00417755,\n",
       "        0.00887728, 0.00417755, 0.00104439, 0.00156658, 0.00052219,\n",
       "        0.00626632, 0.00156658, 0.        , 0.00052219, 0.00104439,\n",
       "        0.00261097, 0.00052219, 0.00104439, 0.00156658, 0.00104439,\n",
       "        0.01096606, 0.00313316, 0.        , 0.        , 0.0078329 ,\n",
       "        0.        ],\n",
       "       [0.        , 0.00104439, 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.00156658, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.00417755, 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.00052219, 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.00678851, 0.00208877, 0.        , 0.        , 0.00574413,\n",
       "        0.00052219, 0.        , 0.00626632, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.00469974, 0.        ,\n",
       "        0.00313316, 0.        , 0.00417755, 0.        , 0.        ,\n",
       "        0.        , 0.00678851, 0.        , 0.        , 0.        ,\n",
       "        0.00052219],\n",
       "       [0.02558747, 0.00104439, 0.        , 0.        , 0.00208877,\n",
       "        0.        , 0.        , 0.00156658, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.00156658, 0.        ,\n",
       "        0.        , 0.        , 0.00835509, 0.        , 0.        ,\n",
       "        0.        , 0.00469974, 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.00626632, 0.        , 0.00626632, 0.        , 0.02506527,\n",
       "        0.01984334, 0.00887728, 0.00261097, 0.00574413, 0.00104439,\n",
       "        0.00626632, 0.        , 0.        , 0.        , 0.00522193,\n",
       "        0.01253264, 0.0078329 , 0.0073107 , 0.00208877, 0.00208877,\n",
       "        0.03185379, 0.00104439, 0.        , 0.        , 0.01462141,\n",
       "        0.00261097],\n",
       "       [0.00052219, 0.01984334, 0.        , 0.        , 0.00052219,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.00208877, 0.        , 0.        ,\n",
       "        0.        , 0.00104439, 0.        , 0.        , 0.00052219,\n",
       "        0.        ],\n",
       "       [0.00469974, 0.        , 0.        , 0.        , 0.00156658,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.00104439, 0.        ,\n",
       "        0.        , 0.        , 0.00052219, 0.        , 0.        ,\n",
       "        0.        , 0.00104439, 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.00261097, 0.03498695, 0.        , 0.        , 0.00313316,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.01096606, 0.        , 0.00052219, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.00156658, 0.        , 0.        , 0.        , 0.00574413,\n",
       "        0.00417755, 0.01096606, 0.00156658, 0.00104439, 0.00208877,\n",
       "        0.00469974, 0.        , 0.        , 0.0078329 , 0.        ,\n",
       "        0.00208877, 0.00417755, 0.00261097, 0.00104439, 0.        ,\n",
       "        0.00365535, 0.        , 0.        , 0.        , 0.00678851,\n",
       "        0.        ],\n",
       "       [0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.00104439, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        ],\n",
       "       [0.        , 0.00052219, 0.00261097, 0.        , 0.        ,\n",
       "        0.        , 0.00052219, 0.00313316, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.01671018, 0.00574413,\n",
       "        0.00104439, 0.        , 0.00678851, 0.        , 0.        ,\n",
       "        0.        , 0.00052219, 0.        , 0.        , 0.01096606,\n",
       "        0.        ],\n",
       "       [0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.00208877, 0.        ,\n",
       "        0.        , 0.        , 0.00313316, 0.        , 0.        ,\n",
       "        0.        , 0.        , 0.        , 0.        , 0.        ,\n",
       "        0.        ]])"
      ]
     },
     "execution_count": 85,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "generate_freq_matrix(ciphertext,swap_chars(key,\"z\",\"k\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 87,
   "id": "8b42a14d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 87,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.all( fswap(generate_freq_matrix(ciphertext,key),\"z\",\"k\")==generate_freq_matrix(ciphertext,swap_chars(key,\"z\",\"k\")) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "284f06b1",
   "metadata": {},
   "outputs": [],
   "source": [
    "key = string.ascii_lowercase\n",
    "freq_mat = generate_freq_matrix(ciphertext,key)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "id": "0682004a",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "56.1 µs ± 1.95 µs per loop (mean ± std. dev. of 7 runs, 10000 loops each)\n"
     ]
    }
   ],
   "source": [
    "%%timeit\n",
    "c_score=get_score(freq_mat,twogram)\n",
    "a,b=np.random.choice(list(string.ascii_lowercase),size=2)\n",
    "putative_key = swap_chars(key,a,b)\n",
    "# putative_freq_mat = fswap(freq_mat,a,b)\n",
    "# putative_score = get_score(putative_freq_mat,twogram)\n",
    "# if putative_score < c_score:\n",
    "#     key=putative_key\n",
    "#     c_score = putative_score\n",
    "#     print(key)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7d6adf2c",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "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.9.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}