ForwardBisimulation.h

Go to the documentation of this file.

 
/*
 * This file is part of Algorithms library toolkit.
 * Copyright (C) 2017 Jan Travnicek (jan.travnicek@fit.cvut.cz)
 
 * Algorithms library toolkit is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 
 * Algorithms library toolkit is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 
 * You should have received a copy of the GNU General Public License
 * along with Algorithms library toolkit.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#pragma once
 
#include <alib/set>
#include <alib/map>
 
#include <automaton/TA/DFTA.h>
#include <automaton/TA/NFTA.h>
 
#include <automaton/TA/UnorderedDFTA.h>
#include <automaton/TA/UnorderedNFTA.h>
 
#include <automaton/FSM/DFA.h>
#include <automaton/FSM/NFA.h>
 
#include <automaton/TA/ArcFactoredDeterministicZAutomaton.h>
#include <automaton/TA/ArcFactoredNondeterministicZAutomaton.h>
 
namespace automaton {
 
namespace properties {
 
class ForwardBisimulation {
    template < class StateType >
    static ext::set < ext::pair < StateType, StateType > > initial ( const ext::set < StateType > & states, const ext::set < StateType > & finals ) {
        ext::set < ext::pair < StateType, StateType > > init;
 
        for ( const StateType & a : states ) {
            for ( const StateType & b : states ) {
                if ( finals.count ( a ) == finals.count ( b ) ) {
                    init.insert ( ext::make_pair ( a, b ) );
                    init.insert ( ext::make_pair ( b, a ) );
                }
            }
        }
 
        return init;
    }
 
public:
 
    template < class SymbolType, class StateType >
    static ext::set < ext::pair < StateType, StateType > > forwardBisimulation ( const automaton::DFA < SymbolType, StateType > & fta );
 
    template < class SymbolType, class StateType >
    static ext::set < ext::pair < StateType, StateType > > forwardBisimulation ( const automaton::NFA < SymbolType, StateType > & fta );
 
    template < class SymbolType, class StateType >
    static ext::set < ext::pair < StateType, StateType > > forwardBisimulation ( const automaton::ArcFactoredDeterministicZAutomaton < SymbolType, StateType > & afdza );
 
    template < class SymbolType, class StateType >
    static ext::set < ext::pair < StateType, StateType > > forwardBisimulation ( const automaton::ArcFactoredNondeterministicZAutomaton < SymbolType, StateType > & afnza );
 
    template < class SymbolType, class StateType >
    static ext::set < ext::pair < StateType, StateType > > forwardBisimulation ( const automaton::DFTA < SymbolType, StateType > & fta );
 
    template < class SymbolType, class StateType >
    static ext::set < ext::pair < StateType, StateType > > forwardBisimulation ( const automaton::NFTA < SymbolType, StateType > & fta );
 
    template < class SymbolType, class StateType >
    static ext::set < ext::pair < StateType, StateType > > forwardBisimulation ( const automaton::UnorderedDFTA < SymbolType, StateType > & fta );
 
    template < class SymbolType, class StateType >
    static ext::set < ext::pair < StateType, StateType > > forwardBisimulation ( const automaton::UnorderedNFTA < SymbolType, StateType > & fta );
};
 
template < class SymbolType, class StateType >
ext::set < ext::pair < StateType, StateType > > ForwardBisimulation::forwardBisimulation ( const automaton::DFA < SymbolType, StateType > & fta ) {
    return forwardBisimulation ( automaton::NFA < SymbolType, StateType > ( fta ) );
}
 
template < class SymbolType, class StateType >
ext::set < ext::pair < StateType, StateType > > ForwardBisimulation::forwardBisimulation ( const automaton::NFA < SymbolType, StateType > & fta ) {
    ext::set < ext::pair < StateType, StateType > > forwardBisimulation = initial ( fta.getStates ( ), fta.getFinalStates ( ) );
 
    bool changed;
    do {
        changed = false;
 
        for ( const StateType & p : fta.getStates ( ) ) {
            for ( const StateType & q : fta.getStates ( ) ) {
                if ( ! forwardBisimulation.contains ( ext::make_pair ( p, q ) ) )
                    continue;
 
                for ( const std::pair < const ext::pair < StateType, SymbolType >, StateType > & pTransition : fta.getTransitions ( ) ) {
                    if ( pTransition.first.first == p ) {
                        bool exists = false;
                        for ( const std::pair < const ext::pair < StateType, SymbolType >, StateType > & qTransition : fta.getTransitions ( ).equal_range ( std::make_pair ( q, pTransition.first.second ) ) ) {
                            if ( forwardBisimulation.contains ( ext::make_pair ( pTransition.second, qTransition.second ) ) ) {
                                exists = true;
                                break;
                            }
                        }
 
                        if ( ! exists ) {
                            forwardBisimulation.erase ( ext::make_pair ( p, q ) );
                            forwardBisimulation.erase ( ext::make_pair ( q, p ) );
                            changed = true;
                        }
                    }
                }
            }
        }
    } while ( changed );
 
    return forwardBisimulation;
}
 
template < class SymbolType, class StateType >
ext::set < ext::pair < StateType, StateType > > ForwardBisimulation::forwardBisimulation ( const automaton::ArcFactoredDeterministicZAutomaton < SymbolType, StateType > & afdza ) {
    return forwardBisimulation ( automaton::ArcFactoredNondeterministicZAutomaton < SymbolType, StateType > ( afdza ) );
}
 
template < class SymbolType, class StateType >
ext::set < ext::pair < StateType, StateType > > ForwardBisimulation::forwardBisimulation ( const automaton::ArcFactoredNondeterministicZAutomaton < SymbolType, StateType > & afnza ) {
    ext::set < ext::pair < StateType, StateType > > forwardBisimulation = initial ( afnza.getStates ( ), afnza.getFinalStates ( ) );
 
    bool changed;
    do {
        changed = false;
 
        for ( const StateType & p : afnza.getStates ( ) ) {
            for ( const StateType & q : afnza.getStates ( ) ) {
                if ( ! forwardBisimulation.contains ( ext::make_pair ( p, q ) ) )
                    continue;
 
                for ( const std::pair < const ext::variant < SymbolType, ext::pair < StateType, StateType > >, StateType > & pTransition : afnza.getTransitions ( ) ) {
                    if ( pTransition.first.template is < ext::pair < StateType, StateType > > ( ) && pTransition.first.template get < ext::pair < StateType, StateType > > ( ).first == p ) {
                        bool exists = false;
                        for ( const std::pair < const ext::variant < SymbolType, ext::pair < StateType, StateType > >, StateType > & qTransition : afnza.getTransitions ( ) ) {
                            if ( qTransition.first.template is < ext::pair < StateType, StateType > > ( ) && qTransition.first.template get < ext::pair < StateType, StateType > > ( ).first == q && pTransition.first.template get < ext::pair < StateType, StateType > > ( ).second == qTransition.first.template get < ext::pair < StateType, StateType > > ( ).second ) {
                                if ( forwardBisimulation.contains ( ext::make_pair ( pTransition.second, qTransition.second ) ) ) {
                                    exists = true;
                                    break;
                                }
                            }
                        }
 
                        if ( ! exists ) {
                            forwardBisimulation.erase ( ext::make_pair ( p, q ) );
                            forwardBisimulation.erase ( ext::make_pair ( q, p ) );
                            changed = true;
                        }
                    }
                    if ( pTransition.first.template is < ext::pair < StateType, StateType > > ( ) && pTransition.first.template get < ext::pair < StateType, StateType > > ( ).second == p ) {
                        bool exists = false;
                        for ( const std::pair < const ext::variant < SymbolType, ext::pair < StateType, StateType > >, StateType > & qTransition : afnza.getTransitions ( ) ) {
                            if ( qTransition.first.template is < ext::pair < StateType, StateType > > ( ) && qTransition.first.template get < ext::pair < StateType, StateType > > ( ).second == q && pTransition.first.template get < ext::pair < StateType, StateType > > ( ).first == qTransition.first.template get < ext::pair < StateType, StateType > > ( ).first ) {
                                if ( forwardBisimulation.contains ( ext::make_pair ( pTransition.second, qTransition.second ) ) ) {
                                    exists = true;
                                    break;
                                }
                            }
                        }
 
                        if ( ! exists ) {
                            forwardBisimulation.erase ( ext::make_pair ( p, q ) );
                            forwardBisimulation.erase ( ext::make_pair ( q, p ) );
                            changed = true;
                        }
                    }
                }
            }
        }
    } while ( changed );
 
    return forwardBisimulation;
}
 
template < class SymbolType, class StateType >
ext::set < ext::pair < StateType, StateType > > ForwardBisimulation::forwardBisimulation ( const automaton::DFTA < SymbolType, StateType > & fta ) {
    return forwardBisimulation ( automaton::NFTA < SymbolType, StateType > ( fta ) );
}
 
template < class SymbolType, class StateType >
ext::set < ext::pair < StateType, StateType > > ForwardBisimulation::forwardBisimulation ( const automaton::NFTA < SymbolType, StateType > & fta ) {
    ext::set < ext::pair < StateType, StateType > > forwardBisimulation = initial ( fta.getStates ( ), fta.getFinalStates ( ) );
 
    bool changed;
    do {
        changed = false;
 
        for ( const StateType & p : fta.getStates ( ) ) {
            for ( const StateType & q : fta.getStates ( ) ) {
                if ( ! forwardBisimulation.contains ( ext::make_pair ( p, q ) ) )
                    continue;
 
                for ( const std::pair < const ext::pair < common::ranked_symbol < SymbolType >, ext::vector < StateType > >, StateType > & pTransition : fta.getTransitions ( ) ) {
                    for ( size_t i = 0; i < pTransition.first.second.size ( ); ++ i ) {
                        if ( pTransition.first.second [ i ] == p ) {
                            ext::vector < StateType > copy = pTransition.first.second;
                            copy [ i ] = q;
 
                            bool exists = false;
                            for ( const std::pair < const ext::pair < common::ranked_symbol < SymbolType >, ext::vector < StateType > >, StateType > & qTransition : fta.getTransitions ( ).equal_range ( std::make_pair ( pTransition.first.first, std::move ( copy ) ) ) ) {
                                if ( forwardBisimulation.contains ( ext::make_pair ( pTransition.second, qTransition.second ) ) ) {
                                    exists = true;
                                    break;
                                }
                            }
 
                            if ( ! exists ) {
                                forwardBisimulation.erase ( ext::make_pair ( p, q ) );
                                forwardBisimulation.erase ( ext::make_pair ( q, p ) );
                                changed = true;
                            }
                        }
                    }
                }
            }
        }
    } while ( changed );
 
    return forwardBisimulation;
}
 
template < class SymbolType, class StateType >
ext::set < ext::pair < StateType, StateType > > ForwardBisimulation::forwardBisimulation ( const automaton::UnorderedDFTA < SymbolType, StateType > & fta ) {
    return forwardBisimulation ( automaton::UnorderedNFTA < SymbolType, StateType > ( fta ) );
}
 
template < class SymbolType, class StateType >
ext::set < ext::pair < StateType, StateType > > ForwardBisimulation::forwardBisimulation ( const automaton::UnorderedNFTA < SymbolType, StateType > & fta ) {
    ext::set < ext::pair < StateType, StateType > > forwardBisimulation = initial ( fta.getStates ( ), fta.getFinalStates ( ) );
 
    bool changed;
    do {
        changed = false;
 
        for ( const StateType & p : fta.getStates ( ) ) {
            for ( const StateType & q : fta.getStates ( ) ) {
                if ( ! forwardBisimulation.contains ( ext::make_pair ( p, q ) ) )
                    continue;
 
                for ( const std::pair < const ext::pair < common::ranked_symbol < SymbolType >, ext::multiset < StateType > >, StateType > & pTransition : fta.getTransitions ( ) ) {
                    for ( const StateType & state : pTransition.first.second ) {
                        if ( state == p ) {
                            ext::multiset < StateType > copy = pTransition.first.second;
                            copy.erase ( copy.find ( state ) );
                            copy.insert ( q );
 
                            bool exists = false;
                            for ( const std::pair < const ext::pair < common::ranked_symbol < SymbolType >, ext::multiset < StateType > >, StateType > & qTransition : fta.getTransitions ( ).equal_range ( std::make_pair ( pTransition.first.first, std::move ( copy ) ) ) ) {
                                if ( forwardBisimulation.contains ( ext::make_pair ( pTransition.second, qTransition.second ) ) ) {
                                    exists = true;
                                    break;
                                }
                            }
 
                            if ( ! exists ) {
                                forwardBisimulation.erase ( ext::make_pair ( p, q ) );
                                forwardBisimulation.erase ( ext::make_pair ( q, p ) );
                                changed = true;
                            }
                        }
                    }
                }
            }
        }
    } while ( changed );
 
    return forwardBisimulation;
}
 
} /* namespace properties */
 
} /* namespace automaton */