#include "stdio.h"
#include "stdlib.h"

unsigned min(unsigned a, unsigned b) { return a<b ? a : b; }

/**
* Détermine le cout minimal du cout de production et du cout de stockage
* satisfaisant la demande des commandes passées sur les n prochains jours.
* n     : Nombre de jours
* Emax  : Stockage maximale possible                        max(E)
* d[]   : Demande selon le jour i                           Di
* q[]   : Quantité de production maximale selon le jour i   Qi
* t[]   : Cout de stockage unitaire selon le jour i         Ti
* c[]   : Cout de production unitaire selon le jour i       Ci
*/
int cout_min(int n, int d[], int q[], int t[], int c[]) {
    int Emax = 0, e, k, x, ee=0;
    for (k=1;k<=n;k++) Emax+=q[k];
    Emax=35;
    unsigned int f[n+2][Emax+1];

    for(e=0;e<=Emax;e++) f[n+1][e]=1000;  
    f[n+1][0]=0;

    for(e=0;e<=Emax;e++) f[n][e]= e>=d[n]? t[n]*(e-d[n]) :d[n]-e<=q[n]? c[n]*(d[n]-e) : 1000;
    for(k=n-0;printf("\n"),k>0;k--)
    for(e=0;e<=Emax;e++)
    { unsigned z=1000;
      for(x=0;x<=q[k];x++) if(ee=e+x-d[k],ee>=0 && ee<=Emax) z=min(z,c[k]*x+t[k]*ee+f[k+1][ee]);
      f[k][e] = z;
      printf("%5d",z);
    }
    for (k = 1, e = 0; k <= n; k++) {
        unsigned int z = f[k][e];
        for(x=0;x<=q[k];x++) if(ee=e+x-d[k],ee>=0 && ee<=Emax) if(z==c[k]*x+t[k]*ee+f[k+1][ee]) break;
        printf("Au jour %d on produit %2d le stock devient %2d=%2d+%2d-%2d. Le cout est %2d=%2d+%2d*%2d+%2d*%2d\n",
                         k,             x,                 ee,  e,  x, d[k],     f[k][e],f[k+1][ee],c[k],x,t[k],ee);
        e=ee;
    }
    return f[1][0];
}


int main(int argc, char const *argv[]) {
    int n = 6;
    int d[] = {0,  5,  6,  6,  4,  5,  7};
    int q[] = {0, 10, 10, 10, 10, 10, 10};
    int c[] = {0,  3,  1,  2,  3,  1,  5};
    int t[] = {0,  2,  2,  1,  2,  2, 10};
    int cout = cout_min(n, d, q, t, c);
    printf("Cout min %d\n", cout);
    return 0;
}