Cod sursa(job #420311)
#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-1; i>=0; i--) {
for(j=0; j<=L; j++) {
for(p=0; p<=timp/A[i]; p++) {
if(j+p>=mla && M[i][j]+P[i+1][p]>=mlb) {
La[i+1]=p; Lb[i+1]=P[i+1][p];
mla=mla-p; mlb=mlb-P[i+1][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;
}