MINOKAWA with program lpsolve :
/* Objective function /
min: Varkonor + Vbistot + Vcrokite + Vgneiss + Vspod + Vochre;

/ Variable bounds */
Varkonor = 16Qarkonor;
Vbistot  = 16Qbistot;
Vcrokite = 16Qcrokite;
Vgneiss  =  5Qgneiss;
Vspod    = 16Qspod;
Vochre   =  8Qochre;

Tfromarkonor = 22000Qarkonor;
Mfromarkonor = 2500Qarkonor;
Cfromarkonor = 320Qarkonor;

Pfrombistot  = 12000Qbistot;
Cfrombistot  = 100Qbistot;
Zfrombistot  = 450Qbistot;

Tfromcrokite = 21000Qcrokite;
Nfromcrokite = 760Qcrokite;
Zfromcrokite = 135Qcrokite;

Pfromgneiss  = 2200Qgneiss;
Mfromgneiss  = 2400Qgneiss;
Ifromgneiss  = 300Qgneiss;

Tfromspod    = 56000Qspod;
Pfromspod    = 12050Qspod;
Mfromspod    = 2100Qspod;
Ifromspod    = 450Qspod;

Tfromochre   = 10000Qochre;
Ifromochre   = 1600Qochre;
Nfromochre   = 120Qochre;

T: 0.85Tfromarkonor+0.85Tfromcrokite+0.85Tfromspod+0.85Tfromochre >= 71267300;
P: 0.85Pfrombistot+0.85Pfromgneiss+0.85Pfromspod >= 17526635;
M: 0.85Mfromarkonor+0.85Mfromgneiss+0.85Mfromspod >= 6486403;
I: 0.85Ifromgneiss+0.85Ifromspod+0.85Ifromochre >= 1057880;
N: 0.85Nfromcrokite+0.85Nfromochre >= 303763;
Z: 0.85Zfrombistot+0.85Zfromcrokite >= 107829;
C: 0.85Cfromarkonor+0.85Cfrombistot >= 49239;