#ifndef DA_H_INCLUDED
#define DA_H_INCLUDED

#include <iostream>
#include <vector>
#include <string>

#include <stdio.h>
#include <cmath>

#include <algorithm>

#include <stdlib.h>
#include <fstream>
#include "../instance/indices.h"
#include "../instance/cinstance.h"
#include "../instance/cinstancereduite.h"
#include "cgraphe.h"
#include "cplan.h"
#include "clique.h"
#include <cplexcall.h>

using namespace std;
class Cplex;

class DualAscent {
    public :

    // Matrice A de taille 260 + 51 X nombre de colonnes
    vector< vector<double> * > A ;
    vector<vector<double> > recharge; //recharge[col][k] = quantité refuelée à la k-ème campagne, sur la col-ième colonne
    double a(int i, int j) const { return (*A[i])[j];}

    // vecteur de cout des colonnes
    vector<double> cost ;

    // fuel remaining
    vector<double> reste; // reste[col] = fuel restant en fin de scénar

    // vecteur de demande de production + affectation des colonnes
    vector<double> b ;
    
    // vecteur de production MAX remplie par toutes les centrales nucléaires si elles fonctionnaient toutes à PMAX + affectation colonne
    vector<double> b_PMAX ;

    // mêmes informations pour la colonne themique
    vector<double> *    therm_A ;
    vector<double>      therm_recharge ;
    double              therm_cost ;
    double              therm_reste ;

    // fonction
    DualAscent(vector<double>  dem, powerplants2_t nb_centrales)
    {
        b = dem ;
        b_PMAX = dem;
        // Ajout Affectation sur b et b_PMAX
        b_PMAX.resize( dem.size() + nb_centrales , 1 );
        b.resize( dem.size() + nb_centrales , 1 );

        double coutTherm = 0;
        double b_PMAX_coutTherm =0;

        for (week_t w=(week_t) 0;w<H();w=(week_t) (w+1) )
        {
        	// Remplissage de b_PMAX
        	double sumPMAXweek = 0;
        	for (powerplants2_t i = (powerplants2_t) 0 ; i < nb_centrales ; i = (powerplants2_t) (i+1) )
        		sumPMAXweek += CInstanceReduite::PMAX(i,w);
        	// std::cout << "sumPMAXweek = " << sumPMAXweek << std::endl;
        	b_PMAX[w] =sumPMAXweek;
        	
        	// remplissage de coutTherm et de b_PMAX_coutTherm
            coutTherm += b[w]*CInstanceReduite::ThCout(w);
            // std::cout<<"w = "<<w << " & b_PMAX.size() = "<<b_PMAX.size() << std::endl; 
			b_PMAX_coutTherm += b_PMAX[w]*CInstanceReduite::ThCout(w);
        }
        
        
        std::vector<CPlan::CArc> vide;
#ifdef PMAX_POUR_DEM
        ajoutcolonne(&b_PMAX, b_PMAX_coutTherm, vide, 0);
#else
        ajoutcolonne(&b, coutTherm, vide, 0);
#endif
    }

    // destructeur
    ~DualAscent() ;

    protected:
        void ajoutcolonne( vector<double> * colonne, double cou, std::vector<CPlan::CArc>& arcs, double resteCol );

   public:
        vector<CPlan *> m_plans ;
        void ajoutcolonne(CPlan* cplan);

        /// Fonction à appeler pour le calcul des Ui Vi
        void boucle(vector<double>& u, vector <double> & lambda, double& BorneSup, const char*x=NULL)  const;

    private:

        vector<double> calculePoids() const;

        vector<double> DAcalc (const vector<double> & ,vector<double> &, const vector<double> &) const;

        vector<double> sousgrad(vector<double> &, const vector<double> &,const vector<double> &,
                                const vector<double> &,const vector<double> &) const;

        // Recherche de borne sup upper bound
        bool testUB () const ;

        double recherche_UB () const ;
		  double recherche_UB2 (const vector<double> &q) const ;

        /**
         * Supprime la colonne correspondant à la centrale thermique.
         * À n'appeler qu'une unique fois, juste après la génération de colonne,
         * et juste avant la préparation des données pour CPlex
         *
         * Appelé par la classe CPlex
         */
    public:
        void vireColonneThermique() ;
        void remetColonneThermique();


        //temp : méthode de debug
        void testL(const vector<double>& L, const string& mes) const
        {
        	for ( size_t i=0 ; i < L.size() ; i++ )
        	{
        		if ( (L[i]==INFINI) || L[i]!=L[i] )
        		{
        			cerr<<mes<<endl;
        			cerr<<"L["<<i<<"] = "<<L[i]<<endl;
        			exit(-1);
        		}
        	}
        }

};

#endif // DA_H_INCLUDED