http://acm.hdu.edu.cn/showproblem.php?pid=6804
题解的方法是取T=sqrt(n+m) * 1000 * 2的容量,然后随机打乱两个集合的物品顺序,然后去更新这个有限制大小的背包,直接在[-T,T]直接更新,期望总能在这一部分碰撞找到最优解。
然而这题过了一车人,154个估计得有140个1e10暴力跑过去的,我加了一大把优化还跑步过去,反而最简单的能跑过,反向优化
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
mt19937 rnd(time(NULL));
const int maxl=2010;
const int maxm=maxl*100,mid=maxm/2;
const ll inf=1ll<<60;
int n,m;ll ans;
struct node
{
int t,w,v;
}a[maxl];
ll f[2][maxm];
inline void prework()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n+m;i++)
{
scanf("%d%d",&a[i].w,&a[i].v);
if(i<=n) a[i].t=1;else a[i].t=2;
}
shuffle(a+1,a+1+n+m,rnd);
}
inline umx(ll &x,ll y)
{
if(x<y) x=y;
}
inline void mainwork()
{
int l=mid,r=mid,d,lim=sqrt(n+m)*1000*2,nl,nr;
f[0][mid]=0;
for(int i=1;i<=n+m;i++)
{
d=i&1;
if(a[i].t==1)
{
nr=min(r+a[i].w,mid+lim);
for(int j=l;j<=r;j++)
f[d][j]=f[d^1][j];
for(int j=r+1;j<=nr;j++)
f[d][j]=-inf;
for(int j=l;j<=nr-a[i].w;j++)
umx(f[d][j+a[i].w],f[d^1][j]+a[i].v);
r=nr;
}
else
{
nl=max(l-a[i].w,mid-lim);
for(int j=l;j<=r;j++)
f[d][j]=f[d^1][j];
for(int j=nl;j<=l-1;j++)
f[d][j]=-inf;
for(int j=nl+a[i].w;j<=r;j++)
umx(f[d][j-a[i].w],f[d^1][j]+a[i].v);
l=nl;
}
}
ans=f[(n+m)&1][mid];
}
inline void print()
{
printf("%lld\n",ans);
}
int main()
{
int t;
scanf("%d",&t);
for(int i=1;i<=t;i++)
{
prework();
mainwork();
print();
}
return 0;
}