Cod sursa(job #420312)

Utilizator vladiiIonescu Vlad vladii Data 18 martie 2010 19:33:05
Problema Lapte Scor 30
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.77 kb
#include <iostream>
using namespace std;

int N, L;
int A[101], B[101], La[101], Lb[101];
int P[101][101], M[101][101];

int Verifica(int timp) {
    int i, j, p;
    for(i=1; i<=N; i++) {
         for(j=0; j<=timp/A[i]; j++) {
              P[i][j]=(timp-j*A[i])/B[i];
         }
    }
    for(i=0; i<=timp/A[1]; i++) {
         if(i>=L) { M[1][L]=max(M[1][L], P[1][i]); 
              if(M[1][L]>=L) { return 1; }
              continue; 
         }
         M[1][i]=P[1][i];
    }
    for(i=2; i<=N; i++) {
         for(j=0; j<=L; j++) {
              for(p=0; p<=timp/A[i]; p++) {
                   //M[i-1][j], P[i][p]
                   if(p+j>=L) { M[i][L]=max(M[i][L], M[i-1][j]+P[i][p]); 
                        if(M[i][L]>=L) { return 1; }
                        continue; 
                   }
                   M[i][j+p]=max(M[i][j+p], M[i-1][j]+P[i][p]);
              }
         }
    }
    return 0;
}

int main() {
    FILE *f1=fopen("lapte.in", "r"), *f2=fopen("lapte.out", "w");
    int i, j, p, sol=0, mla, mlb, ok=0;
    fscanf(f1, "%d%d", &N, &L);
    for(i=1; i<=N; i++) {
         fscanf(f1, "%d%d", &A[i], &B[i]);
    }
    int ls=1, ld=101;
    while(ls<=ld) {
         int mij=(ls+ld)/2;
         if(Verifica(mij)) {
              //e posibil, deci caut un timp mai mic
              sol=mij;
              ld=mij-1;
         }
         else {
              //timp mai mare
              ls=mij+1;
         }
         for(i=1; i<=N; i++) {
              for(j=0; j<=L; j++) {
                   P[i][j]=0; M[i][j]=0;
              }
         }
    }
    fprintf(f2, "%d\n", sol);
    //reconstituiesc solutia
    int timp=sol;
    for(i=1; i<=N; i++) {
         for(j=0; j<=timp/A[i]; j++) {
              P[i][j]=(timp-j*A[i])/B[i];
         }
    }
    for(i=0; i<=timp/A[1]; i++) {
         if(i>=L) { M[1][L]=max(M[1][L], P[1][i]); continue; }
         M[1][i]=P[1][i];
    }
    for(i=2; i<=N; i++) {
         for(j=0; j<=L; j++) {
              for(p=0; p<=timp/A[i]; p++) {
                   if(p+j>=L) { M[i][L]=max(M[i][L], M[i-1][j]+P[i][p]); continue; }
                   M[i][j+p]=max(M[i][j+p], M[i-1][j]+P[i][p]);
              }
         }
    }
    mla=L; mlb=L; A[0]=1;
    for(i=N; i>=1; i--) {
         for(j=0; j<=L; j++) {
              for(p=0; p<=timp/A[i]; p++) {
                   if(j+p>=mla && M[i-1][j]+P[i][p]>=mlb) {
                        La[i]=p; Lb[i]=P[i][p];
                        mla=mla-p; mlb=mlb-P[i][p];
                        ok=1; break;
                   }
              }
              if(ok) { ok=0; break; }
         }
    }
    for(i=1; i<=N; i++) {
         fprintf(f2, "%d %d\n", La[i], Lb[i]);
    }
    fclose(f1); fclose(f2);
    return 0;
}