#ifndef DA_H_INCLUDED
#define DA_H_INCLUDED

#include <iostream>
#include <vector>
#include <string.h>

#include <stdio.h>
#include <math.h>

#include <algorithm>

#include <stdlib.h>
#include <fstream>
#include "indices.h"
#include "cinstance.h"
#include "cinstancereduite.h"
#include "cplan.h"
#include "clique.h"

using namespace std;

class DualAscent {
    public :

    // Matrice A de taille 260 + 51 X nombre de colonnes
    vector< vector<double> * > A ;
    double a(int i, int j) const { return (*A[i])[j];}

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

    // vecteur de demande de production + affectation des colonnes
    vector<double> b ;

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

        double coutTherm = 0;

        for (week_t i=(week_t) 0;i<H();i=(week_t) (i+1) )
        {
            coutTherm += b[i]*CInstanceReduite::ThCout(i);
        }
        ajoutcolonne(&b, coutTherm);
    }

    // destructeur
    ~DualAscent()
    {
//        for ( size_t i=0; i<A.size(); i++) delete A[i];
    }

    protected:
    void ajoutcolonne( vector<double> * colonne, double cou )
    {
        A.push_back(colonne);
        cost.push_back(cou);
    }

   public:
    void ajoutcolonne(CPlan* cplan){
        Clique::ajoutcolonne(cplan, A.size());
        ajoutcolonne(&cplan->m_colonne, cplan->m_cout);
    }

    void print();

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

    private:

        vector<double> calculePoids() const;

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

        vector <int> EcritureCplex(const vector <double> &, const char *) const;

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

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

        // operateur binaire
        bool epsiqual( const double &, const double &) const;

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

#endif // DA_H_INCLUDED