From c22feadcd68a2702846bc5d9542f171b64d0ecd2 Mon Sep 17 00:00:00 2001
From: "Ira W. Snyder" <devel@irasnyder.com>
Date: Fri, 19 Oct 2007 10:43:52 -0700
Subject: [PATCH] Initial commit

Signed-off-by: Ira W. Snyder <devel@irasnyder.com>
---
 Makefile     |  29 +++
 Problem2.in  |  55 ++++++
 Problem3.out | Bin 0 -> 13256 bytes
 aes.cpp      | 523 +++++++++++++++++++++++++++++++++++++++++++++++++++
 aes.hpp      |  56 ++++++
 aes_test.cpp | 146 ++++++++++++++
 problem2.cpp |  41 ++++
 problem3.cpp |  40 ++++
 8 files changed, 890 insertions(+)
 create mode 100644 Makefile
 create mode 100644 Problem2.in
 create mode 100644 Problem3.out
 create mode 100644 aes.cpp
 create mode 100644 aes.hpp
 create mode 100644 aes_test.cpp
 create mode 100644 problem2.cpp
 create mode 100644 problem3.cpp

diff --git a/Makefile b/Makefile
new file mode 100644
index 0000000..7e77240
--- /dev/null
+++ b/Makefile
@@ -0,0 +1,29 @@
+CC=g++
+CFLAGS=-march=athlon-xp -O2
+
+all: problem2 problem3
+
+aes_test: aes.o aes_test.o
+	$(CC) $(CFLAGS) -o $@ $^
+
+problem2: aes.o problem2.o
+	$(CC) $(CFLAGS) -o $@ $^
+
+problem3: aes.o problem3.o
+	$(CC) $(CFLAGS) -o $@ $^
+
+clean:
+	rm -f *.o aes_test problem2 problem3
+
+.cpp.o:
+	$(CC) $(CFLAGS) -c -o $@ $<
+
+run: run_p2 run_p3
+
+run_p2: problem2
+	@cat Problem2.in | ./$<
+
+run_p3: problem3
+	@./$<
+
+.PHONY: run all clean
diff --git a/Problem2.in b/Problem2.in
new file mode 100644
index 0000000..5afe35d
--- /dev/null
+++ b/Problem2.in
@@ -0,0 +1,55 @@
+#PN*3kUAprThcbB]
+                                      .M
+                                 .:AMMO:
+                        .:AMMMMMHIIIHMMM.
+             ....   .AMMMMMMMMMMMHHHMHHMMMML:AMF"
+             .:MMMMMLAMMMMMMMHMMMMMMHHIHHIIIHMMMML.
+                 "WMMMMMMMMMMMMMMMMMMH:::::HMMMMMMHII:.
+            .AMMMMMMMHHHMMMMMMMMMMHHHHHMMMMMMMMMAMMMHHHHL.
+          .MMMMMMMMMMHHMMMMMMMMHHHHMMMMMMMMMMMMMHTWMHHHHHML
+         .MMMMMMMMMMMMMMMMMMMHHHHHHHHHMHMMHHHHIII:::HMHHHHMM.
+         .MMMMMMMMMMMMMMMMMMMMMMHHHHHHMHHHHHHIIIIIIIIHMHHHHHM.
+         MMMMMMMMMMMMMMMMMHHMMHHHHHIIIHHH::IIHHII:::::IHHHHHHHL
+         "MMMMMMMMMMMMMMMMHIIIHMMMMHHIIHHLI::IIHHHHIIIHHHHHHHHML
+         .MMMMMMMMMMMMMM"WMMMHHHMMMMMMMMMMMLHHHMMMMMMHHHHHHHHHHH
+        .MMMMMMMMMMMWWMW""YYHMMMMMMMMMMMMF""HMMMMMMMMMHHHHHHHH.
+        .MMMMMMMMMM W" V                         W"WMMMMMHHHHHHHHHH
+       "MMMMMMMMMM".                                 "WHHHMH"HHHHHHL
+       MMMMMMMMMMF  .                                         IHHHHH.
+       MMMMMMMMMM .                                  .        HHHHHHH
+       MMMMMMMMMF. .                               .  .       HHHHHHH.
+       MMMMMMMMM .     ,AWMMMMML.              ..    .  .     HHHHHHH.
+     :MMMMMMMMM".  .  F"'    'WM:.         ,::HMMA, .  .      HHHHMMM
+     :MMMMMMMMF.  . ."         WH..      AMM"'     "  .  .    HHHMMMM
+      MMMMMMMM . .     ,;AAAHHWL"..     .:'                   HHHHHHH
+      MMMMMMM:. . .   -MK"OTO L :I..    ...:HMA-.             "HHHHHH
+ ,:IIIILTMMMMI::.      L,,,,.  ::I..    .. K"OTO"ML           'HHHHHH
+ LHT::LIIIIMMI::. .      '""'.IHH:..    .. :.,,,,           '  HMMMH: HLI'
+ ILTT::"IIITMII::.  .         .IIII.     . '""""             ' MMMFT:::.
+ HML:::WMIINMHI:::.. .          .:I.     .   . .  .        '  .M"'.....I.
+ "HWHINWI:.'.HHII::..          .HHI     .II.    .  .      . . :M.',, ..I:
+  "MLI"ML': :HHII::...        MMHHL     :::::  . :..      .'.'.'HHTML.II:
+   "MMLIHHWL:IHHII::....:I:" :MHHWHI:...:W,,"  '':::.      ..'  ":.HH:II:
+     "MMMHITIIHHH:::::IWF"    """T99"'  '""    '.':II:..'.'..'  I'.HHIHI'
+       YMMHII:IHHHH:::IT..     . .   ...  . .    ''THHI::.'.' .;H.""."H"
+         HHII:MHHI"::IWWL     . .     .    .  .     HH"HHHIIHHH":HWWM"
+          """ MMHI::HY""ML,          ...     . ..  :"  :HIIIIIILTMH"
+               MMHI:.'    'HL,,,,,,,,..,,,......,:" . ''::HH "HWW
+               'MMH:..   . 'MMML,: """MM""""MMM"      .'.IH'"MH"
+                "MMHL..   .. "MMMMMML,MM,HMMMF    .   .IHM"
+                  "MMHHL    .. "MMMMMMMMMMMM"  . .  '.IHF'
+                    'MMMML    .. "MMMMMMMM"  .     .'HMF
+                     HHHMML.                    .'MMF"
+                    IHHHHHMML.               .'HMF"
+                    HHHHHHITMML.           .'IF..
+                    "HHHHHHIITML,.       ..:F...
+                     'HHHHHHHHHMMWWWWWW::"......
+                       HHHHHHHMMMMMMF"'........
+                        HHHHHHHHHH............
+                          HHHHHHHH...........
+                           HHHHIII..........
+                            HHIII..........
+                             HII.........
+                              "H........
+                                ......
+
diff --git a/Problem3.out b/Problem3.out
new file mode 100644
index 0000000000000000000000000000000000000000..276f784aa7705c3e4d764a58657efab215860fc2
GIT binary patch
literal 13256
zcmbuGQ*>or)3#$f>DacN?%1}Sbj%&wwr$&X(y`g;*mlSK|Kr^U&v?goSmUm`s%p)J
zy=JX#<3hviVJhe5N#|tZtZ2YTYo(&V1YpqU7;nG6F|NYTDE-y^1Ru1Nw?U1f0LkxL
zh{diOV*KY_Kqq6TACG$HL9~f9$zkxmh9;gsOrTX+7V>Arkq_o6u)fs&Gr^2PqBy<8
zlDzIjvj4r{&*1x+#hI(NUAmKuLFg;RBa{|ETFu;Cj$e7clcV_e0U|dExe^#}pk{>G
z?+aj_>1;qfcrUcXf{s(IQeo2gxU{RZ%i`>Sgg`SacY*X9us_URKhEvdz)8n49lOVt
zkflfvQD2Z+(BL^+#rnXV#|Ln~kWybhdWwElMja0^IchjOj<W^B+9BwUm!_16o@SaR
z-uRb{x9YcN4^TM!Aq<ko!1*2|B(=VQFCM)}d%;z(ySC7!g3Y1U7*L5Hb7dTWpX>wo
z5Bc)v_r&PRYOOHHJ<SRarU_s!H5_$F;^HiGO`3Ztqf|*Bn;D7Q`pjhP6CWO>$(b@f
z|5o^s4+o~PZ*VrxQX~a}o)~Oo^9&^K{zqesA8<WBxL}G(cm2Yb80aH{k7jF_evYQS
zC!--V;8m*ZN=gv;ARvRI<8LL1RVs=zz7<Yy{I@sTb#72F{_>|xpB$RNyqggHm(OOd
z=-<Yb)CbG}9|cVz`!aKu)N|TkZ~|o|FRI>kmD1ul_u;0RZAvtBF=)uJ&k8hnnhK$t
zS^MBH1Yq8S-JNLu!)kF)Z?jz6Jv+3iy$=gPRj&ewK<RkuOICz0H^s6bnMR78J(9~L
zH2|fI^t6otTy5mRAtn-#ToMG#lh?S(m@*>fUSPQm2n%1ZWdb-~an(`o{q>}HgB))g
z74l|2X3WoJPO#{*)kZGSNlLz?oMFVBV){t!tr*!i0Oq+pR+6Sf62IZdV%zq6Eyb3D
zfR1l;ux-F9NV5Bbyd!44_9%BQC-c`yZA0*hYe!jb*-ByJYVueqs~16fFYp8N4it7M
zynhN%ru8^4fu3hJ1M9(NW?hATH(Ea&41*LzVHTGsex_*7uJ)Bns~G!Gt;PBNLg!(1
zm57h}FbMPq1S)bU^dl2wXK5dtGm$zQWP||oQV%>vioeZee}0kS5g;(P0H5!`{+(in
zx;7MYM$>eAFYGrSMF?@$m*eKY^+qU|oOa5%iH9MLd1-zJ&WgS=JO(4aaRJ8{(|RrJ
z$;n@n^rdh(RdYZWG2)}XZ*GT_pN7(Am+CA6`yT*QWbD><tW{opb7PyTty7|+dzrF=
zHQY}&`zUcjEAVgoUg1D2($z=E+X{zV5`g1%^M4alOYLFlrQ>!{8I3g@W>%*bu(VzJ
zgJ#vW+YoUA`!DjapPx-KEg=Wv^Fa~MuF0krciO^<U!><whF1}NBEDNK6vA792$LJW
zM<Y<r@W23aBxS#c@Jvp>mVB>x8qf=eoPhT&h!q3xkbjX#b_+>s`>LdLRRKfk_<|s#
z(*e$x1ATG|D8nGP#jar5yWPQ7j;42C8LJwF(nM*-Y-hzX^~*i>!k8Xlv$1v|Y{gfI
zlVBzAjbDf<F=FbeVq;25%C~u)U;azDIsu;JmnG68Vg(~|!fmg{wZhpDK5J?XoDbMu
z-i*JHg)X5Tc!ki*jWr(=Aty7Xl`0;xi`&mtEar2rB6{<ca)Es*WeD}EcO0Z(-A@vA
z?E1qP@9C2>&96I)7XtfHCw0|lE<=$FP4TXkelJ;jru7m@kuMYMnE}ZHoG(svK6tof
zY1afPB5-2|tw+9Z(g;k28+zJz$J@(pi>?zN^A3x7=JSf>MgVnDs7#(cm>F#ww5vOm
z_#Q2tCcdwXdHWgc$Bv{1L_Iknab^k)F|Kb%1)6hFun{Zyv=VSWc*bR|EfJ+GH#pGV
zq8+6f>cwY9)%~1@wSfR6x(!qU1Q^G8Ze~B{9y$@MEIU&?e$bgbUW9tZC~Ln}`i(a4
z1TDj$qV06J^1iAj3ZmngYZM&x^==iAYtGjkC`q#n;Cu}u;HdrHVt8}g)7*DT=;s>6
zIug|kRHLEddd(T|0G&Xn==;+n8LK%vHP=6lV3e;Wew~pxiGSvrbU8{~R)j;%-wFG2
z;zM7@8A+lAMTZ?*dDlcsx9!jI*N|<0sMi{ms|)L$PNKhE^<jDFQGKf7wRTX9nc0-q
zDM(l}w$fagKJPWDC8o}es(e`S8~TbwHdzi-Nv6bA*nWd<FrzV1IX0sko{*+2&-%Cv
z54k)k|J7RnUkLuL68rqdiM|HhkBSAHQ2far3}ZBaz3IEu3!AKKAxXBdDdOf(qje!_
z1C<q`Xw0m`yq^oWaBlyA!9M+*v)N*@r>tsv=wbkjyVWEGX60q?Eib={%SA;xB2PzN
z#1>3&Jpa~h{`es-v}LsZyJIWu;V@gUR^=Y9O^@VBFpLdz#$gT4F5&R#t>T`0=aTl4
zdjc2hE5%0v!`L(PjfQX=?Bx$82%J3x+82CAqZ+FdwR?f%`w<&8?U)ke!*o&1g(&tE
z&Q18c@_`L`SJaYkg36#HnpA_<rUmD>=5MnYW99EpY4q0VJ&IgH+7|*GQ#muD0uP%0
z!g!@G*9Nz9MhU}GWqEo&9gs1)@<D4$rb6n%$?^etNPzXFhtRpQr0Zc7i>ns#+9P&K
z%rA*>n7d;WV=$QK$AzYuzi7>2)JWmS9~W^;!;T&_Fx|CT57~`j?#$8reeOdUv46)s
z1BMeo2^sv=^>Zk`rTbv+A!zUEh-h*M+^Egk<)`#<qxi|;fw7Dqc)Z}`?c6d;t6cet
z*T<p%;vC92Y!}5iI5~BMmg;@4!>DG&dkWbTBCJG}daDI#xB;tWJ1g3@#eWjBZL4d2
z5^(ozT?jRS*^(5G=+JK!zXo1k-$=R&{nal2{6c!RHLj*YT@Who%m@m4Nt;8K%~bFO
z1-87E@NOobu;%_n!IAgmvUW%g%>E|NJ-%l`kUgy~M>S!(VUkLs1EH3)p4Fr@#6uyM
z#v%k|hq{qnMS9^r7Y(e^jWaq>QCA1&M)wfG^314+f{jj{1@-ZehhNQv1#SG$<pSdh
zHy|4mLDcHv){e!0V;$@u%_Dy)Z3cdU{6Hj`%X)^X&c^6sCD5#DI`Qi-N5qRBv+lfo
zPVv$a6`mJ>(xpYNDYDETbM~32Ubz1Ws(dTTviS3U%4E?`l~aEDvtRr=t^#zmn0uiS
z61e{sb4bw{kn~$_cjszMXL#0xW(f51Jk?Fx@19cUU7uga0;X`QLsnYEVAYu64&QVv
z`m6L=slT5F-CWkQ_>sNLuiByEkRS`rw|YFITTn{3V~w-XD_UsysEjc~R}rCgX?a@a
z>$P+Tv5PfFex3gc=W4mJ+n@M6OF<wFfqdu@454CjbYKpFh5Dfzy>Lhyxzv8!^vKn1
zj|LFfTRd7)qpQ<J{kls5cb&*V_9eEx4AN|Xi5`y<8^tIEDX2UNoq)mPGl{&(R8D=x
zmkyPt4uRW^wM_zN#2CKeNg-DR?iW<%YVU3pAz#(sjgyW2taw|>vgPb)(ldxsu=iM7
zh%>z%>Ky0pwLPwK^v_)cFCL7%5#=fBR2vBhev2yJ^3IK_$@}El%MDG~R2r7_V`l_{
zxv?cztH7Xgu~wLnrN8w#RHS5r6z4~peVrO6>mG3a_PifNV{TwQ=%wC)y6*8gm>$DA
zq{9B~MNG4vAEIrf3~9WE<@QeZ&3q13N-fm52{qUFq#Bv5?Qv^?`^XS<B5Kdhg#m%h
zPjh-8TLk`&JlAJgdub9Rzn&~*euap=iqE1f)E3M7m02UeXB!B@r?RC(3-P9pIVRnK
zi!{EU*_nNJdD~wiK9d1Y%Eer@scWM_0a^QAv7@kX$>&qZL$d-%ybjGP19a(fdrQCr
zz~-@qPc7Jvf<%~f(J@e0S=h^r9*BN~n{?vjqVqP?IjGMrX)KC{CNZ$NdB+^}l7mNT
zcH(8V@4+#XjOa`DFa-Oj_AJ^9jOF%sF&ne^mwai>9lpZNSH3RYVcyYyLu{kHuMK0d
zobymIcZIsE=icqe{0mCVV=U(+D4=6}&Mf(GD05)2mF6{t%m+oJ%Zg%F!;)D4*|R0y
zBd}TXeZuOVrUFz~VY$+fPh$BPzZDhBlC`AFIdeGa&2zow8gB5WqNpGM`Hi1XV1&bY
zBg;89!G!p$9(Bbf66q16v9U=YCQ|v28tf#08EyONB<coUb1qHHrw2<9^l6rNH@Jph
zIJDc|FuafGIPDyi?ET0&W^E*0xnIynr?rzo5~*;8agCJ~)khZW4AQMxnz|KV+8%DU
z6#owCTD*S2txT-+ceNo;boxRvxK7y=3sTKGAwg0WYu~n7uzT4bDcG9?m(JbM7YF)O
zS{TaNNSRwnO`J*ytjD+K@1S{}V~dresZt|dBh#{ER6ient`h_kuf<Fue9HjwamQRx
zT2XrurXb8e%m=(xsMa4#KT@)n{k;ppH`giB@TkQ3+6{C2TtUh;J#%6Ui5F(-fqgv?
zyCm>g;$u90%D+Z|_~-rcTo|UJs|v*douY#vb91;BBT9EmNB_XeJa_DpXrykK7({mr
zhX#`fuuT(Ujw(PB#Fw`8Wi*Bvkn(s&C=1Wjo%X}M;pL(4a!1)j`ItqgOMC9}-Y-IU
z<7T3LAM*8DYj}_Hz%{EN3B&zS#Q$;n=7sQkqA_M&X{NBZQXGw6gS!<eIJHWU7fkEA
z$Ng*?mLTU_0Nh4MPNuS%L_=JoZI~GP{aG$$esN|{-yn~5RIPZUu|&RQp4JNFyKrk}
z$74~%%y#%Gc5q!;!!S=>9R(&5doGwdWb?J#EQ8V~#A|NOyn<Q2o6}{)d%=yvoOtgo
zO-qKU70C(XV<8;@>a+CHTaD3pMC)Me!R5338l%eS#>|`+<X@sic|!#HP}v_srLbs6
zO~?AUO+9W!U=v>U%r`?m*Z>g*{;2JAtyhR_#k7wIrpf^bX|Fj-#y@5@pb~`<fkTdm
zqz2K653;IE^s%3YStle%h3ZWbEdYY;$1lL8mCb@0f%R#9zJ4)CLv9DIOBQi+aRy$6
zT^#k68f43U>Bi=02~-;y@cks(Wni)q*8Pqp%vZLzoQ?kX<)Wds>`v_|YuK?d2$+Zk
z-6D^J;5&ku2hW6>m`3a!c6LCpWl#HS<AIo+I$#NWs{ioS_1P;?j&D6GoZh#vA8Q{H
zgJHtC*>=WBsCcA@1qFvK(}oPUrRwK4v@JcTnkL-SMlQ*hln;d6gaAzJ!q6-7E;l>|
z>SN~%jIT7kupVmz5tb2ME6m;;$A{i}`?wv887I^Bk-2C1W^tE&2BA(ugd@FdK_ABf
zwJ|#gr!F#yHL-wePSlWg^M+po+QhF;aqZNPRV3V!@Q%}txW1~`nNQ(D>-f=%^RA`0
zs)RhI0+_Ta&ce8E{3bdvSj`OG%?giVgTGI=@C9)vUw};#b%xu<JeF6_j^o;y?4J{=
zX;#uWyaudRjkyK+jpOFKOYrgH#BKcfKU%OyD3x)9DeV`2^KHsjci2u<U2vLu`%9|e
zR51jI?`~x+STM()j6nGwE8kt98HJ~_twE!F-}LE!_rt99yfJYBUXPt<oH#5T&mo-e
zqnwWJMfTiNOk3vTC<XvaQpLc1h+s3iD>_!)DH2z(`Jt#MzEY(GMi`Ghi%`DB6{;J;
z&>4rEg2;-a(tbctZP}U&LM1Ri1KP8%CJEw;u!P<1jit57G2%`7=;!I403x3)^a<Lj
zhA2n3K=f}4aO3Q=`LW6gte&53qlc9Xva&0EHrK%~s2oujGE1!JxZjEg?HO2{f%!wc
zC@A{VZ%vtk9SYb7WydYJBC)=Evxl*o4XR6#IFP{g3mcb2B3Vr^;qyKx%-q+$fS7ny
z;j79%r*hSVldeaN1M{QMtI==$KGZ*;a}it;mKf9*lVG!6`pxfx*-qMzHKX@E9WNiq
z;xeiYhwtM<V|9ZIUbja_K{_cUsi{xL<l|!DWhUJ9>nVt9xiu=+Jd=R=%Nc0BdN;U|
zxMQ?un6qG@XZy(6awBv5JT+$KzuBDJ8>6DAN#C14gI!>r@jz#jUMjvtVTSY!i{~(4
z4CTi^1M`B(lJQ)S&bHP*O7Uur(>4MHsI;hZpqaGLGG%Q!Hh!-X@mz`c@R356!#&lw
zh^0)?Qgyv+Lt$uHs8P2P=sDiS<V<axYW>E7#36rD0^bGZ6LLshid1K2vR~~dHxE$+
znRa=hd(MPQ_FG4n!b9E`yC|>p)uvHKNV{trHgTpmmoY3KOG~M~$;*%`;r4-^0`FIl
zM4KVQ^}K!+VMwwdqPLBB+;6fMdo$LgMqnPL7bZmrUEEbBZPvJ^?~QR2`WrDaXGLt2
zquB2oA8WdjJCK#zkCK3&!_!0ql9fwPPtyL`-&)m-QXT^rVr|=_X&J|(W0I^cJ?C2|
z4N>3)seubOyB+(na9;Z&#uU!gwBG4F%6#b)A;q0VRzmtM$hrZ^s~@<3V7tvJyDmSY
z3ADcw&NNQqR^%unXhDFbUcR_hqU8MOV|+81naqC@KQts3YOuRh*h287e-4r>o;iyd
zF(JYm(GZw7%0R5*NO?ZDPyNY>FU{2ai}&GqICWtFXMfi-WoDK5=p@Ak-$!X^gkRd0
zj_iQugz94*7I800Xkww)0PJH4m^Y^65+{__<vj(LmDw1{B3H8*rW_SaI#3W-45Dx2
zk)nTEUaUSl`?i?u9@0<S;8~XsxI58;D%lT#W{~7emjmVj+jiR>oC96oii}+kZ>C@q
zlfFj5QSe`aq@0d#Jq=A$(i%5kM-Aq!P7Ujgy{I`#1%Ts$ZFl7Op<_q$ZmCjGMpszf
zHB$n+vJA})C0V!P9ne=^wL@yQIlNly>WSTx><bH$YT$TaJ4N}OOqLS4dANaVR2K#l
zBB^5b95Qtqe4&Aks!A})Fv^l$_fEbJ<@G&(JXEW#2{^t^xu^sA0USpr{(HBYhMsFo
zo4cVN0e<)^JYUafTxFzZ99S#E4nX3^S8m0kgT{nZfmY*iUpWs#6Z+GVTA3{{4>nSw
z<a6f&A%`M}9Q87Gkn#MSA#GVqdHg4hB(M4icMnPD2B|WZRw~$f2^x!1Q*6`?XC_i0
z+z+UV>bA*7b>Qm<q;sLF$!j85@6U#4XB0_4@cj<E^jPv(Yi+TQ%_=HBzP0%e*C8MC
z$&WICV>N5pvnaxsOWNM;vyfSpucR4O0jOn<rCzEF4at3H4sk6L@&j0TviI~?<G7_R
zhlP(?Nb&ISRwU%>H?}{kMhX_xKGPxNBc~xff_q#Xw!&dd6=@E-3BI-B6`^W8SkqJ4
z2Bp!rO>GxX)K`p$>_2hil!hbuS9~dF%QMG_T9%a86hkFUH0%g_<Zw-xeYdkdrUJ63
z!Wx&ted#yAx4ZX5i)d+W<fT;9iWRQz=UK){LcIkwSyANzXjRmo?x^1EFDg^@t@?9}
z?K(Z<rH($*8T4{nZ(C<=aT3S+YyRpX`h()Jp|_&$?gVB!yyM2$5Xj#Z+C0Jk=<pKc
zfNjXOjsE*>L*^O2BSGAU9}_V_<lg0~6S_5Ias_-ND3@h6hl{$MV7yQ8J$_X7<4-Qw
zDA+tS5x}lcm>|7r@E{7G!u;K4Uav<dlSjRTG_Xx+8tnF`=Hj?q&kjS<76oS#TqsAR
zs?ohia47th${y#c!(_KH4d$UX9^U2Ji+msJ!$jlX1l;90I=+icb#vy0aQxoB7m{0m
zNijUzJY;IyXUuh6UA-S=>r-0RlC@#9W@aogs32@ZzOVIRL;No(7^X|YY01zTs7-Rr
zjsCVVVzkOi#ol3eS|1re_=Z@{L*Sn30M4x#w;FSLlE=^;Yc~N5dqtjk=!?H~X*q6*
zFjhrTe_{sBl}GfgigjjQiDMx=fuC-`a*^QgnMZn<SH?B!fwb4oF5=IF95DNz*55&e
zMer*!cj4t*SF@rh@Nhp_$K-evF<0=s2J08yNh7#a3u<}lnut~KM;8B5AP9A!wasrP
znNdjAh8?%fUXBRtml3}y2$0c177(wtsjR3EgMRy|m-eQ;3=x<Wg96L0W)5wH)jLf+
ztLgYqBT(>G1AP5pl|_4@5K1{<J}o>!=kb`eK=#w4Hta?DS)|kX<92?`Og0Sf^5>^2
zUkXQ1vV+%Cx-%KEbhbPSlmEr5ShkW|gi>s1jm4O!2S=3k)xU}W=2y^WnI91Sj6<h$
z3$+844YU<+SWZ)q&wg~76Q*lKfyF}?4rip=e(w*rh>aw;x^el0<Oi$~8+%4=V4e92
zC82oFQ^G4%4%=^*<~)Tz=mYZ$<*2A1<~KTBU3grdkz59Sk#ljl(bja29yU@%vN6z6
z?g!y(3kM%1Th2CHALp7w?%ckyQo<GWWg^p?(nN|opjhZcD4WH(Dj}3ir`TM0R_$Fz
zV%;7zt+99hd2*TwnwL=(faDP3QxpML+x6_L8~7-FW1@RXY;b?VTRMC#VGHaN3l|77
ziKKi6z0MoGc{-^?!@`4F<|sv<0bO@G-8qsItjY44oDimLi5+LDeE8n1-@Zcvb1+g~
z{R;l+N_91`mBXeU`ut+mIjlu#q!6O|aD1rA5Fg(C<$637_pX^kdR-cc+s!G@&Nb44
zsYl&5119T|``qJoEDKkyuSMRt#!Lbg{5kRn?nF+A#a4-^RdV?Wq5oqV$qayLP(G9#
zbNKVkS3@aveF<^*2UT5FvSD`0v*V|zv3%OTYv3%hwqU}PuUPFd^4hgf>pXOaS+@h;
z=G`RU$YR{=inXRJBTSmXr=)4n-fl<J<@1!6Qm1bG->LKsVV$#8zBC@!&+v+5Sj>tw
z<Wlz{;PnXvC`{5dtw+sAsrUQx4o*hYA$Io@8Z+`O@ud8M$Y_a)U)W4psmK~QrXr*1
zDRLNA_U4YfM7aM-fc7V2P#dEfT+Z^5^hp;x50~`56vq%_ymu61IgUvE8P)427w12c
zrZ<~Y@k!5W5U>&h@t}i6-vD~3GNHS?aU7A=Xi(*R1K8LhH;pc6I;0!3SNs4Iqb~@C
zyQW=$$|0A?hfLDtAHJj_LezN)Tt3`WgM+Sp-gu{j-f_@ovq_yNjMPP4gk&`KpXPTn
zVDlykVx-{P!3YmqD%;-2te&b9*vaASUS<3~Tmt2s1y|`?o3aM{ylCGr^)OLr1C@v=
zuVi$iug14=eHNsYUq)PCyf=#=@c*?)Ghh<u;39c)wvfp%PWJ21r#?anrQ<_KB|0TS
z)SLwl<4ajbbkAFjjFn*x3jGMkqZ8`LzBL;Y7VBfM<-3_LIRz8u%G1jhdcSt(%Hhl<
zU4!PcpzjQ4LRv%@BrPjjYq|aWgXjS)x`fw>t-w+$U80?T4c_BG|G<pI&5z$Uj{yPj
zd9P=Qd2(pu^D0=_+<8p{HFPX9Wm7fF)?;p_q<+sh3CR76t5cv8p#vSDHV$x!8k#Of
z!F$Ok+$0wB?1mRk^Ps81&;r9qQS~loyj@aPD1%idZFV<0rK$+uTLBr=e@~Y}c)JA0
zSMG(qz2c|Jrqa$1lYJ41)n0cKZH0f=G$~${%hcz)a_aI%#iIDiF8cgz>`_m|PBs8O
zo_~z@b$($69VT82mIg0@?KbaaP2*(PkL$KUQEE_Pxtpdp9@ZPuJ>sL$=9gHhL~<U}
zUtqA=5Q)`g*PF(loBgDBk|fjjLJG`#&tAKvi`QAvHDOT2M^Hc;J=hnvQO#51Lwlr#
zSq!u0z+|NFDmKGur+A$9(UZaJdL@v!&g%5`*fl0r)T#%-04Oa@g~|Tn?HSo)21vvN
zw^sVD^Fke*jMc{6Y{+^yfbr<8_k0s1&QXVlYfIr0Z3Mg*DkJtPA9MH?mNb)jLYW<T
zhy(X0gZMb28NOnNhe^?q0L8f5_5y>BZmi?M*spt_R;!Ef?Uy1x`efb9j|#9Im)uH!
z85#WCr~GWxdIhb==qtX-L$sq&ONw8X_Z>y0{=fGcdQZ%~2Bl%b4*Xn6sDB<6I#UiT
zON6;nHy7!eU1G>ou~)7Iq~`cAYNPLkzuSpjRUn+Iq(hTguOTOw&_0Dk?ST!I1(M~K
z(-w<GPfsC#_E9&^N*jo$vM9=k(djYdcz{sX9h7`7p(HZ8%?or@I;n?Rm6WOpsNxZw
zu|r5c9rKkS1x0U8OJ<wEeHi0LuYp!!=tqSxz_!Y!g<TPE({g3>guyEl`5`Qfz>xjQ
z?wu-_Co;~;O64dU4<!fJa0+!O|Kp0TE7ww>$=ky#1w4QBIbblFZ}}aCV`anG1y<?O
z?<{g>AojGSLVBd1Dwr?eN*GF@S(r#n$2K|^Ym{l66u@m!oAYMD7H7G$g>D)J1<B#H
z^>@Ah`v!A^s!Ya{96*IY?QPo${Kf*kyd{YW_<rVuT#ehqNo0`aaHnx$LFNR!f3$AN
zCn2xvMVcb>cLsBz6wwA-(YyRdZS~=1L}f4K?FaDu2lEUj_nO8LcQ2tfj!`vgj)>!l
zv|6*iZ!wy}1D%=5mMoYkwf^<E!xA(y41$Ut`pFa@o?4_x7cFXni5s0$a0xe3;=lm=
zNp!x34rw&{j+ex8T*o4$iH_d=+sZqV;fNQZiZ;K%oK)oE?FLyo66wZfoZzT9(oAlC
z9uw~U(P7^mJioh8*qGyDFeY<s*;OCy#Ufg7iyG^_{IVDC`{uh`Y|wL7^~BT3x5)%$
z-LeL;Bbvj9aE#63X)gP{874_GGN{WJ!zx?*^LV;HLURit&Asx)qIq)PD-XI<{{qhs
zAiRwyc(^(rUKczMhv}fc0e|mN`n@E|3%4k0$mIF>el>QQM-o|JVEtF~{EfQzFy-M`
zX$%G5K0?dy9WYzOo#r*4)XFxhza-;fxr}fLk{k60tYlZ9uEKqV^O3SsEc({C){76Y
zsX_vEqx<8H57y0^f*(Dlp`-K~Lya+E=N@T?KLe?LJp=?Ssm+gu3B_ry0N3Lr27drX
zZtTzw-oPEYL6%$37dn7#Opvj%-aow0Rk_)}?2gsqAmp~`8MhdC<w!smEx`e6y57F9
z*riLh;#9|MnK|2yZb&^KcjK<V6Nc{e8C_=Q{=I~%yRQv@b2(pYMgA#us-1pC46H`!
zCy!$g5)5Y|Pa7oz{#@PUB>U_Z!;f(qzVG!j4x7e07q=xM6#>)I?&GIT@jO=5c*<@J
z;k|1p`WSjVI11N9>y_=q4RF|;%!S}C;Z-lzjCm)(=P$ewanNUL^CNnmt-vDJJ8ZPr
zrGq)hC&EiNB;=T+(4Q1YdD=a{%p9=eMkV%@24v0<y)8&6t8xIr-s}Qg@2=BO9{N^D
z&cT_mFfGQlxLoD4%zhXHtbn!XAF+zk(h%D&A|LzehsM$UD#W#CG777+(V2`(aZ@;g
zIJxM1z2O2r_JJi<YiznSCq>NcMwzaY)`(DS$&$fJncNM{LHP(*Fe>q3zurVQ2(984
zgGxz!Dhh##y;gPa@=Hbtw<~IJXX8*^tDvCsJ=5eLicUUhhm7WI1RLqXIa*Lzq;bIg
z#<kJ*9f(Q^%fqp^SO50mFUvcP&&ml$>PU<iiGOIN=yM!Tq*A<Uyd+ec&w+NXbq63S
zpXsQOsk`>0#8q_T0WzP3yud%U@b}mRN|BkvuSsLHp}8%BWIDW0J<uhTG-wvY36DYE
z-2v+Lv&pUu3QmS#QdLq}Bf=i@skSlF!1Ft^TASebfvVoYZ+gXFp3EL7soT|-vhjg7
z8jXu-pt@#YdqvV2XY0LGGybR-MbI&_gB0r6#uC{0lhhI{S{-u$^ryoYE?D(>{8q}a
zu9AP&*GOP_N0})^w;!a=I!}(!8q435?H=b)lWotT=POW|PRRW1n5MAtA(={SuF@?v
z$~V>ZftJnsU2v4l&5S?mLH>6|B}_%x1GRfYZo`(s^AR$)A2sK$W&!1=IrhvXdWphG
z&Y2%~Es_e^<q*fR0mt_C$&(&NmY}3pH5I($Kkw&4apX0V{Kn88!b=d^)ceZZ97B-t
zjOL())14ps_m+1hr?)#b>Hel(1TAXtPMxR{wbMi`EMie8F0*Hx`9<8NlVB`0%tp~>
zcL}}&Jt-r=>#I;^wi5bc5VQ44nPaDu2oV`c_aackMzWl0z#qp`iJ!}|y4cH$jk@v`
z6~o7ZYs82@B>e7{dPr~Uvzc{Q2)Y>)s}uIvN+5cKD<~dD*Nq2OUsHigPR<@$d*_F?
zX{Afyi=T*5wmy>>oE-VrX?%SgXTQ`p1Q1#P&j_fmn+L8ZHJNEyL_^C!VSyysSVJka
ziC3?ok4s@}bGi=$N(_=sj=hf?Fj!_59<?jdTSPPk-dh*~Q&?BPUb|jW`zmer=Jc{?
zkPMPWC-_4#wu~qonD5yMmV6WV)*LGSz^$t@*b2NK_4W-88a&T>p7<$IsR(E0^#E?I
zt$*v*fIC3kJq7`gt<y-$L>)pQb+-N!tfQOhF|hsD%)0m!gfUh4G+rmSH1PRc;S-=J
zQrFRgbSQg(CzB)sT+i_DF?MF9>VR$#(~3TK<Uok+sB0MMnvpIDEG{j{);Z;n1Rm)5
z8s!ht@HnbD+7f#7)GYs^3o~_0E5-`P=Y!wqir6s9c~*>aM={qdULCe_LluAYntv$i
z<Ez$_SB<Y-j9b+q)<rc!>GFhRhEv(FwffMKjRXtj7*J`sTzX#EKzC18eu_5SIXJJ=
zMD9qsW%?ImBN3*8yK|QNm8tCr(}>gFvWCOS<0sbbydnYDFFFLX(KBsR--3snwpAC2
z2&@OG8%e*_ghB~1DhyYmL7vq7`3(p#&F13XEOA~KU1SGbHUrdozbz_B(PLhJg@XXc
zr*)XpQ`=WMvzmM!>erh$0P7dy-(J?66f)kaMJ=vLnjAOKgZAm)D%?zdY|ccrmRqn~
zzL}5<Vp<FgD)QUIvNHqEvyLSvRU#I-cOplRvfEAE5`g*CII9C2II}v7EV6sStQ#H|
z<S`K*3@?pD8&5KIA<2A_v_ZehJ-KyUMTd?tW&LjaAEmQnl`Y7ayfVkqB*mS`4^`)5
zjb5FveTXLdvW8lL`7<)wfe&B4^1WPI+1<aAMxCW`HOr$c9OjmEi$3W@Z!PX+?(p!X
zbhL^TP7Pb8zikN=XZ|(9v5x3dmWs5Q2+n-QXK89pPza9t-6w9*GMF1y999*IZ6a1M
zm>qMbjt9Kn+Xf#m3z$|LAm`MouT7E~!1Z8X8fKfO&{<E#K3Z2ARH?nizV8MvLPQ=6
z0Y-|ZuPqn4WIg^wZ|k$|VtA@S_@V&E2fiapiio-mN!!6_>(L}G0P8$5gn(Z`c10Z5
z9&F}x0olNM<o5(URT3!Fu7|lZ`L4i!_ERtQdbAhU^OGVmIqSmTkHYFmU+lL`)*Hna
z;Jz~uYT!MlMDZaEw5$y@zqn7l`NrKG+7I{5Z|j1LvO^((<ZMQsTERr>a*8fA4ZT=&
zVf-_{C29AeI(@P#8>I}4e(e*9-~R*FqO2UsX*osj6eC8^>UiK;;6;#^G}bT@G0jFz
z_c-;K;8NqZ%fi@NI1+*{w6ctHFmyxLO0LdhGdsS-o3hwTz%{#qd6mwm18+1C1>7%q
z_<Bc5mK_nR$f;=JxY8i7?oE`&<3Q1hC!R8DIpP^`2ds;_w|-aM_ntM@bjY$^^!%qE
ztb|-4XA(X$*6A3>v2xyjzGv0_-ZF0n+f|Q%#3Na=K@SPMpNprGm>5hF^H{6x|J>qO
zB@LWE*<p}AeQWrr6bcOYAjhTvSQjK=9o#DR_X{>`y3VMIQx6$bfbv1APD@ozzp~Ok
zpf5vf-@g=^a@DF4W<3@RnE{Rmwtu2K_4x;)75klt+-Lb=jsoke?Hs<alhKJYb)>b#
zyVHysN_U=WC{!Gh3=o={n^{=|ti2-Ad))!uWDpL2xs`i?<NwPhwH)~4ix^VMf?;Hm
ziS$4Qe7p)8?DIP}V2#YA0DlL@@c)_r{1+ei_j7n2c8d|=@zYuClFZqk<8-i#{N=ol
z>c+ZVmZP7VI86V=`!9W?+w-Z3>Dv(K>H*ZnfweKPzZ&jI!oK2l9D8V{A`blWKhINO
ze_-37tr|WHli4m=^Y$zqr1=&$-HU+yj}!e?co%Ny?@%9F=(1h&_tIMr*asshVlF=5
z`2W@B$XwhqD6-K&Z{jW4f;k5Vj^_g0nbk;5|DkUs85W8&)&ti6%O<?s(W!Wym&B4)
zmQAMKQm#haWvt{e%#sjw!Z{@mF<De(36i*eo^91u#;agj_-DTQUwmMlq3<mZ?H0lC
zs&sbPE|UTiSQpqg-d6wnMPs^u0h1&8VfBv=^nVIq{SDeR_^moQP-Ur`-Q0H{5|Hf9
ztP7TWTS_SvMF1lT#gp?uer;iDTxMXn`2T;o|DVq{jPC$aY2`%v3Wi6%!pr)9=8eGh
z{jWCg_t3;6M3K}{+HDD4!asu3gr5Fa-hcTMdR$=An*JePdCQ>LH;yOyFaQ6m?$r*l
nGgy4O3@WJ9;bPYIKVOgkl|MCC1_@3c4W4>JV-;y#7xDiA3vybK

literal 0
HcmV?d00001

diff --git a/aes.cpp b/aes.cpp
new file mode 100644
index 0000000..1aa92c4
--- /dev/null
+++ b/aes.cpp
@@ -0,0 +1,523 @@
+#include "aes.hpp"
+
+/* static function prototypes */
+static byteArray word2bytes (word input);
+static word bytes2word (byte b0, byte b1, byte b2, byte b3);
+static void circular_left_shift (byteArray &bytes, int shift_amt);
+static void circular_right_shift (byteArray &bytes, int shift_amt);
+static byte mult (const byte ax, const byte bx);
+static byte xtimes (const byte bx);
+static void printState (byteArray &bytes, std::string name);
+
+AES::AES (const byteArray& key)
+	: Nb(4)					// This is constant in AES
+	, Nk(key.size() / 4)	// This can be either 4, 6, or 8 (128, 192, or 256 bit)
+	, Nr(Nk + Nb + 2)
+	, keySchedule(Nb * (Nr+1), 0x00000000)
+{
+	// Check the arguments
+	if (!(Nk == 4 || Nk == 6 || Nk == 8))
+		throw incorrectKeySizeException();
+
+	// Generate the Key Schedule
+	KeyExpansion (key, keySchedule);
+
+#if 0
+	std::printf ("Key Schedule\n");
+	for (int i=0; i<keySchedule.size()/4; ++i)
+	{
+		for (int j=0; j<4; ++j)
+		{
+			byteArray temp = word2bytes (keySchedule.at(i*4+j));
+			std::printf ("%.2x %.2x %.2x %.2x ", temp[0], temp[1], temp[2], temp[3]);
+		}
+		std::printf ("\n");
+	}
+#endif
+}
+
+byteArray AES::encrypt (const byteArray& plaintext) const
+{
+	// Make sure that the plaintext size is a multiple of 16
+	if (plaintext.size() != 16)
+		throw incorrectTextSizeException ();
+
+	int round;
+	byteArray state (plaintext);
+
+	/* Round 0 */
+	//std::printf ("Round 0\n");
+	//printState (state, "input");
+	AddRoundKey (state, GetRoundKey (0));
+
+	/* Round 1 to Nr-1 */
+	for (round=1; round<Nr; ++round)
+	{
+		//std::printf ("Round %d\n", round);
+		//printState (state, "start");
+		SubBytes (state);
+		//printState (state, "sbyte");
+		ShiftRows (state);
+		//printState (state, "srows");
+		MixColumns (state);
+		//printState (state, "mcols");
+		AddRoundKey (state, GetRoundKey (round));
+	}
+
+	/* Round Nr */
+	//std::printf ("Round %d\n", round);
+	//printState (state, "start");
+	SubBytes (state);
+	//printState (state, "sbyte");
+	ShiftRows (state);
+	//printState (state, "srows");
+	AddRoundKey (state, GetRoundKey (round));
+
+	return state;
+}
+
+byteArray AES::decrypt (const byteArray& ciphertext) const
+{
+	// Make sure that the plaintext size is a multiple of 16
+	if (ciphertext.size() != 16)
+		throw incorrectTextSizeException ();
+
+	int round = Nr;
+	byteArray state (ciphertext);
+
+	/* Round Nr-1 */
+	AddRoundKey (state, GetRoundKey (round));
+
+	/* Round Nr-2 to 1 */
+	for (round=Nr-1; round>0; --round)
+	{
+		InvShiftRows (state);
+		InvSubBytes (state);
+		AddRoundKey (state, GetRoundKey (round));
+		InvMixColumns (state);
+	}
+
+	/* Round 0 */
+	InvShiftRows (state);
+	InvSubBytes (state);
+	AddRoundKey (state, GetRoundKey (round));
+
+	return state;
+}
+
+void AES::KeyExpansion (const byteArray& key, wordArray& w) const
+{
+	const static word Rcon[] = {
+			0x00000000,
+			0x01000000,
+			0x02000000,
+			0x04000000,
+			0x08000000,
+			0x10000000,
+			0x20000000,
+			0x40000000,
+			0x80000000,
+			0x1b000000,
+			0x36000000,
+	};
+
+
+	int i;
+	word temp;
+
+	/* Copy the key bits into the beginning of the word array */
+	for (i=0; i<Nk; ++i)
+		w[i] = bytes2word (key[i*4+0], key[i*4+1], key[i*4+2], key[i*4+3]);
+
+	for (i=Nk; i < (Nb * (Nr+1)); ++i)
+	{
+		temp = w[i-1]; // copy the previous word into temp
+
+		if (i % Nk == 0)
+			temp = SubWord (RotWord (temp)) ^ Rcon[i/Nk];
+		else if (Nk > 6 && i % Nk == 4)
+			temp = SubWord (temp);
+
+		w[i] = w[i-Nk] ^ temp;
+	}
+}
+
+void AES::SubBytes (byteArray &state) const
+{
+	int i;
+	static const byte sbox[] = {
+			/* 0     1     2     3     4     5     6     7     8     9     A     B     C     D     E     F */
+			0x63, 0x7c, 0x77, 0x7b, 0xf2, 0x6b, 0x6f, 0xc5, 0x30, 0x01, 0x67, 0x2b, 0xfe, 0xd7, 0xab, 0x76,
+			0xca, 0x82, 0xc9, 0x7d, 0xfa, 0x59, 0x47, 0xf0, 0xad, 0xd4, 0xa2, 0xaf, 0x9c, 0xa4, 0x72, 0xc0,
+			0xb7, 0xfd, 0x93, 0x26, 0x36, 0x3f, 0xf7, 0xcc, 0x34, 0xa5, 0xe5, 0xf1, 0x71, 0xd8, 0x31, 0x15,
+			0x04, 0xc7, 0x23, 0xc3, 0x18, 0x96, 0x05, 0x9a, 0x07, 0x12, 0x80, 0xe2, 0xeb, 0x27, 0xb2, 0x75,
+			0x09, 0x83, 0x2c, 0x1a, 0x1b, 0x6e, 0x5a, 0xa0, 0x52, 0x3b, 0xd6, 0xb3, 0x29, 0xe3, 0x2f, 0x84,
+			0x53, 0xd1, 0x00, 0xed, 0x20, 0xfc, 0xb1, 0x5b, 0x6a, 0xcb, 0xbe, 0x39, 0x4a, 0x4c, 0x58, 0xcf,
+			0xd0, 0xef, 0xaa, 0xfb, 0x43, 0x4d, 0x33, 0x85, 0x45, 0xf9, 0x02, 0x7f, 0x50, 0x3c, 0x9f, 0xa8,
+			0x51, 0xa3, 0x40, 0x8f, 0x92, 0x9d, 0x38, 0xf5, 0xbc, 0xb6, 0xda, 0x21, 0x10, 0xff, 0xf3, 0xd2,
+			0xcd, 0x0c, 0x13, 0xec, 0x5f, 0x97, 0x44, 0x17, 0xc4, 0xa7, 0x7e, 0x3d, 0x64, 0x5d, 0x19, 0x73,
+			0x60, 0x81, 0x4f, 0xdc, 0x22, 0x2a, 0x90, 0x88, 0x46, 0xee, 0xb8, 0x14, 0xde, 0x5e, 0x0b, 0xdb,
+			0xe0, 0x32, 0x3a, 0x0a, 0x49, 0x06, 0x24, 0x5c, 0xc2, 0xd3, 0xac, 0x62, 0x91, 0x95, 0xe4, 0x79,
+			0xe7, 0xc8, 0x37, 0x6d, 0x8d, 0xd5, 0x4e, 0xa9, 0x6c, 0x56, 0xf4, 0xea, 0x65, 0x7a, 0xae, 0x08,
+			0xba, 0x78, 0x25, 0x2e, 0x1c, 0xa6, 0xb4, 0xc6, 0xe8, 0xdd, 0x74, 0x1f, 0x4b, 0xbd, 0x8b, 0x8a,
+			0x70, 0x3e, 0xb5, 0x66, 0x48, 0x03, 0xf6, 0x0e, 0x61, 0x35, 0x57, 0xb9, 0x86, 0xc1, 0x1d, 0x9e,
+			0xe1, 0xf8, 0x98, 0x11, 0x69, 0xd9, 0x8e, 0x94, 0x9b, 0x1e, 0x87, 0xe9, 0xce, 0x55, 0x28, 0xdf,
+			0x8c, 0xa1, 0x89, 0x0d, 0xbf, 0xe6, 0x42, 0x68, 0x41, 0x99, 0x2d, 0x0f, 0xb0, 0x54, 0xbb, 0x16,
+	};
+
+	for (i=0; i<state.size(); ++i)
+		state[i] = sbox[state[i]];
+}
+
+void AES::InvSubBytes (byteArray& state) const
+{
+	if (state.size() != Nb * 4)
+		throw badStateArrayException ();
+
+	int i;
+	static const byte inv_sbox[] = {
+			/* 0     1     2     3     4     5     6     7     8     9     A     B     C     D     E     F */
+			0x52, 0x09, 0x6a, 0xd5, 0x30, 0x36, 0xa5, 0x38, 0xbf, 0x40, 0xa3, 0x9e, 0x81, 0xf3, 0xd7, 0xfb,
+			0x7c, 0xe3, 0x39, 0x82, 0x9b, 0x2f, 0xff, 0x87, 0x34, 0x8e, 0x43, 0x44, 0xc4, 0xde, 0xe9, 0xcb,
+			0x54, 0x7b, 0x94, 0x32, 0xa6, 0xc2, 0x23, 0x3d, 0xee, 0x4c, 0x95, 0x0b, 0x42, 0xfa, 0xc3, 0x4e,
+			0x08, 0x2e, 0xa1, 0x66, 0x28, 0xd9, 0x24, 0xb2, 0x76, 0x5b, 0xa2, 0x49, 0x6d, 0x8b, 0xd1, 0x25,
+			0x72, 0xf8, 0xf6, 0x64, 0x86, 0x68, 0x98, 0x16, 0xd4, 0xa4, 0x5c, 0xcc, 0x5d, 0x65, 0xb6, 0x92,
+			0x6c, 0x70, 0x48, 0x50, 0xfd, 0xed, 0xb9, 0xda, 0x5e, 0x15, 0x46, 0x57, 0xa7, 0x8d, 0x9d, 0x84,
+			0x90, 0xd8, 0xab, 0x00, 0x8c, 0xbc, 0xd3, 0x0a, 0xf7, 0xe4, 0x58, 0x05, 0xb8, 0xb3, 0x45, 0x06,
+			0xd0, 0x2c, 0x1e, 0x8f, 0xca, 0x3f, 0x0f, 0x02, 0xc1, 0xaf, 0xbd, 0x03, 0x01, 0x13, 0x8a, 0x6b,
+			0x3a, 0x91, 0x11, 0x41, 0x4f, 0x67, 0xdc, 0xea, 0x97, 0xf2, 0xcf, 0xce, 0xf0, 0xb4, 0xe6, 0x73,
+			0x96, 0xac, 0x74, 0x22, 0xe7, 0xad, 0x35, 0x85, 0xe2, 0xf9, 0x37, 0xe8, 0x1c, 0x75, 0xdf, 0x6e,
+			0x47, 0xf1, 0x1a, 0x71, 0x1d, 0x29, 0xc5, 0x89, 0x6f, 0xb7, 0x62, 0x0e, 0xaa, 0x18, 0xbe, 0x1b,
+			0xfc, 0x56, 0x3e, 0x4b, 0xc6, 0xd2, 0x79, 0x20, 0x9a, 0xdb, 0xc0, 0xfe, 0x78, 0xcd, 0x5a, 0xf4,
+			0x1f, 0xdd, 0xa8, 0x33, 0x88, 0x07, 0xc7, 0x31, 0xb1, 0x12, 0x10, 0x59, 0x27, 0x80, 0xec, 0x5f,
+			0x60, 0x51, 0x7f, 0xa9, 0x19, 0xb5, 0x4a, 0x0d, 0x2d, 0xe5, 0x7a, 0x9f, 0x93, 0xc9, 0x9c, 0xef,
+			0xa0, 0xe0, 0x3b, 0x4d, 0xae, 0x2a, 0xf5, 0xb0, 0xc8, 0xeb, 0xbb, 0x3c, 0x83, 0x53, 0x99, 0x61,
+			0x17, 0x2b, 0x04, 0x7e, 0xba, 0x77, 0xd6, 0x26, 0xe1, 0x69, 0x14, 0x63, 0x55, 0x21, 0x0c, 0x7d,
+	};
+
+	for (i=0; i<state.size(); ++i)
+		state[i] = inv_sbox[state[i]];
+}
+
+
+void AES::ShiftRows (byteArray& state) const
+{
+	if (state.size() != Nb * 4)
+		throw badStateArrayException ();
+
+	int r, c;
+	byteArray temp (Nb, 0);
+
+	for (r=0; r<state.size()/4; ++r)
+	{
+		// Copy into temp
+		for (c=0; c<Nb; ++c)
+			temp.at(c) = state.at ((c*state.size()/4)+r);
+
+#if 0
+		std::printf ("before cls(%d)=", r);
+		for (c=0; c<Nb; ++c)
+			std::printf ("%.2x", temp.at(c));
+		std::printf (" -- ");
+#endif
+
+		// CLS 0, 1, 2, 3
+		circular_left_shift (temp, r);
+
+#if 0
+		std::printf ("after cls(%d)=", r);
+		for (c=0; c<Nb; ++c)
+			std::printf ("%.2x", temp.at(c));
+		std::printf ("\n");
+#endif
+
+		// Copy back to state matrix
+		for (c=0; c<Nb; ++c)
+			state.at((c*state.size()/4)+r) = temp.at(c);
+	}
+}
+
+void AES::InvShiftRows (byteArray& state) const
+{
+	if (state.size() != Nb * 4)
+		throw badStateArrayException ();
+
+	int r, c;
+	byteArray temp (Nb, 0);
+
+	for (r=0; r<4; ++r)
+	{
+		// Copy into temp
+		for (c=0; c<temp.size(); ++c)
+			temp[c] = state[(c*4)+r];
+
+		// CRS 0, 1, 2, 3
+		circular_right_shift (temp, r);
+
+		// Copy back to state matrix
+		for (c=0; c<temp.size(); ++c)
+			state.at((c*4)+r) = temp.at(c);
+	}
+}
+
+void AES::MixColumns (byteArray& state) const
+{
+	if (state.size() != Nb * 4)
+		throw badStateArrayException ();
+
+	const static byte transform[] = {
+			0x02, 0x03, 0x01, 0x01,
+			0x01, 0x02, 0x03, 0x01,
+			0x01, 0x01, 0x02, 0x03,
+			0x03, 0x01, 0x01, 0x02,
+	};
+
+	int r, c, i, j;
+	byteArray temp (Nb, 0);
+	byteArray result (Nb, 0);
+	byte total;
+
+	for (r=0; r<4; ++r)
+	{
+		/* Get this column */
+		for (c=0; c<Nb; ++c)
+			temp[c] = state[(r*4)+c];
+
+		/* Do the Multiply */
+		for (i=0; i<4; ++i)
+		{
+			result[i] = 0x00;
+
+			for (j=0; j<4; ++j)
+				result[i] = result[i] ^ mult (transform[i*4+j], temp[j]);
+		}
+
+		/* Copy back into state matrix */
+		for (c=0; c<Nb; ++c)
+			state[(r*4)+c] = result[c];
+	}
+}
+
+void AES::InvMixColumns (byteArray& state) const
+{
+	if (state.size() != Nb * 4)
+		throw badStateArrayException ();
+
+	const static byte transform_inv[] = {
+			0x0e, 0x0b, 0x0d, 0x09,
+			0x09, 0x0e, 0x0b, 0x0d,
+			0x0d, 0x09, 0x0e, 0x0b,
+			0x0b, 0x0d, 0x09, 0x0e,
+	};
+
+	int r, c, i, j;
+	byteArray temp (Nb, 0);
+	byteArray result (Nb, 0);
+	byte total;
+
+	for (r=0; r<4; ++r)
+	{
+		/* Get this column */
+		for (c=0; c<Nb; ++c)
+			temp[c] = state[(r*4)+c];
+
+		/* Do the Multiply */
+		for (i=0; i<4; ++i)
+		{
+			result[i] = 0x00;
+
+			for (j=0; j<4; ++j)
+				result[i] ^= mult (transform_inv[(i*4)+j], temp[j]);
+		}
+
+		/* Copy back into state matrix */
+		for (c=0; c<Nb; ++c)
+			state[(r*4)+c] = result[c];
+	}
+}
+
+word AES::SubWord (const word& input) const
+{
+	byteArray bInput = word2bytes (input);
+
+	SubBytes (bInput);
+
+	return bytes2word (bInput[0], bInput[1], bInput[2], bInput[3]);
+}
+
+word AES::RotWord (const word& input) const
+{
+	byteArray bInput = word2bytes (input);
+
+	circular_left_shift (bInput, 1);
+
+	return bytes2word (bInput[0], bInput[1], bInput[2], bInput[3]);
+}
+
+wordArray AES::GetRoundKey (const int round) const
+{
+	wordArray temp (4, 0);
+
+	temp[0] = keySchedule.at (round*Nb + 0);
+	temp[1] = keySchedule.at (round*Nb + 1);
+	temp[2] = keySchedule.at (round*Nb + 2);
+	temp[3] = keySchedule.at (round*Nb + 3);
+
+#if 0
+	std::printf ("ksch%d   ", round);
+	for (int i=0; i<4; ++i)
+	{
+		byteArray btemp = word2bytes (temp[i]);
+		std::printf ("%.2x%.2x%.2x%.2x", btemp[0], btemp[1], btemp[2], btemp[3]);
+	}
+	std::printf ("\n");
+#endif
+
+	return temp;
+}
+
+void AES::AddRoundKey (byteArray& state, const wordArray& w) const
+{
+	int i, j;
+
+	for (i=0; i<w.size(); ++i)
+	{
+		byteArray wBytes = word2bytes (w[i]);
+
+		for (j=0; j<Nb; ++j)
+		{
+			//std::printf ("state.at(%d) ^= wBytes.at(%d) -- %.2x ^ %.2x = %.2x\n", i*Nb+j, j, state.at (i*Nb+j), wBytes.at(j), state.at(i*Nb+j) ^ wBytes.at(j));
+			state.at(i*Nb+j) ^= wBytes.at(j);
+		}
+	}
+}
+
+/******************************************************************************
+ *                              STATIC FUNCTIONS                              *
+ ******************************************************************************/
+
+static word bytes2word (const byte b0, const byte b1, const byte b2, const byte b3)
+{
+	word output;
+	output = (0x00000000)  | b0;
+	output = (output << 8) | b1;
+	output = (output << 8) | b2;
+	output = (output << 8) | b3;
+
+	return output;
+}
+
+static byteArray word2bytes (const word input)
+{
+	byteArray output (4, 0x00);
+	output[0] = (input & 0xff000000) >> 24;
+	output[1] = (input & 0x00ff0000) >> 16;
+	output[2] = (input & 0x0000ff00) >> 8;
+	output[3] = (input & 0x000000ff) >> 0;
+
+	return output;
+}
+
+static int ring_mod (const int number, const int mod_amt)
+{
+	int temp = number;
+
+	while (temp < 0)
+		temp += mod_amt;
+
+	return temp;
+}
+
+/* ROL all of the bytes in @bytes by @shift_amt */
+static void circular_left_shift (byteArray &bytes, int shift_amt)
+{
+	int i;
+	byteArray temp (bytes.size(), 0);
+
+#if 0
+	std::printf ("BEFORE CLS(%d): ", shift_amt);
+	for (i=0; i<bytes.size(); ++i)
+		std::printf ("%.2x", bytes[i]);
+	std::printf ("\n");
+#endif
+
+	for (i=0; i<temp.size(); ++i)
+	{
+		int tindex = i;
+		int bindex = (i+shift_amt) % bytes.size();
+		//std::printf ("temp[%d] = bytes[%d] = %.2x\n", tindex, bindex, bytes[bindex]);
+		temp[i] = bytes[(i+shift_amt) % bytes.size()];
+	}
+
+#if 0
+	std::printf ("AFTER CLS(%d): ", shift_amt);
+	for (i=0; i<temp.size(); ++i)
+		std::printf ("%.2x", temp[i]);
+	std::printf ("\n");
+#endif
+
+	for (i=0; i<bytes.size(); ++i)
+		bytes.at(i) = temp.at(i);
+}
+
+/* ROR all of the bytes in @bytes by @shift_amt */
+static void circular_right_shift (byteArray &bytes, int shift_amt)
+{
+	int i;
+	byteArray temp (bytes.size(), 0);
+
+	for (i=0; i<temp.size(); ++i)
+		temp[i] = bytes[ring_mod (i-shift_amt, bytes.size())];
+
+	for (i=0; i<bytes.size(); ++i)
+		bytes[i] = temp[i];
+}
+
+static byte xtimes (const byte bx)
+{
+	const byte mx = 0x1b; /* x^8 + x^4 + x^3 + x + 1 */
+
+	/* See Notes Pg 36. This is if b7 == 1 */
+	if (bx & 0x80)
+		return (bx << 1) ^ mx;
+
+	/* This is if b7 == 0 */
+	return (bx << 1);
+}
+
+static byte mult (const byte ax, const byte bx)
+{
+	int i;
+	byte xibx = bx;
+	byte ai;
+	byte total = 0x00;
+
+	for (i=0; i<8; ++i)
+	{
+		/* Find a0 through a7 */
+		ai = ax & (1 << i);
+
+		/* If ai is not zero, add it into the total */
+		if (ai)
+			total ^= xibx;
+
+		/* Update x^i * b(x) */
+		xibx = xtimes (xibx);
+	}
+
+	return total;
+}
+
+static void printState (byteArray &bytes, std::string name)
+{
+	int i;
+
+	std::cout << name << ":  ";
+	for (i=0; i<16; ++i)
+		std::printf ("%.2x", bytes.at(i));
+
+	std::printf ("\n");
+}
+
+
+/* vim: set ts=4 sts=4 sw=4 noet tw=112 nowrap: */
diff --git a/aes.hpp b/aes.hpp
new file mode 100644
index 0000000..43808d4
--- /dev/null
+++ b/aes.hpp
@@ -0,0 +1,56 @@
+#ifndef AES_HPP
+#define AES_HPP
+
+#include <vector>
+#include <cstdio>
+#include <string>
+#include <iostream>
+using namespace std;
+
+// Byte and Byte Array types
+typedef unsigned char byte;
+typedef std::vector<byte> byteArray;
+
+// Word and Word Array types
+typedef unsigned int word;
+typedef std::vector<word> wordArray;
+
+// Exceptions
+class incorrectKeySizeException { };
+class incorrectTextSizeException { };
+class badStateArrayException { };
+
+class AES
+{
+	public:
+		AES (const byteArray& key);
+
+		byteArray encrypt (const byteArray& plaintext) const;
+		byteArray decrypt (const byteArray& ciphertext) const;
+
+	private:
+		unsigned int Nb;
+		unsigned int Nk;
+		unsigned int Nr;
+
+		wordArray keySchedule;
+
+		void KeyExpansion (const byteArray& key, wordArray& w) const;
+		word SubWord (const word& input) const;
+		word RotWord (const word& input) const;
+		wordArray GetRoundKey (const int round) const;
+
+		void SubBytes (byteArray& state) const;
+		void ShiftRows (byteArray& state) const;
+		void MixColumns (byteArray& state) const;
+		void AddRoundKey (byteArray& state, const wordArray& w) const;
+
+		void InvSubBytes (byteArray& state) const;
+		void InvShiftRows (byteArray& state) const;
+		void InvMixColumns (byteArray& state) const;
+
+};
+
+#endif /* AES_HPP */
+
+/* vim: set ts=4 sts=4 sw=4 noet tw=112: */
diff --git a/aes_test.cpp b/aes_test.cpp
new file mode 100644
index 0000000..c3a0543
--- /dev/null
+++ b/aes_test.cpp
@@ -0,0 +1,146 @@
+#include "aes.hpp"
+#include <iostream>
+#include <fstream>
+
+static void print_by_16s (const byteArray b, std::string name)
+{
+	int i, j;
+
+	std::cout << "Printing: " << name << "\n" << std::endl;
+
+	for (i=0; i<b.size()/16; ++i)
+	{
+		for (j=0; j<16; ++j)
+			std::printf ("%.2x ", b.at(i*16+j));
+
+		std::printf ("\n");
+	}
+}
+
+
+int main (int argc, char *argv[])
+{
+
+//#define DEBUG 1
+//#define MY_TEST2 1
+//#define PROBLEM2 1
+//#define PROBLEM3 1
+
+#if PROBLEM2
+	int i;
+
+	/* Read the Key */
+	byteArray key;
+
+	for (i=0; i<16; ++i)
+		key.push_back (getchar());
+
+	/* Create the AES object */
+	AES aes (key);
+
+	/* Read all characters of input */
+	byteArray input;
+
+	for (char c=getchar(); c != EOF; c=getchar())
+	{
+		input.push_back (c);
+
+		if (input.size() == 16)
+		{
+			byteArray ciphertext = aes.encrypt (input);
+			byteArray plaintext  = aes.decrypt (ciphertext);
+
+			/* Print it */
+			for (i=0; i<plaintext.size(); ++i)
+				std::printf ("%c", plaintext.at(i));
+
+			input.clear();
+		}
+
+	}
+#endif // PROBLEM2
+
+#if PROBLEM3
+	std::ifstream fin ("Problem3.out", ios::binary);
+	int i;
+
+	byteArray key;
+
+	for (i=0; i<24; ++i)
+		key.push_back (fin.get());
+
+	/* Create the AES Object */
+	AES aes (key);
+
+	byteArray input;
+
+	/* Read the input and decrypt */
+	for (char c=fin.get(); !fin.eof(); c=fin.get())
+	{
+		input.push_back (c);
+
+		if (input.size() == 16)
+		{
+			byteArray plaintext = aes.decrypt (input);
+
+			/* Print it */
+			for (i=0; i<plaintext.size(); ++i)
+				std::printf ("%c", plaintext.at(i));
+
+		input.clear();
+		}
+	}
+#endif //PROBLEM3
+
+#if MY_TEST
+	const int inputSize = 4096;
+
+	byteArray key (16, 2);
+	byteArray input (inputSize, 0);
+	int i, j;
+
+	for (i=0; i<inputSize; ++i)
+		input.at (i) = i;
+
+	AES aes (key);
+
+	print_by_16s (input, "Input");
+
+	byteArray ciphertext = aes.encrypt (input);
+
+	print_by_16s (ciphertext, "Ciphertext");
+
+	byteArray plaintext  = aes.decrypt (ciphertext);
+
+	print_by_16s (plaintext, "Plaintext");
+#endif
+
+#if MY_TEST2
+	int i;
+
+	byteArray key (16, 0);
+
+	for (i=0; i<16; ++i)
+		key.at (i) = i;
+
+	print_by_16s (key, "Key");
+
+	byteArray plaintext (16, 0);
+
+	for (i=0; i<16; ++i)
+		plaintext.at (i) = (i << 4) | i;
+
+	print_by_16s (plaintext, "Plaintext");
+
+	AES aes (key);
+
+	byteArray ciphertext = aes.encrypt (plaintext);
+
+	print_by_16s (ciphertext, "Ciphertext");
+#endif // MY_TEST2
+
+	return 0;
+
+}
+
+/* vim: set ts=4 sts=4 sw=4 noet tw=120 nowrap: */
diff --git a/problem2.cpp b/problem2.cpp
new file mode 100644
index 0000000..83b53c2
--- /dev/null
+++ b/problem2.cpp
@@ -0,0 +1,41 @@
+#include "aes.hpp"
+#include <cstdio>
+
+int main (int argc, char *argv[])
+{
+	int i;
+
+	/* Read the Key */
+	byteArray key;
+
+	for (i=0; i<16; ++i)
+		key.push_back (getchar());
+
+	/* Create the AES object */
+	AES aes (key);
+
+	/* Read all characters of input */
+	byteArray input;
+
+	for (char c=getchar(); c != EOF; c=getchar())
+	{
+		input.push_back (c);
+
+		if (input.size() == 16)
+		{
+			byteArray ciphertext = aes.encrypt (input);
+			byteArray plaintext  = aes.decrypt (ciphertext);
+
+			/* Print it */
+			for (i=0; i<plaintext.size(); ++i)
+				std::printf ("%c", plaintext.at(i));
+
+			input.clear();
+		}
+
+	}
+
+	return 0;
+}
+
+/* vim: set ts=4 sts=4 sw=4 noet tw=120 nowrap: */
diff --git a/problem3.cpp b/problem3.cpp
new file mode 100644
index 0000000..df343bd
--- /dev/null
+++ b/problem3.cpp
@@ -0,0 +1,40 @@
+#include "aes.hpp"
+#include <fstream>
+#include <cstdio>
+
+int main (int argc, char *argv[])
+{
+	std::ifstream fin ("Problem3.out", ios::binary);
+	int i;
+
+	byteArray key;
+
+	for (i=0; i<24; ++i)
+		key.push_back (fin.get());
+
+	/* Create the AES Object */
+	AES aes (key);
+
+	byteArray input;
+
+	/* Read the input and decrypt */
+	for (char c=fin.get(); !fin.eof(); c=fin.get())
+	{
+		input.push_back (c);
+
+		if (input.size() == 16)
+		{
+			byteArray plaintext = aes.decrypt (input);
+
+			/* Print it */
+			for (i=0; i<plaintext.size(); ++i)
+				std::printf ("%c", plaintext.at(i));
+
+		input.clear();
+		}
+	}
+
+	return 0;
+}
+
+/* vim: set ts=4 sts=4 sw=4 noet tw=120 nowrap: */
-- 
2.25.1