#include "selecteur_1.h"

#include <cinstance.h>
#include <cinstancereduite.h>
#include <cplan.h>

Selecteur_1::Selecteur_1()
 : Selecteur()
{
}


Selecteur_1::~Selecteur_1()
{
}


void Selecteur_1::setproblemdata()
{
    /* Création des vecteurs de contraintes de clique */
    m_cliques.clear();
    ajoute_contraintes ( C14, C21 );

    /* définition des zones */
    m_cplex.addColonne( "colonnes nucléaires",      CPlan::getNbPlans() );
    m_cplex.addColonne( "production thermique",     H()                             );

    m_cplex.addLigne( "production=demande",         H()                             );
    m_cplex.addLigne( "unicité colonne nucléaire",  I()                             );
    m_cplex.addLigne( "contraintes de clique",      m_cliques.size()                );

    /// l'objectif est de minimiser une combinaison linéaire des éléments de zobj (en ligne)
    ///     donc de trouver des valeurs pour NUMCOLS variables
    // coeff: coût des colonnes nucléaires (centrale nucléaire)
    // variables associées : binaires (0 ou 1)
    m_cplex.setColonne( "colonnes nucléaires" );
    for ( size_t p=0; p<CPlan::getNbPlans(); p++)
        m_cplex.addToObj( p,
                          0,
                          1,
                          CPlexModel::Binary,
                          CPlan::getPlan(p)->m_cout.v );

    // coeff: coût de la production thermique unitaire par semaine
    // variables associées : énergie produite, continue (entre pmin et pmax)
    m_cplex.setColonne( "production thermique");
    for ( week_t h=(week_t) 0; h<H(); ++h)
        m_cplex.addToObj( h,
                          CInstanceReduite::ThProdMin(h).v,
                                  CInstanceReduite::ThProdMax(h).v,
                                          CPlexModel::Continuous,
                                          CInstanceReduite::ThCout(h).v );


    /// remplissage de Z (contraintes "horizontales")

    /// satisfaction de la demande
    m_cplex.setLigne( "production=demande" );

    // partie droite : demande, opérateur =
    for (week_t h=(week_t) 0; h<H(); ++h)
        m_cplex.addToRHS( h, CPlexModel::Equals, CInstanceReduite::getEDEM_min()[h].v );

    // production des centrales nucléaires par semaine
    m_cplex.setColonne( "colonnes nucléaires" );
    for (size_t p=0; p<CPlan::getNbPlans(); p++)
        for (week_t h=(week_t) 0; h<H(); ++h)
            m_cplex.addToZ( h, p, CPlan::getPlan(p)->m_colonne[h].v );

    // production de la centrale thermique par semaine (matrice identité)
    m_cplex.setColonne( "production thermique" );
    for (week_t h=(week_t)0; h<H(); ++h)
        m_cplex.addToZ( h, h, 1 );

    /// sélection de la colonne nucléaire
    m_cplex.setLigne( "unicité colonne nucléaire" );

    // partie droite : 1, opérateur =
    for (powerplants2_t i=(powerplants2_t) 0; i<I(); ++i)
        m_cplex.addToRHS( i, CPlexModel::Equals, 1 );
    // partie gauche : 1 sur la centrale nucléaire correspondant à chaque colonne
    m_cplex.setColonne( "colonnes nucléaires" );
    for (size_t p=0; p<CPlan::getNbPlans(); p++)
        m_cplex.addToZ( CPlan::getPlan(p)->m_i, p, 1 );

    /// contraintes de cliques
    m_cplex.setLigne( "contraintes de clique" );
    m_cplex.setColonne( "colonnes nucléaires" );

    // partie droite : partie droite de la contrainte de clique
    // partie gauche : partie gauche de la contrainte de clique
    size_t c=0 ;
    for (std::map<std::vector<double>, double>::iterator it=m_cliques.begin();
         it!=m_cliques.end(); c++, it++)
    {
        m_cplex.addToRHS( c, CPlexModel::LowerThan, it->second );
        for (size_t p=0; p<CPlan::getNbPlans(); p++)
            m_cplex.addToZ( c, p, it->first[p] );
    }

}  /* END setproblemdata */