洛谷 普及/提高- DP复习

P2370 yyy2015c01 的 U 盘

        排序的背包问题。

	bool cmp(thing a,thing b){
	    return a.w<b.w; //按体积由小到大
    }
    
    sort(a+1,a+1+n,cmp);
	for(int i=1;i<=n;i++){
		for(int j=s;j>=a[i].w;j--){
			if(dp[j-a[i].w]+a[i].v>dp[j]){
				dp[j]=dp[j-a[i].w]+a[i].v;	
			}	
			if(dp[s]>=p){
				printf("%d",a[i].w);
				return 0; 
			}
		}
	}

---

P3133 [USACO16JAN]Radio Contact G

        模拟性DP，难度不大，就是比较复杂。

int cal(int a,int b){
	int dx=step1[a][0]-step2[b][0];
	int dy=step1[a][1]-step2[b][1];
	return dx*dx+dy*dy;
}//计算距离

	step1[0][0]=fx,step1[0][1]=fy;
	step2[0][0]=bx,step2[0][1]=by;
	for(int i=1;i<=l1;i++){
		if(fj[i]=='W')
			step1[i][0]=step1[i-1][0]-1,step1[i][1]=step1[i-1][1];
		if(fj[i]=='S')
			step1[i][0]=step1[i-1][0],step1[i][1]=step1[i-1][1]-1;
		if(fj[i]=='E')
			step1[i][0]=step1[i-1][0]+1,step1[i][1]=step1[i-1][1];
		if(fj[i]=='N')
			step1[i][0]=step1[i-1][0],step1[i][1]=step1[i-1][1]+1;
	}
	for(int i=1;i<=l2;i++){
		if(be[i]=='W')
			step2[i][0]=step2[i-1][0]-1,step2[i][1]=step2[i-1][1];
		if(be[i]=='S')
			step2[i][0]=step2[i-1][0],step2[i][1]=step2[i-1][1]-1;
		if(be[i]=='E')
			step2[i][0]=step2[i-1][0]+1,step2[i][1]=step2[i-1][1];
		if(be[i]=='N')
			step2[i][0]=step2[i-1][0],step2[i][1]=step2[i-1][1]+1;
	}
	memset(dp,0x3f,sizeof(dp));
	dp[0][0]=0;
	for(int i=1;i<=l1;i++){
		dp[i][0]=dp[i-1][0]+cal(i,0);
	}
	for(int i=1;i<=l2;i++){
		dp[0][i]=dp[0][i-1]+cal(0,i);
	}
	for(int i=1;i<=l1;i++){
		for(int j=1;j<=l2;j++){
			int e=cal(i,j);
			dp[i][j]=min(dp[i][j],min(dp[i-1][j-1],min(dp[i][j-1],dp[i-1][j]))+e);
		}
	}

---

P2049 魔术棋子

        有技巧的DP。需要我们对于取模的操作进行转移。

/*dp[i][j][l]表示在位置(i,j)能不能得到l。*/ 
int main()
{
    read(m),read(n),read(k);
    for(int i=1;i<=m;i++)
    {
        for(int j=1;j<=n;j++)
        {
            read(num[i][j]);
            num[i][j]%=k;    //一开始就mod 
        }
    }
    dp[1][1][num[1][1]]=true;    //初始化位置(1,1)的数 
    for(int i=1;i<=m;i++)
     for(int j=1;j<=n;j++)
      for(int l=0;l<k;l++)    //因为mod k 后得到的数一定小于k，所以从0到k枚举 
       if(!dp[i][j][l*num[i][j]%k])        //没有计算过 
        dp[i][j][l*num[i][j]%k]=dp[i-1][j][l]||dp[i][j-1][l];    
//l*num[i][j]%k表示当前格子数乘从左边或上边传下来的数l再mod k，dp[i-1][j][l]和dp[i][j-1][l]表示在上方或左方能不能得到l 
    for(int i=0;i<k;i++)    //找个数 
    {
        if(dp[m][n][i]) ans++;
    }
    printf("%d\n",ans);
    for(int i=0;i<k;i++)    //找可能得到的数 
    {
        if(dp[m][n][i]) printf("%d ",i);
    }
    return 0;
}

---

P1566 加等式

	for(int i=1;i<=l;i++)
		s2[i]=s1[l-i+1];
	for(int i=1;i<=l;i++){
		for(int j=1;j<=l;j++){
			if(s1[i]==s2[j]){
				dp[i][j]=dp[i-1][j-1]+1;
			}
			else{
				dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
			}
		}
	}

        计数dp，由元素和一定入手，考虑方案数的累加。

---

P1435 [IOI2000]回文字串 / [蓝桥杯2016省]密码脱落

        整出两个一正一反的序列，然后找最长公共子序列。

		sort(a+1,a+n+1);
		int ans=0;f[0]=1;
		for(int i=1;i<=n;i++){
			ans+=f[a[i]];    //选a[i]的方案数
			for(int j=sum;j>=a[i];j--){
				f[j]+=f[j-a[i]];    //不选a[i]的方案数
			}	
		}

---

P1853 投资的最大效益

类似于“纪念品”这道题。

	for(int i=1;i<=n;i++){
		for(int j=1;j<=d;j++){
			for(int k=a[j];k<=s;k++)
				dp[k]=max(dp[k],dp[k-a[j]]+b[j]);
		}
		s+=dp[s];
	}