add picasso writeup

Signed-off-by: Julien CLEMENT <julien.clement@epita.fr>

jujure/content/writeups/fcsc_2023/brachiosaure.md

@@ -1,6 +1,6 @@
 ---
 title: "Inverting a bunch of matrices because reverse engineering or something | Brachiosaure @ FCSC 2023"
-date: "2023-04-30 18:00:00"
+date: "2023-05-01 18:00:00"
 author: "Juju"
 tags: ["Reverse", "Writeup", "fcsc"]
 toc: true
@@ -30,6 +30,9 @@ Author: `\J`
 [brachiosaure](/brachiosaure/brachiosaure)
 
 
+# Writeup
+
+
 ## Overview
 
 As we can guess from the challenge description, this binary will take images
@@ -192,7 +195,7 @@ tricks up my sleeve to solve this, which I finally almost all threw by the
 window when I discovered how easy the solver actually is.
 
 I still think that even though completely over engineered and full of dead ends,
-my though process may be of value for someone I guess.
+my thought process may be of value for someone I guess.
 
 ### Recontextualizing
 
@@ -215,7 +218,7 @@ we may resize the images of the QR codes to add some data.
 
 ### Dead end
 
-So then this is where I had a really dumb idea. I though that maybe I could
+So then this is where I had a really dumb idea. I thought that maybe I could
 append carefully crafted columns and empty lines to the first matrix, and
 carefully crafted line and empty columns to the second matrix so that the result
 of the dot product int overflow back to 0 or 1. Thus outputing 2 new images with the QR code in the top left corner.
@@ -362,13 +365,13 @@ The patched serial png
 {{< image src="/brachiosaure/90QOCSdkzFE3rrYD2GdkrZkh4q_serial.png" style="border-radius: 8px;" >}}
 
 
-As you can notice, there isn't any noise visible by naked eye, this script therefore get the flag 100% of the time.
+As you can notice, there isn't any noise visible by naked eye, this script therefore gets the flag 100% of the time.
 
 Why is it so clean ? Simply because of the strategy we used to compute the
 inverse matrix:
 
 - We add empty matrices: they are black so no noise
 
-- We add identity matrices: they only have the diagonal set to 1 so only a little bit grayer than the black, no noise visible by naked eyes
+- We add identity matrices: they only have the diagonal set to 1 so only a little bit grayer than theso here you go.
 
 - We add the opposite of the matrix, and this is the clean part: our original matrices only hold black and white pixels so respectively `0x0` and `0xff`, so the opposite of `0` is still `0` and the opposite of `0xff` if `1` modulo 256, so like the identity matrix, they are nearly invisible. If you look closely though :eyes: you will see that all white pixels of the QR code were indeed reflected as very faint taint of gray in its inverse matrix on the other image.

jujure/content/writeups/fcsc_2023/picasso.md

@@ -0,0 +1,511 @@
+---
+title: "Reversing and solving nested puzzles | Picasso @ FCSC 2023"
+date: "2023-05-01 18:00:00"
+author: "Juju"
+tags: ["Reverse", "Writeup", "fcsc"]
+toc: true
+---
+
+# Intro
+
+Yes, this is yet another puzzle reversing challenge from `\J` for the FCSC 2023,
+yes I swear I'm also able to flag other types of challenges too, I can't help I
+like algorithms too much. I will try to make lower level content soon :tm:
+:eyes:.
+
+Yes I also know that `face0xff` already did a
+[writeup](https://ctf.0xff.re/2023/fcsc-2023/picasso) on this, but I still want
+to writeup it because I liked the challenge and even though my solution is kind
+of overengineered compared to the one of `face0xff`, I really enjoyed writing
+the solver algorithm for the first puzzle so here you go.
+
+{{< image src="/picasso/meme.jpg" style="border-radius: 8px;" >}}
+
+
+## Challenge description
+
+`reverse` | `482 pts` `8 solves` `:star::star::star:`
+
+```
+Trouvez une entrée permettant d'afficher le message de réussite, et envoyez-la
+sur le service distant pour récupérer le flag.
+
+nc challenges.france-cybersecurity-challenge.fr 2251
+
+SHA256(picasso) =
+f3d6aafce2c069fdf486fce75112ccc09593ddd54452407d4cf066be12daf7fd.
+```
+
+Author: `\J`
+
+## Given files
+
+[picasso](/picasso/picasso)
+
+# Writeup
+
+## Reverse
+
+### Overview
+
+OK, nothing magic you know how this works: it's a puzzle, we need to find the
+input that solves said puzzle. Now the twist is that knowing how `\J` writes his
+challenges, the problem is either NP-complete, obfuscated with 50+ nanomites or
+put together with 12 other puzzles that must be solved together with the same
+input. Make your bet, which one is it this time ?
+
+So I open up my decompiler and I find a single function, the whole thing holds
+in a single function in less than `0x200` bytes of machine code. Let's look this
+up.
+
+{{< code file="/static/picasso/main.c" language="c" >}}
+
+So I cleared everything so you can understand the structure of the code. It is
+really simple:
+
+- First the program copies `0x36` bytes from the init_state array to a local array
+- It then prompts for the password, asking for a maximum of 24 bytes
+- Initialize an index that will be used to iterate over the input in an infinite loop
+- We have 2 clear parts in this loop:
+    - The first one when we are done iterating on the input: We will go back to what this does later but we can clearly see that it performs the final check and prints us the flag if we are good.
+    - The second one is the part that iterates on the input, so the one that is executed first. We can also see that it matches the characters from our input against lower case letters, instantly exiting the program if it isn't one
+
+A first look at the init state shows us that it only contains only numbers
+between 0x0 and 0xf, a good hint that it may be used to index an array of `0x10`
+elements.
+
+### Permutations on the state
+
+Let's start by looking inside that `else` block since it is the one doing stuff
+with our input.
+
+Ok so first thing you need to notice: you may think that characters are matched
+on all lowercase letters but if you look closely, the `'i'` and the `'o'`
+are missing, thus giving us an alphabet of 24 elements.
+
+I know what you are thinking `"But that's also the total size of the input, so
+characters from the input are used to index said input"`, yeah that could have
+been a good idea, but it doesn't do that and it is simply a coincidence.
+
+So anyway we recover the index of the character in the alphabet and divide it
+by 4. The remainder is kept for later and the quotient is used to index
+a permutations matrix, which holds 6 permutation tables of 0x36 elements each.
+
+The permutation matrix is linearized though so I already redefined the dimesions
+so that's why it's showing really nicely in the code below. But to know the
+dimensions I basically saw that the permutation matrix was indexed by `(pos / 4) * 0x36`, since our alphabet holds 24 characters, that gives us `24 / 4 == 6` arrays of size `0x36` in the matrix. We can confirm the 0x36 is we will iterate over `0x36` elements of the permutation array right after this.
+
+We then copy the state and apply the permutations given by the fetched
+permutation array. We repeat this process N times, N being `pos % 4`. We do not
+actually really care about the copy, it is simply there as temporary buffer so
+that permutations don't cancel each other or erase data in the state, it is
+juste an implementation detail.
+
+So to sum this up, each character from our input will chose a permutation array
+to apply to our state and how many times to apply the permutation. They are 6
+permutations table possibles, and we can apply the same permutation between 0
+and 3 times with a single input character.
+
+So basically our input controls how we will shuffle the initial state.
+
+You can also see right bellow the code a sample of the permutations matrix,
+values seem really arbitrary but let's see if you can already get a feeling of
+what this is doing, I will explain it after understanding the entire `if` block
+containing the final check.
+
+{{< code file="/static/picasso/permutations.c" language="c" >}}
+
+And here a sample of the permutation matrix, telling you where each index of the
+state is going to be moved.
+
+{{< code file="/static/picasso/permutations_array.c" language="c" >}}
+
+
+### State checker
+
+OK so the first part shuffled the state, so now I guess this simply checks if
+the resulting shuffled state is valid.
+
+This one may be kind of tricky if you want to understand the implementations
+details, if you don't care you can just skip right when I will recap what the
+puzzle actually is.
+
+So I already renamed everything to keep it clear, I invite you to follow the
+code and its explanations a bit lower in parallel.
+
+
+{{< code file="/static/picasso/slide.c" language="c" >}}
+
+Basically you start with a board represented as an `uint64_t` intialized at
+`0x3da8e0915f2c4b67`. For now, view this board as an array of 16 elements, each element being 4 bits (`4 * 16 = 64`).
+
+As you know 4 bits is a single hexadecimal digit, so actually each hexadecimal
+digit from this board is an element, so the initial board is:
+`[3, d, a, 8, e, 0, 9, 1, 5, f, 2, c, 4, b, 6, 7]`
+
+You may have noticed that every element is unique, they represent the entire
+range of numbers between 0x0 and 0xf, this is no coincidence.
+
+We then create some variables that we will ignore for now, and an index (`i`)
+that will be used to iterate over the state that was shuffled in the first part.
+
+While iterating on the state, we fetch the current byte from the state, extended
+to 64 bits. Remember how we told that state only held values between `0` and
+`0xf` at the beginning ? Good, you are starting to see a pattern.
+
+We create a pointer to a variable I called `allowed_moves`, it is once again a
+two dimensionnal array because we iterate over it in a nested loop. I knew its
+dimensions because this is the part to go to the next sub-array:
+
+```c
+allowed_moves += 0x28;
+```
+
+So I knew each sub-array is `0x28` bytes, and I saw that the index used to
+iterrate the outer loop is the `shifter` variable, initialized to `0x3c` and
+decremented by `4` each loop until it is strictly negative, so a total of `0x10`
+iterations. We then have a matrix containing `0x10` arrays of `0x28`. 
+
+I also knew that the `0x28` bytes corresponded to 5 `uint64_t` since the inner
+loop iterate over an `uint64_t *`
+
+I pasted the entire content of this matrix below so you can have a look, you
+will see that every element of this matrix has a single bit set to 1, and bits that are all sets at a position of `0 % 4` for that matter.
+
+Okay let's take a look at the `shifter`. It is initialized to `0x3c` and is used
+to shift the `board`, before keeping only the lowest weight hex-digit with the
+`& 0xf`. The hex digit is then compared against the byte we fetched from our
+state.  If they differ, we basically skip this loop, decrementing the `shifter`
+by 4 (meaning that on the next iteration we will select the next hex digit of
+the board), and go to the next array in the `allowed_moves` matrix.
+
+So to recap this, if you consider the `board` as an array of 4 bits elements,
+`shifter` is actually an index used to iterate over this array. So we are
+iterating to find at which index is the byte from the `state` inside the
+`board`.  `allowed_moves` will follow the same iterations, meaning that we will
+select an `allowed_moves` sub-array based on where the `state` byte is located
+in the `board`.
+
+Good, so once we find out `state` byte inside `board` we start iterating over
+the `allowed_moves`, basically the move is multiplied by `0xf`, which will effectively set the 4 bits corresponding to the single bit in the initial move, this will now be used as a mask on the board.
+
+So this move simply selects an hex digit from the board actually, if the
+selected byte is NOT 0, then we continue the loop, skipping the move, but if it
+is 0, then we exit the loop. Final check to see if the move itself was not 0, if
+it is then we go back to trying the next `shifter` and `allowed_moves`, if not
+then gg we selected a move. If at the end, no move was selected we go out of the
+loop and basically lose.
+
+So whatever this means, a valid move from a certain position in the board, must
+match with the `0` digit from the `board` being at certain indexes, which, for
+now, since we did dig up the `allowed_moves` matrix yet, feel kind of arbitrary.
+
+Feeling confused ? Don't worry it will soon make sense.
+
+{{< code file="/static/picasso/allowed_moves.c" language="c" >}}
+
+
+Great, we now have a `move` which is the index of the `0` element of the board,
+a `shifter`, being the index of the byte from our `state` in the same board. The
+previous paragraph was meant to make sure this move/shifter association is
+valid.
+
+Let's talk about the `swapper` now, it is composed of 2 parts:
+
+- `c << shifter`: so it will be the representation of the board as if only `c`
+was on it
+
+- `c * move`: the representation of the board as if `c` was at the position
+indicated by the move, and no other element on the board. So it the board, with
+c at the position of the 0 and nothing else on the board
+
+We then compute the next board by `xoring` the current one with the swapper.
+This will have the effect to `xor` the current element of the `state` AND the
+current `0` from the board with the current element of the `state`. This will
+have the effect to actually swap the `0` and our element of the `state`.
+
+And after that we simply go the next `state` byte, until we have done all `0x36`
+of them.
+
+Let's recap:
+
+- each character from the `state` selects the corresponding hex digit in the
+board.
+
+- According to its position in the board, we will try to find a valid move,
+being a move that have `0` as destination, but only certains positions of `0`
+are allowed according to were the hex digit was in the board.
+
+- If a valid move was found, then we swap the `0` and the `state` byte.
+
+FINALLY, there is the final check, after all the moves from the state have been applied, the board must be `0x123456789abcdf0`.
+
+Great, now let's look at the `allowed_moves`.
+
+So what I did to understand what all those valid moves were was to take a pen,
+drew an array of 16 elements, and for each element I drew an arrow to every
+valid moves from that position. So just ignore that there has to be a `0` at
+these positions for now, it's simply a matter of "if there was a 0 there, could
+I move here ?". And it started looking like this for the first few elements:
+
+{{< image src="/picasso/swaps.png" style="border-radius: 8px;" >}}
+
+What you need to notice is that every move can be inverted, so you can go back to your position after moving. MOST tiles can go to the tiles right next to them or to the tiles that are 4 tiles farther. For example, tile at index 1 can go either at 0, 2, or 5. You are probably starting to understand, let's just put
+the final touch: 3 cannot go to 4 and vise-versa.
+
+So this is simply a flattened 4 * 4 grid, where you can only move to the
+adjacent tiles. Let's now replace the indexes I put in the diagram above with
+the actual values of the initial board (`0x3da8e0915f2c4b67`) we have:
+
+{{< image src="/picasso/starting_grid.png" style="border-radius: 8px;" >}}
+
+I know you understood everything but let's just see the target board
+(`0x123456789abcdef0`):
+
+{{< image src="/picasso/target_grid.png" style="border-radius: 8px;" >}}
+
+So we need to get from one to the other in `0x36` moves exactly, each move swapping the 0 with an adjacent tile.
+
+Cool, all of this simply for a slide puzzle. Each element of the state corresponds to the tile that we must push on the `0` (which represent the empty tile in the slide puzzle). Let's go and solve th...
+
+
+### But wait there was a first step
+
+I hope you did not forgot that the state used to select which tiles are moved in
+which order was shuffled by the first step according to our input.
+
+Maybe it's time to understand how this shuffle work.
+
+OK so I got this one super fast actually, let's just get you the context back:
+
+- 6 permutations tables possibles
+- each can be applied `[0;4[` times each time
+- `0x36 == 54` elements in the state
+
+If you are a puzzle fan you already what this is I don't even need to show you
+the permutation matrix again but hey:
+
+```c
+char permutations[0x6][0x36] = 
+{
+    [0x0] = 
+    {
+        [0x00] =  0x09
+        [0x01] =  0x01
+        [0x02] =  0x02
+        [0x03] =  0x0c
+        [0x04] =  0x04
+        [0x05] =  0x05
+        [0x06] =  0x0f
+        [0x07] =  0x07
+        [0x08] =  0x08
+        [0x09] =  0x2d
+        ...
+        [0x34] =  0x34
+        [0x35] =  0x35
+    }
+    [0x1] = 
+    {
+        [0x00] =  0x00
+        [0x01] =  0x01
+        ...
+```
+
+Literally what I did is I saw the first permutation table, noticed that only 1
+element out of 3 were moving (so I though this was a 3 * 3 array) and that the
+first ones were moving 9 tiles further than there original location (so probaly
+a 3 * 3 * 3 cube), oh wait did I say cube ? 6 faces ? can be rotated up to 3 times ? 54 elements ? yeah ok you got it aswell.
+
+It's a rubiks cube.
+
+Each permutation table is a counter-clockwise rotation of a face of the cube.
+
+The process of mapping the cube in its flattened representation and to know
+which permutations corresponds to which face is just teadious so I will juste
+give you this one out.
+
+Do keep in mind that I was not familiar with standard rubiks cube notations, so
+I named a move by the color of the face that is rotating instead of the usual
+U,F,... it felt more intuitive to me.
+
+
+Here are the indexes of each cube tile on the flatenned state:
+
+{{< image src="/picasso/cube_indexes.png" style="border-radius: 8px;" >}}
+
+I attributed colors to faces arbitrarily to respect the usual cube
+configuration, but as long as I keep everythin coherent with this representation
+I will do just fine.
+
+With the above cube, fold it back and see how the first permutation array indeed performs a counter-clockwise rotation of the red face.
+
+Here is the initial state of the cube:
+
+{{< image src="/picasso/cube_init.png" style="border-radius: 8px;" >}}
+
+
+## Solving the slide puzzle
+
+### Recap everything
+
+So know we are done with the reversing part, with our input, the program will:
+
+- Consider each character as a move on a rubiks cube, applying said move
+- Consider the flatenned shuffled cube as a list of slide puzzle moves
+
+Basically we need to solve a rubik's cube in order to ordonate a slide puzzle
+solution.
+
+To solve this, we will need to take one problem at a time, first find a solutions to the slide puzzle, then map the solution to a cube which we will consider as the solved cube, then map the color of the solved cube back to
+the starting cube in order to solve it.
+
+### Naive backtracking
+
+So to start I want to say that I am really impressed that the solver of
+`face0xff` found a solution this fast, I was really convinced that without the
+little tricks that I will show you, it would take too long to search for a 54
+moves solution, but hey seems like it works if you implement every slide puzzle
+heuristics you may think of.
+
+So my approach was a backtracking algorithm, however you can clearly see why
+backtracking 54 moves is not a good solution so we will need to add some
+heuristics to narrow down the search for valid moves.
+
+Basically the algorithm works as follow given a current state of the grid:
+
+- If the puzzle is solved it's won
+- If you have no move left it's lost
+- generate all possible moves
+- try each move by applying the move and recursively calling this procedure
+- if the recursive call returned true then it's won
+- else cancel the move and try the next
+- if no move returned true there is no solution
+
+
+### Slide puzzle heuristics
+
+Okay this algorithm will take way too long, we need to cancel some search paths
+early so that we will not further go down the calls if we can already know that we cannot solve anymore from a given grid.
+
+
+#### Cannot cancel previous move
+
+This heuristic takes the assumption that you cannot play the same move as the
+previous one, because otherwise you will go back to same state as the previous
+move, effectively doing 2 moves for nothing.
+
+Okay this heuristic is not exactly 100% reliable, since actually the puzzle is
+also solvable this way, simply with useless moves. But I was convinced that 54
+was the optimal solution, I don't imagine that `\J` would have asked us to solve
+a puzzle in a non-optimal way.
+
+#### Last move
+
+Since the tile in bottom right needs to be `0`, this means that the last move needs to be `f` or `c`, so if at any point, both of these tiles do not have anymore moves remaining moves while the puzzle isn't already solved then you can already tell it's unsolvable.
+
+Oh I didn't talk about the limited number of moves per tiles ? Don't worry I'll
+explain in a minute.
+
+### Rubik's cube heuristics
+
+These are the heuristics you can determine based on the fact that the moves
+are shuffled from a rubik's cube initial state.
+
+#### Face centers
+
+In a rubik's cube the center tile of each face CANNOT move, this means that since we know the starting position of the cube, we already know that:
+
+- Move `4` will be `f`
+- Move `d` will be `2`
+- Move `16` will be `c`
+- Move `1f` will be `a`
+- Move `28` will be `a`
+- Move `31` will be `4`
+
+#### Move counts and manhattan distance
+
+Of course I kept the best one for the end, this is by far the most effective
+one and it is probably enough to solve the whole thing by itself.
+
+This heuristic is based on the fact that we can only shuffle the cube, we can
+never add or remove moves, therefore we can know how many times a tile will be
+moved. For example there are 4 `1` on the cube, so this means that the `1` tile
+of the slide puzzle must be played exactly 4 times, no more no less.
+
+So now, knowing this, while backtracking, we can keep track of how many times a tile has moved by decrementing its move counter.
+
+Then, we can compute the manhattan distance of the tile to its supposed
+location, which is the minimum number of moves that would be required for the
+tile to go to that location if there were no other tile on the grid. For example
+tile `f` placed in the top left corner will have a manhattan distance of 5 to
+its supposed location.
+
+So whenever the move counter of a tile becomes lower than its manhattant distance to its suppose location, you know that you will never be able to place the tile back to where it is supposed to be.
+
+
+### Solver
+
+So I implemented this in rust because I wanted this to be fast:
+
+{{< code file="/static/picasso/slide_solver.rs" language="rust" >}}
+
+And guess what ? It was super duper fast:
+
+```console
+$ time cargo run --release
+    Finished release [optimized] target(s) in 0.00s
+     Running `target/release/picasso`
+[e, 5, 4, b, f, e, d, 3, 5, 4, e, d, 9, 2, 6, f, b, e, d, 9, 2, 1, c, 7, f, b, e, d, 9, 2, 1, a, 3, 1, 4, 5, 1, 4, 2, 6, a, 3, 4, 2, 6, a, 7, c, 8, 4, 3, 7, b, f]
+
+________________________________________________________
+Executed in   77.37 millis    fish           external
+   usr time   71.83 millis  182.00 micros   71.65 millis
+   sys time    3.93 millis  673.00 micros    3.26 millis
+```
+
+That shit is probably completely overengineered and I don't care, I had fun
+writing it.
+
+I checked manually that this solution was correct and moved to the next step.
+
+
+## Solving the rubik's cube
+
+We now have our solution for the slide puzzle, so let's map it as a solved
+rubik's cube, we know the mapping using the index cube net which gives us which
+tiles of the cube corresponds to which move index on the slide puzzle:
+
+{{< image src="/picasso/solved_cube.png" style="border-radius: 8px;" >}}
+
+Since this is the solved cube, we attribute all the colors of the tiles.
+
+We now need to attribute the corresponding colors to the initial state of the
+cube so that we can input a valid cube to a solver.
+
+To do this, we use the fact that tiles that are connected along an edge or a
+corner will always be linked, so this way we can identify the tile groups on the
+cube and identify their color from the solved cube.
+
+{{< image src="/picasso/cube_init_colors.png" style="border-radius: 8px;" >}}
+
+Then I basically just used [this online
+solver](https://www.grubiks.com/solvers/rubiks-cube-3x3x3/), put manually all
+the colors, clicked on solve, and the site instantly found a solution in 22
+moves. So I noted manually every move, which face was supposed to move and wrote
+a script that mapped the moves back to the expected output of the program.
+
+In the script `('g', 2)` means rotate the green face 2 times, `('o', 3)` means rotate the orange face 3 times counter-clock wise ...
+
+{{< code file="/static/picasso/solve.py" language="python" >}}
+
+And piping this script in the program indeed gets us the flag.
+
+```console
+$ ./solve.py | ./picasso
+Password: Win!
+Send your input on the remote service to retrieve the flag.
+
+$ ./solve.py | nc challenges.france-cybersecurity-challenge.fr 2251
+Password: Win!
+FCSC{235b605a121bdd4b09adc4823bdf0967c446647c1ec69234813068a916fd83a6}
+```

jujure/static/picasso/allowed_moves.c

@@ -0,0 +1,131 @@
+uint64_t allowed_moves[0x10][0x5] = 
+{
+    [0x0] = 
+    {
+        [0x0] =  0x0100000000000000
+        [0x1] =  0x0000100000000000
+        [0x2] =  0x0000000000000000
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0x1] = 
+    {
+        [0x0] =  0x1000000000000000
+        [0x1] =  0x0010000000000000
+        [0x2] =  0x0000010000000000
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0x2] = 
+    {
+        [0x0] =  0x0100000000000000
+        [0x1] =  0x0001000000000000
+        [0x2] =  0x0000001000000000
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0x3] = 
+    {
+        [0x0] =  0x0010000000000000
+        [0x1] =  0x0000000100000000
+        [0x2] =  0x0000000000000000
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0x4] = 
+    {
+        [0x0] =  0x1000000000000000
+        [0x1] =  0x0000010000000000
+        [0x2] =  0x0000000010000000
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0x5] = 
+    {
+        [0x0] =  0x0100000000000000
+        [0x1] =  0x0000100000000000
+        [0x2] =  0x0000001000000000
+        [0x3] =  0x0000000001000000
+        [0x4] =  0x0000000000000000
+    }
+    [0x6] = 
+    {
+        [0x0] =  0x0010000000000000
+        [0x1] =  0x0000010000000000
+        [0x2] =  0x0000000100000000
+        [0x3] =  0x0000000000100000
+        [0x4] =  0x0000000000000000
+    }
+    [0x7] = 
+    {
+        [0x0] =  0x0001000000000000
+        [0x1] =  0x0000001000000000
+        [0x2] =  0x0000000000010000
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0x8] = 
+    {
+        [0x0] =  0x0000100000000000
+        [0x1] =  0x0000000001000000
+        [0x2] =  0x0000000000001000
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0x9] = 
+    {
+        [0x0] =  0x0000010000000000
+        [0x1] =  0x0000000010000000
+        [0x2] =  0x0000000000100000
+        [0x3] =  0x0000000000000100
+        [0x4] =  0x0000000000000000
+    }
+    [0xa] = 
+    {
+        [0x0] =  0x0000001000000000
+        [0x1] =  0x0000000001000000
+        [0x2] =  0x0000000000010000
+        [0x3] =  0x0000000000000010
+        [0x4] =  0x0000000000000000
+    }
+    [0xb] = 
+    {
+        [0x0] =  0x0000000100000000
+        [0x1] =  0x0000000000100000
+        [0x2] =  0x0000000000000001
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0xc] = 
+    {
+        [0x0] =  0x0000000010000000
+        [0x1] =  0x0000000000000100
+        [0x2] =  0x0000000000000000
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0xd] = 
+    {
+        [0x0] =  0x0000000001000000
+        [0x1] =  0x0000000000001000
+        [0x2] =  0x0000000000000010
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0xe] = 
+    {
+        [0x0] =  0x0000000000100000
+        [0x1] =  0x0000000000000100
+        [0x2] =  0x0000000000000001
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+    [0xf] = 
+    {
+        [0x0] =  0x0000000000010000
+        [0x1] =  0x0000000000000010
+        [0x2] =  0x0000000000000000
+        [0x3] =  0x0000000000000000
+        [0x4] =  0x0000000000000000
+    }
+}

jujure/static/picasso/main.c

@@ -0,0 +1,71 @@
+char init_state[0x36] = 
+{
+	0x0d, 0x03, 0x02, 0x0e, 0x0f, 0x09, 0x07, 0x0e, 0x0d, 0x06, 0x03, 0x02, 0x07, 0x02, 0x08, 0x0e,
+	0x07, 0x04, 0x04, 0x03, 0x0d, 0x04, 0x0c, 0x03, 0x0f, 0x02, 0x0a, 0x05, 0x09, 0x01, 0x06, 0x0a,
+	0x0b, 0x02, 0x05, 0x0c, 0x0e, 0x0b, 0x0d, 0x01, 0x0a, 0x01, 0x05, 0x09, 0x01, 0x0f, 0x06, 0x0e,
+	0x04, 0x04, 0x0b, 0x04, 0x0b, 0x0f, 0x00
+};
+
+
+int32_t main(int32_t argc, char** argv, char** envp)
+{
+    char* const init_state_ptr = &init_state;
+    char state[0x36];
+    char* state_ptr = &state;
+    for (int64_t size = 0x36; size != 0; size -= 1)
+    {
+        *state_ptr = *init_state_ptr;
+        state_ptr += 1;
+        init_state_ptr += 1;
+    }
+
+    printf("Password: ", init_state_ptr);
+    fflush(stdout);
+
+    char input[0x18];
+    scanf("%24s", &input);
+
+    int64_t state_index = 0;
+
+    while (true)
+    {
+        if (strlen(&input) <= state_index)
+        {
+            // Cutting everything that happens
+            ...
+            if (check == OK)
+            {
+                puts("Win!");
+                FILE file = fopen("flag.txt", &rb);
+                if (file != NULL)
+                {
+                    char flag[0x46];
+                    fread(&flag, 1, 0x46, file);
+                    fclose(file);
+                    puts(flag);
+                    return 0;
+                }
+                else
+                    puts("Send your input on the remote se…");
+            }
+            else
+                puts("Nope!");
+            break;
+        }
+        else
+        {
+            char* pos_ptr = strchr("abcdefghjklmnpqrstuvwxyz", input[state_index]);
+            if (pos_ptr != NULL)
+            {
+                // Do stuff iterating on the state
+                ...
+                state_index += 1;
+                continue;
+            }
+        }
+        puts("Nope!");
+        break;
+    }
+    exit(1);
+    /* no return */
+}

jujure/static/picasso/permutations.c

@@ -0,0 +1,37 @@
+else
+{
+    char* pos_ptr = strchr("abcdefghjklmnpqrstuvwxyz", input[state_index]);
+    if (pos_ptr != NULL)
+    {
+        int64_t pos = pos_ptr - "abcdefghjklmnpqrstuvwxyz";
+        int32_t remainder = pos % 4;
+        char* permutation_array = &permutations[pos / 4];
+        while (true)
+        {
+            if (remainder < 1)
+                break;
+            remainder -= 1;
+
+            char state_copy[0x36];
+            char* state_copy_ptr = &state_copy;
+
+            char* state_ptr_1 = &state;
+
+            int64_t i = 0x36;
+            for (; i != 0; i -= 1)
+            {
+                *state_copy_ptr = *state_ptr_1;
+                state_copy_ptr += 1;
+                state_ptr_1 += 1;
+            }
+            do
+            {
+                char elem = state_copy[permutation_array[i]];
+                state[i] = elem;
+                i += 1
+            } while (i != 0x36);
+        }
+        state_index += 1;
+        continue;
+    }
+}

jujure/static/picasso/permutations_array.c

@@ -0,0 +1,56 @@
+char permutations[0x6][0x36] = 
+{
+    [0x0] = 
+    {
+        [0x00] =  0x09
+        [0x01] =  0x01
+        [0x02] =  0x02
+        [0x03] =  0x0c
+        [0x04] =  0x04
+        [0x05] =  0x05
+        [0x06] =  0x0f
+        [0x07] =  0x07
+        [0x08] =  0x08
+        [0x09] =  0x2d
+        ...
+        [0x34] =  0x34
+        [0x35] =  0x35
+    }
+    [0x1] = 
+    {
+        [0x00] =  0x00
+        [0x01] =  0x01
+        [0x02] =  0x2a
+        [0x03] =  0x03
+        ...
+        [0x34] =  0x34
+        [0x35] =  0x11
+    [0x2] = 
+    {
+        [0x00] =  0x02
+        [0x01] =  0x05
+        ...
+        [0x35] =  0x35
+    }
+    [0x3] = 
+    {
+        [0x00] =  0x00
+        [0x01] =  0x01
+        ...
+        [0x35] =  0x33
+    }
+    [0x4] = 
+    {
+        [0x00] =  0x00
+        [0x01] =  0x01
+        ...
+        [0x35] =  0x35
+    }
+    [0x5] = 
+    {
+        [0x00] =  0x21
+        [0x01] =  0x1e
+        ...
+        [0x35] =  0x14
+    }
+}

jujure/static/picasso/slide.c

@@ -0,0 +1,67 @@
+if (strlen(&input) <= state_index)
+{
+    int64_t board = 0x3da8e0915f2c4b67;
+    int64_t next_board;
+    int64_t i = 0;
+    uint64_t move;
+    while (true)
+    {
+        uint64_t c = state[i];
+        void* allowed_moves = &allowed_moves;
+        int32_t shifter = 0x3c;
+        do
+        {
+            if (c == ((board >> shifter) & 0xf))
+            {
+                uint64_t* moves = allowed_moves;
+                do
+                {
+                    move = *moves;
+
+                    if (move == 0)
+                        break;
+
+                    moves += 1;
+                } while (((move * 0xf) & board) != 0);
+
+                if (move != 0)
+                    break;
+            }
+            shifter -= 4;
+            allowed_moves += 0x28;
+        } while (shifter != 0xfffffffc);
+
+        if (c != ((board >> shifter) & 0xf) || move == 0)
+            break;
+
+        int64_t swapper = ((c * move) ^ (c << shifter));
+        next_board = (board ^ swapper);
+        if (board == swapper)
+            break;
+
+        i += 1;
+
+        if (i == 0x36)
+            break;
+
+        board = next_board;
+    }
+    if (next_board == 0x123456789abcdef0)
+    {
+        puts("Win!");
+        FILE file = fopen("flag.txt", &rb);
+        if (file != NULL)
+        {
+            char flag[0x46];
+            fread(&flag, 1, 0x46, file);
+            fclose(file);
+            puts(flag);
+            return 0;
+        }
+        else
+            puts("Send your input on the remote se…");
+    }
+    else
+        puts("Nope!");
+    break;
+}

jujure/static/picasso/slide_solver.rs

@@ -0,0 +1,193 @@
+use std::collections::HashMap;
+
+struct Board {
+    board: [usize; 0x10],
+    moves: Vec<usize>,
+    max_moves: usize,
+    move_number: usize,
+
+    empt_idx: usize,
+
+
+    solved_grid: [usize; 0x10],
+
+    fixed_moves: HashMap<usize, usize>,
+
+    move_limit: [isize; 0x10]
+}
+
+impl Board {
+    pub fn new(board: [usize; 0x10]) -> Self {
+        let empt_idx: usize = board.iter().position(|&r| r == 0).unwrap().try_into().unwrap();
+        let all_moves: [usize; 0x36] = [
+            0x0d, 0x03, 0x02, 0x0e, 0x0f, 0x09, 0x07, 0x0e, 0x0d, 0x06, 0x03, 0x02, 0x07, 0x02, 0x08, 0x0e,
+            0x07, 0x04, 0x04, 0x03, 0x0d, 0x04, 0x0c, 0x03, 0x0f, 0x02, 0x0a, 0x05, 0x09, 0x01, 0x06, 0x0a,
+            0x0b, 0x02, 0x05, 0x0c, 0x0e, 0x0b, 0x0d, 0x01, 0x0a, 0x01, 0x05, 0x09, 0x01, 0x0f, 0x06, 0x0e,
+            0x04, 0x04, 0x0b, 0x04, 0x0b, 0x0f
+        ];
+        let solved_grid: [usize; 0x10] = [1,2,3,4,5,6,7,8,9,0xa,0xb,0xc,0xd,0xe,0xf,0];
+
+        let mut fixed_moves: HashMap<usize, usize> = HashMap::new();
+
+        for i in vec![4, 0xd, 0x16, 0x1f, 0x28, 0x31] {
+            fixed_moves.insert(i, all_moves[i]);
+        }
+
+        let mut move_limit: [isize; 0x10] = [0; 0x10];
+
+        for mv in all_moves.iter() {
+            move_limit[*mv] += 1;
+        }
+
+        Board {
+            board: board,
+            moves: Vec::new(),
+            max_moves: 0x36,
+            move_number: 0,
+            empt_idx: empt_idx,
+            solved_grid: solved_grid,
+            fixed_moves: fixed_moves,
+            move_limit: move_limit,
+        }
+    }
+
+    pub fn solved(&self) -> bool {
+        self.board == self.solved_grid
+    }
+
+    pub fn get_move(&self, off: Vec<isize>) -> Option<usize> {
+        let empt_idx: isize = self.empt_idx.try_into().unwrap();
+        let empt_i: isize = empt_idx / 4;
+        let empt_j: isize = empt_idx % 4;
+        let src_i: isize = empt_i + off[0];
+        let src_j: isize = empt_j + off[1];
+        if src_i < 4 && src_i >= 0 && src_j < 4 && src_j >= 0 {
+            let idx: usize = (src_i * 4 + src_j).try_into().unwrap();
+            return Some(self.board[idx])
+        }
+        None
+    }
+
+    pub fn generate_moves(&self) -> Vec<usize> {
+        let mut moves: Vec<usize> = Vec::new();
+
+        for off in vec![vec![-1, 0], vec![1, 0], vec![0, -1], vec![0, 1]] {
+            let maybe_move = self.get_move(off);
+            if let Some(mv) = maybe_move {
+                moves.push(mv);
+            }
+        }
+
+        moves
+    }
+
+    pub fn apply_move(&mut self, mv: usize, undo: bool) {
+        let i: usize = self.board.iter().position(|&x| x == mv).unwrap();
+        let tmp: usize = self.board[i];
+        self.board[i] = self.board[self.empt_idx];
+        self.board[self.empt_idx] = tmp;
+        self.empt_idx = i;
+        if undo {
+            self.move_number -= 1;
+            self.moves.pop();
+            self.move_limit[mv] += 1;
+        }
+        else {
+            self.move_number += 1;
+            self.moves.push(mv);
+            self.move_limit[mv] -= 1;
+        }
+    }
+
+    pub fn manhattan_distance(&self, mv: usize) -> isize {
+        let mv_idx: isize = self.board.iter().position(|&x| x == mv).unwrap().try_into().unwrap();
+        let mv_i: isize = mv_idx / 4;
+        let mv_j: isize = mv_idx % 4;
+
+        let supposed_idx: isize = (mv - 1).try_into().unwrap();
+        let supposed_i: isize = supposed_idx / 4;
+        let supposed_j: isize = supposed_idx % 4;
+
+        let manhattan: isize = (supposed_j - mv_j).abs() + (supposed_i - mv_i).abs();
+
+        manhattan
+    }
+
+    pub fn solve(&mut self) -> bool {
+
+        // GG
+        if self.solved() {
+            return true;
+        }
+        
+        // No more moves
+        if self.move_number >= self.max_moves {
+            return false;
+        }
+
+        // Last move need to be 0xc or 0xf
+        if self.move_limit[0xc] <= 0 && self.move_limit[0xf] <= 0 {
+            return false;
+        }
+
+        let last_move: Option<usize> = self.moves.last().copied();
+
+        if let Some(mv) = last_move {
+            // Tile cannot possibly go back to its supposed location
+            let man: isize = self.manhattan_distance(mv);
+            if self.move_limit[mv] < man {
+                return false
+            }
+        }
+
+        // Generate possible moves from current grid
+        let mut moves: Vec<usize> = self.generate_moves();
+
+        // Cannot apply same move twice
+        if let Some(mv) = last_move {
+            moves = moves.iter().filter(|&x| *x != mv).cloned().collect();
+        }
+
+        // Filter fixed moves from the cube center tiles
+        if self.fixed_moves.contains_key(&self.move_number) {
+            let mv: usize = *self.fixed_moves.get(&self.move_number).unwrap();
+            moves = moves.iter().filter(|&x| *x == mv).cloned().collect();
+        }
+
+
+        for mv in moves {
+            // Apply move
+            self.apply_move(mv, false);
+            // Recursive call
+            let solved: bool = self.solve();
+            // GG
+            if solved {
+                return true;
+            }
+            // Undo move
+            self.apply_move(mv, true);
+        }
+
+        false
+    }
+}
+
+
+fn solve() {
+    let arr: [usize; 0x10] = [
+        3, 0xd, 0xa, 8,
+        0xe, 0, 9, 1,
+        5, 0xf, 2, 0xc,
+        4, 0xb, 6, 7
+    ];
+
+    let mut board: Board = Board::new(arr);
+
+    board.solve();
+
+    println!("{:x?}", board.moves);
+}
+
+fn main() {
+    solve()
+}

jujure/static/picasso/solve.py

@@ -0,0 +1,55 @@
+#!/usr/bin/env python3
+
+import os
+
+alphabet = b"abcdefghjklmnpqrstuvwxyz"
+
+moves = [
+    ('g', 2),
+    ('o', 3),
+    ('g', 3),
+    ('b', 2),
+    ('y', 2),
+    ('g', 3),
+    ('w', 1),
+    ('o', 3),
+    ('r', 3),
+    ('y', 1),
+    ('r', 3),
+    ('w', 1),
+    ('b', 1),
+    ('o', 2),
+    ('g', 1),
+    ('y', 2),
+    ('g', 1),
+    ('y', 2),
+    ('r', 2),
+    ('g', 2),
+    ('w', 2),
+    ('o', 2),
+]
+
+sequence = b""
+
+def face_id(face):
+    res = 0
+    if face == 'w':
+        res = 2
+    if face == 'b':
+        res = 4
+    if face == 'o':
+        res  = 1
+    if face == 'r':
+        res = 0
+    if face == 'y':
+        res = 3
+    if face == 'g':
+        res = 5
+
+    return res
+
+for move in moves:
+    i = face_id(move[0]) * 4 + move[1]
+    sequence += alphabet[i].to_bytes(1, "little")
+
+os.write(1, sequence + b'\n')