Клеточный автомат блуждания частицы по поверхности фуллерена C60

class TForm1 : public TForm
{
__published: // IDE-managed Components
 
    TPanel * shapes_pnl;
    TButton * btn_step;
 
    void __fastcall btn_stepClick( TObject * Sender );
 
private: // User declarations
 
        public : // User declarations
 
    void __fastcall cell_mouse_down( TObject * Sender, TMouseButton Button, TShiftState Shift, int X, int Y );
    __fastcall TForm1( TComponent * Owner );
    __fastcall ~TForm1( );
} ;

#include <vcl.h>
#pragma hdrstop
 
#include <list>
#include <vector>
#include <map>
#include <math.h>
#include <algorithm>
#include <StrUtils.hpp>
#include <Windows.h>
#include <Math.hpp>
 
#include <ShellAnimations.hpp>
 
#include "Unit1.h"
// ---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 * Form1;
 
class t_cell_shape : public TShape
{
    void __fastcall Paint( )
    {
        Brush->Color = charged ? ( gp == 3 ? clAqua : clSkyBlue ) : ( rp == 5 ? clGreen : clRed );
        TShape::Paint( );
        Canvas->TextOut( Width / 2 - 15, 5, Name );
        Canvas->TextOut( Width / 2 - 30, 20, "rp: " + IntToStr( rp ) + "; gp: " + IntToStr( gp ) );
 
        String _neibs;
        for ( t_neibs::iterator it( neibs.begin( ) ), end( neibs.end( ) ); it != end; ++it )
        {
            _neibs += IntToStr( it->first ) + " ";
        }
        Canvas->TextOut( 5, 35, "neibs: " + _neibs );
        Canvas->TextOut( 10, 50, "charge_id: " + IntToStr( charge_id ) );
    }
 
public:
    bool charged, just_received;
    int rp, gp; // receive_priority, give_priority
    int charge_id;
 
    typedef std::map < int, t_cell_shape * > t_neibs;
 
    t_neibs neibs;
 
    void __fastcall move_charge( t_cell_shape * sh )
    {
        sh->charged = true;
        sh->charge_id = charge_id;
        sh->just_received = true;
        charged = false;
        charge_id = 0;
    }
 
    __fastcall t_cell_shape( TComponent * Owner ) :
        TShape( Owner ),
        charge_id( 0 ),
        just_received( false )
    {
    }
} ;
 
t_cell_shape * shapes[ 12 ];
 
// ---------------------------------------------------------------------------
__fastcall TForm1::TForm1( TComponent * Owner ) :
    TForm( Owner )
{
    t_cell_shape * sh;
    int w( shapes_pnl->Width / 4 - 30 );
 
    std::srand( std::time( 0 ) );
 
    for ( int i( 0 ), i_limit( 12 ); i < i_limit; ++i )
    {
        sh = shapes[ i ] = new t_cell_shape( this );
 
        sh->Parent = shapes_pnl;
        sh->Name = "cell_" + IntToStr( i );
        sh->Shape = stCircle;
        sh->Width = w;
        sh->Height = w;
        sh->Pen->Style = psClear;
        sh->OnMouseDown = &cell_mouse_down;
    }
    for ( size_t i( 0 ), i_limit( 3 ); i < i_limit; ++i )
    {
        for ( size_t k( 0 ), k_limit( 4 ); k < k_limit; ++k )
        {
            sh = ( *( shapes + i * 4 + k ) );
            sh->Left = 10 + ( ( i + 1 ) % 2 ? w / 2 : 0 ) + k * w + k * 10;
            sh->Top = 10 + i * w;
 
            sh->rp = rand( ) % 100 > 50 ? 5 : 6;
            sh->gp = rand( ) % 100 > 50 ? 3 : 4;
        }
    }
    for ( int i( 0 ), i_limit( 12 ); i < i_limit; ++i )
    {
        for ( int k( 0 ), k_limit( 12 ); k < k_limit; ++k )
        {
            if ( ( k != i ) &&
                InRange( shapes[ k ]->Left, shapes[ i ]->Left - w - 20, shapes[ i ]->Left + w + 20 ) &&
                InRange( shapes[ k ]->Top, shapes[ i ]->Top - w - 20, shapes[ i ]->Top + w + 20 ) )
            {
                shapes[ i ]->neibs[ k ] = shapes[ k ];
            }
        }
    }
}
 
// ---------------------------------------------------------------------------
__fastcall TForm1::~TForm1( )
{
 
}
 
// ---------------------------------------------------------------------------
 
// ---------------------------------------------------------------------------
void __fastcall TForm1::cell_mouse_down( TObject * Sender, TMouseButton Button, TShiftState Shift, int X, int Y )
{
    static int charge_id( 1 );
    t_cell_shape * sh( static_cast < t_cell_shape * > ( Sender ) );
    sh->charged = !sh->charged;
    sh->charge_id = charge_id++ ;
    sh->Repaint( );
}
 
// ---------------------------------------------------------------------------
void __fastcall TForm1::btn_stepClick( TObject * Sender )
{
    static int step( 0 );
    static_cast < TButton * > ( Sender )->Caption = "Шаг " + IntToStr( ++step );
    bool passed( false ), uncharged_exists( false );
    t_cell_shape::t_neibs::iterator last_uncharged;
 
    for ( int i( 0 ), i_limit( 12 ); i < i_limit; ++i, passed = false, uncharged_exists = false )
    {
        for ( t_cell_shape::t_neibs::iterator it( shapes[ i ]->neibs.begin( ) ), end( shapes[ i ]->neibs.end( ) );
            shapes[ i ]->charged && !shapes[ i ]->just_received && ( it != end ); ++it )
        {
            if ( !it->second->charged )
            {
                uncharged_exists = true;
                last_uncharged = it;
                if ( it->second->rp == 5 )
                {
                    shapes[ i ]->move_charge( it->second );
                    passed = true;
                }
            }
        }
        if ( ( shapes[ i ]->gp == 3 ) && uncharged_exists && !passed )
        {
            shapes[ i ]->move_charge( last_uncharged->second );
        }
    }
    for ( int i( 0 ), i_limit( 12 ); i < i_limit; ++i )
    {
        shapes[ i ]->rp = rand( ) % 100 > 50 ? 5 : 6;
        shapes[ i ]->gp = rand( ) % 100 > 50 ? 3 : 4;
        shapes[ i ]->just_received = false;
        shapes[ i ]->Repaint( );
    }
}