#include <iostream>
#include <thread>
#include <chrono>
#include <cmath>
#include <vector>
 
using namespace std;
 
// Function to simulate a layer of the fusion pellet
void printLayer(int layer, int maxLayer) {
    int width = maxLayer*2+1;
    for(int y=0; y<width; y++) {
        for(int x=0; x<width; x++) {
            int dist = abs(x - maxLayer) + abs(y - maxLayer);
            if(dist == layer) cout << "*";          // Fusion halo
            else if(dist < layer) cout << "O";      // Fuel pellet inner
            else cout << " ";                        // Empty space
        }
        cout << endl;
    }
}
 
// Function to print alpha particles as arrows
void printAlphaParticles() {
    vector<string> arrows = { " ^ ", " > ", " v ", " < "};
    cout << "Alpha particles emitted: ";
    for(auto a: arrows) cout << a;
    cout << endl;
}
 
int main() {
    int maxLayer = 5;      // Max fusion halo layers
    int delayMs = 500;     // Delay between steps
 
    cout << "=== Antimatter-Catalyzed Fusion Simulation ===\n";
    cout << "Energy Units: Antimatter (~1.88 GeV), Alpha (~3-8 MeV)\n\n";
 
    for(int layer=0; layer<=maxLayer; layer++) {
        cout << "Step " << layer+1 << ": Fusion halo expanding...\n";
        printLayer(layer, maxLayer);
        printAlphaParticles();
        cout << "\n";
        this_thread::sleep_for(chrono::milliseconds(delayMs));
    }
 
    cout << "Thrust / Energy Output achieved!\n";
    cout << "Antimatter energy released: ~1.88 GeV\n";
    cout << "Alpha particle energy released: ~" << 3 + rand()%6 << " MeV each\n";
    cout << "=== Simulation Complete ===\n";
 
    return 0;
}
 

I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8dGhyZWFkPgojaW5jbHVkZSA8Y2hyb25vPgojaW5jbHVkZSA8Y21hdGg+CiNpbmNsdWRlIDx2ZWN0b3I+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKLy8gRnVuY3Rpb24gdG8gc2ltdWxhdGUgYSBsYXllciBvZiB0aGUgZnVzaW9uIHBlbGxldAp2b2lkIHByaW50TGF5ZXIoaW50IGxheWVyLCBpbnQgbWF4TGF5ZXIpIHsKICAgIGludCB3aWR0aCA9IG1heExheWVyKjIrMTsKICAgIGZvcihpbnQgeT0wOyB5PHdpZHRoOyB5KyspIHsKICAgICAgICBmb3IoaW50IHg9MDsgeDx3aWR0aDsgeCsrKSB7CiAgICAgICAgICAgIGludCBkaXN0ID0gYWJzKHggLSBtYXhMYXllcikgKyBhYnMoeSAtIG1heExheWVyKTsKICAgICAgICAgICAgaWYoZGlzdCA9PSBsYXllcikgY291dCA8PCAiKiI7ICAgICAgICAgIC8vIEZ1c2lvbiBoYWxvCiAgICAgICAgICAgIGVsc2UgaWYoZGlzdCA8IGxheWVyKSBjb3V0IDw8ICJPIjsgICAgICAvLyBGdWVsIHBlbGxldCBpbm5lcgogICAgICAgICAgICBlbHNlIGNvdXQgPDwgIiAiOyAgICAgICAgICAgICAgICAgICAgICAgIC8vIEVtcHR5IHNwYWNlCiAgICAgICAgfQogICAgICAgIGNvdXQgPDwgZW5kbDsKICAgIH0KfQoKLy8gRnVuY3Rpb24gdG8gcHJpbnQgYWxwaGEgcGFydGljbGVzIGFzIGFycm93cwp2b2lkIHByaW50QWxwaGFQYXJ0aWNsZXMoKSB7CiAgICB2ZWN0b3I8c3RyaW5nPiBhcnJvd3MgPSB7ICIgXiAiLCAiID4gIiwgIiB2ICIsICIgPCAifTsKICAgIGNvdXQgPDwgIkFscGhhIHBhcnRpY2xlcyBlbWl0dGVkOiAiOwogICAgZm9yKGF1dG8gYTogYXJyb3dzKSBjb3V0IDw8IGE7CiAgICBjb3V0IDw8IGVuZGw7Cn0KCmludCBtYWluKCkgewogICAgaW50IG1heExheWVyID0gNTsgICAgICAvLyBNYXggZnVzaW9uIGhhbG8gbGF5ZXJzCiAgICBpbnQgZGVsYXlNcyA9IDUwMDsgICAgIC8vIERlbGF5IGJldHdlZW4gc3RlcHMKCiAgICBjb3V0IDw8ICI9PT0gQW50aW1hdHRlci1DYXRhbHl6ZWQgRnVzaW9uIFNpbXVsYXRpb24gPT09XG4iOwogICAgY291dCA8PCAiRW5lcmd5IFVuaXRzOiBBbnRpbWF0dGVyICh+MS44OCBHZVYpLCBBbHBoYSAofjMtOCBNZVYpXG5cbiI7CgogICAgZm9yKGludCBsYXllcj0wOyBsYXllcjw9bWF4TGF5ZXI7IGxheWVyKyspIHsKICAgICAgICBjb3V0IDw8ICJTdGVwICIgPDwgbGF5ZXIrMSA8PCAiOiBGdXNpb24gaGFsbyBleHBhbmRpbmcuLi5cbiI7CiAgICAgICAgcHJpbnRMYXllcihsYXllciwgbWF4TGF5ZXIpOwogICAgICAgIHByaW50QWxwaGFQYXJ0aWNsZXMoKTsKICAgICAgICBjb3V0IDw8ICJcbiI7CiAgICAgICAgdGhpc190aHJlYWQ6OnNsZWVwX2ZvcihjaHJvbm86Om1pbGxpc2Vjb25kcyhkZWxheU1zKSk7CiAgICB9CgogICAgY291dCA8PCAiVGhydXN0IC8gRW5lcmd5IE91dHB1dCBhY2hpZXZlZCFcbiI7CiAgICBjb3V0IDw8ICJBbnRpbWF0dGVyIGVuZXJneSByZWxlYXNlZDogfjEuODggR2VWXG4iOwogICAgY291dCA8PCAiQWxwaGEgcGFydGljbGUgZW5lcmd5IHJlbGVhc2VkOiB+IiA8PCAzICsgcmFuZCgpJTYgPDwgIiBNZVYgZWFjaFxuIjsKICAgIGNvdXQgPDwgIj09PSBTaW11bGF0aW9uIENvbXBsZXRlID09PVxuIjsKCiAgICByZXR1cm4gMDsKfQo=