又是一道揹包......
這裏寫的比較好:http://blog.csdn.net/huyuncong/article/details/6884369
題意不復雜:給定n個揹包,m個物品,物品大小不能用dp過,但是有一個很詭異的性質:對於任意兩個物品,一定滿足一個物品的大小是另一個物品大小的整數倍,求最多可以在揹包中放入幾個物品。
這裏的關鍵條件:倍數關係,是需要使用進制來考慮的,把物品出現過大小當成進制,最後把揹包大小用這種進制來表示,然後從低位往高位貪心即可(說的是相當不清晰啊),這裏hyc寫得還可以,就進上面那個地址看算了。
program lmd;
var
fod,num,sys,bag:array[0..1000000] of int64;
have:array[0..1000000]of int64;
tot,i,j,bj,ans,n,m:longint;
procedure inf;
begin
assign(input,'rumor.in');
assign(output,'rumor.out');
reset(input); rewrite(output);
end;
procedure ouf;
begin
close(input); close(output);
end;
procedure init;
begin
read(n,m);
for i:=1 to n do read(bag[i]);
for i:=1 to m do read(fod[i]);
end;
procedure sort(l,r:longint);
var i,j,x,tmp:longint;
begin
i:=l; j:=r; x:=fod[(l+r) shr 1];
repeat
while fod[i]<x do inc(i);
while fod[j]>x do dec(j);
if i<= j then
begin
tmp:=fod[i]; fod[i]:=fod[j]; fod[j]:=tmp;
inc(i); dec(j);
end;
until i>j;
if i<r then sort(i,r);
if l<j then sort(l,j);
end;
procedure prepare;
begin
sort(1,m);
tot:=0;
for i:=1 to m do
if fod[i]<>fod[i-1] then
begin
inc(tot); sys[tot]:=fod[i];
end;
fillchar(num,sizeof(num),0);
sys[tot+1]:=maxlongint;
for i:=1 to n do
begin
for j:=1 to tot do
if sys[j]>bag[i] then break
else
inc(num[j],(bag[i] mod sys[j+1]) div sys[j]);
end;
for i:=1 to tot do
have[i]:=sys[i]*num[i];
end;
procedure main;
begin
bj:=1; ans:=0;
for i:=1 to m do
begin
while ((have[bj]<fod[i]) or (sys[bj]<fod[i])) and (bj<tot) do
inc(bj);
if (have[bj]<fod[i]) or (sys[bj]<fod[i]) then break;
dec(have[bj],fod[i]); inc(ans);
end;
write(ans);
end;
begin
inf;
init;
prepare;
main;
ouf;
end.