【XJOI娱乐赛】C.病毒研究

题目

http://115.236.49.52:83/contest/1351

思路

比赛的时候状态设错了……

设f[l][r]为现在可以确定活性在l到r之间，需要的期望代价乘上(r-l+1)=dist的值。

转移：f[l][r]=min(f[l-w[i]][r-w[i]]+v[i]*dist)

复杂度O(a[n]²m)

我们先对所有操作跑一次完全背包，这样我们就可以使用一次操作来替代连续的多次操作。进行了完全背包以后，我们获得了新的vi和wi 。

这样我们可以强制使转移跨状态

使用前缀和优化，复杂度O（a[n]m）

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=2077;
const ll inf=1e14;
int n,m,a[N],S[N],nxt[N];
ll sum[N],f[N][2],g[N];
ll aii(int l,int r)
{
	if(r<=a[1]) return 0;
	int pl=nxt[l],pr=S[r];
	if(pl>pr) return inf;
	ll yjy=sum[pr]-sum[pl];
	yjy+=f[l][1]; 
	yjy+=f[r][0];
	return yjy;
}
ll work(int l,int r)
{
	if(r<=a[1]) return 0;
	ll yjy=inf;
	for(int i=1; i<=l-1; i++)
		yjy=min(yjy,aii(l-i,r-i)+g[i]*(r-l+1));
	return yjy;
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		memset(S,0,sizeof S);
		memset(nxt,0,sizeof nxt);
		memset(f,0,sizeof f);
		memset(g,0,sizeof g);
		memset(sum,0,sizeof sum);
		scanf("%d%d",&n,&m);
		for(int i=1; i<=n; i++) scanf("%d",&a[i]);
		for(int i=1; i<=n; i++)
			for(int j=a[i-1]+1; j<=a[i]; j++)
			{
				S[j]=i-1; nxt[j]=i;
			}
		for(int i=1; i<=a[n]; i++) g[i]=inf;
		for(int i=1; i<=m; i++)
		{
			int v,c; 
			scanf("%d%d",&v,&c);
			for(int j=c; j<=a[n]; j++)
				g[j]=min(g[j],g[j-c]+v);
		}
		for(int i=1; i<=n; i++)
		{
			for(int j=a[i-1]+1; j<=a[i]; j++) f[j][0]=work(a[i-1]+1,j);
			for(int j=a[i-1]+1; j<=a[i]; j++) f[j][1]=work(j,a[i]);
			sum[i]=sum[i-1]+f[a[i]][0];
		}
		if(sum[n]>=inf) printf("-1\n");
		else printf("%lld\n",sum[n]);
	}
}