21训练 dp专题题解

传送门

A：观察图形可发现，第i格是由第i - 1格和第i - 2 格转移过来，即 f[i]=f[i-1]+f[i-2];
代码：

#include<iostream>

using namespace std;
typedef long long ll;

ll f[100];
int main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	f[1]=1,f[2]=1;
	for(int i=3;i<=50;i++) f[i]=f[i-1]+f[i-2];
	int t;
	cin>>t;
	while(t--){
		int n,m;
		cin>>n>>m;
		m-=n-1;
		cout<<f[m]<<"\n";
	}
}

B 由贪心可知 价格最大的最后买，这时要注意要想买价格最大的，必须要有5块钱 所以直接可以转化成背包问题 但是为了最后买价格最大的 总体积要提前 - 5，要注意一下如果钱数小于5 一个都不能买；
代码：

#include<bits/stdc++.h>

using namespace std;
int dp[1010],v[1010];
int main(){
	ios::sync_with_stdio(false);
	int n;
	while(cin>>n,n){
		for(int i=1;i<=n;i++) cin>>v[i];
		int m;
		cin>>m;
		sort(v+1,v+1+n);
		if(m<5){
			cout<<m<<"\n";
			continue;
		}
		for(int i=0;i<=m;i++) dp[i]=0;
		for(int i=1;i<n;i++){
			for(int j=m-5;j>=v[i];j--){
				dp[j]=max(dp[j],dp[j-v[i]]+v[i]);
			}
		}
		cout<<m-dp[m-5]-v[n]<<"\n";
	}
	
} 

C、直接从下往上遍历即可，也可以从上往下遍历。
代码：

#include<bits/stdc++.h>

using namespace std;
int a[1010][1010];
int main(){
	ios::sync_with_stdio(false);
	int t;
	cin>>t;
	while(t--){
		int n;
		cin>>n;
		for(int i=1;i<=n;i++){
			for(int j=1;j<=i;j++) cin>>a[i][j];
		}
		
		for(int i=n-1;i>=1;i--){
			for(int j=1;j<=i;j++){
				a[i][j]=max(a[i+1][j],a[i+1][j+1])+a[i][j];
			}
		}
		cout<<a[1][1]<<"\n";
		
	}
	
} 

D、经典导弹问题 思考由O(n^2)的写法为什么能转移到这个O(nlogn)的写法，这里直接上O(nlogn)的代码
代码：

#include<bits/stdc++.h>

using namespace std;
int a[1000010],q[1000010];
int main(){
	ios::sync_with_stdio(false);
	int n;
	while(cin>>n){
		for(int i=1;i<=n;i++) cin>>a[i];
		q[1]=a[1];
		int len=1;
		for(int i=1;i<=n;i++){
			if(q[len]<a[i]) q[++len]=a[i];
			else{
				int l=1,r=len;
				while(l<r){
					int mid=l+r>>1;
					if(q[mid]>=a[i]) r=mid;
					else l=mid+1;
				}
				q[l]=a[i];
			}
		}
		cout<<len<<"\n";
	}
	
} 

E 搞不懂这个题为啥没人写 ，代码略，再看看能写不

F、第一秒可以拿到 位置4 5 6的馅饼 所以初始化就是

        f[1][4]=a[1][4];
        f[1][5]=a[1][5];
        f[1][6]=a[1][6];

因为每一秒都可以拿到上一秒的 第 i 或者 i - 1或者i - 2位置的馅饼，再处理一下0 和10的位置，状态转移方程就是：

if(j>0) f[i][j]=max(f[i-1][j],f[i-1][j-1])+a[i][j];
if(j<10) f[i][j]=max(f[i][j],max(f[i-1][j],f[i-1][j+1])+a[i][j]);

代码：

#include<bits/stdc++.h>

using namespace std;
typedef long long ll;
const int N=100010;
int f[N][12],a[N][12];


int main(){
	ios::sync_with_stdio(false);
	
	int n;
    while(cin>>n,n!=0){
        memset(a,0,sizeof a);
        memset(f,0,sizeof f);
        
        for(int i=1;i<=n;i++){
            int t,x;
            cin>>x>>t;
            a[t][x]++;
        }
        f[1][4]=a[1][4];
        f[1][5]=a[1][5];
        f[1][6]=a[1][6];
        for(int i=2;i<=100000;i++){
            for(int j=0;j<=10;j++){
                if(j>0) f[i][j]=max(f[i-1][j],f[i-1][j-1])+a[i][j];
                if(j<10) f[i][j]=max(f[i][j],max(f[i-1][j],f[i-1][j+1])+a[i][j]);
            }
        }
        int res=0;
        for(int i=0;i<=10;i++) res=max(res,f[100000][i]);
        cout<<res<<"\n";
    }
}

G、可以先预处理出来 只管第 n - 1和第 n 位的关系（即第n位的颜色不等于第n-1位的颜色） 不管第n位和第 1 位的关系（即第n位的颜色不等于第n-1位的颜色），每次求答案时候 在处理第n位和第 1 位的关系
dp[i][j][k]表示第一位颜色是i 第 j 位颜色为k的方案数
代码：

#include<iostream>
#include<algorithm>

using namespace std;
typedef long long ll;

ll dp[4][60][4];

int main() {
	ios::sync_with_stdio(false);
	int n;
	dp[1][1][1]=1;
	dp[2][1][2]=1;
	dp[3][1][3]=1;
	for(int i=2; i<=50; i++) {
		for(int j=1; j<=3; j++) {
			for(int k=1; k<=3; k++) {
				for(int z=1; z<=3; z++) {
					if(z==k) continue;
					dp[j][i][k]+=dp[j][i-1][z];
				}
			}
		}
	}
	while(cin>>n) {


		ll res=0;
		for(int i=1; i<=3; i++) {
			for(int j=1; j<=3; j++) {
				if(i==j) continue;
				res+=dp[i][n][j];
			}
		}
		if(n==1) res=3;
		cout<<res<<"\n";
	}
}

J、多重背包模板题，不再进行分析，直接上代码：

#include<bits/stdc++.h>

using namespace std;

int n;

typedef long long ll;
ll w[100010],v[100010],f[100010];
void Test(){
	memset(f,0,sizeof f);
	ll a,b,m,num;
	cin>>n>>m;
	int cnt=1;
	for(int i=1;i<=m;i++){
		cin>>a>>b>>num;
		ll ans=1;
		while(num>=ans){
			v[cnt]=ans*a;
			w[cnt++]=ans*b;
			num-=ans;
			ans*=2;
		}
		if(num!=0) v[cnt]=num*a,w[cnt++]=num*b;
	}
	for(int i=1;i<cnt;i++){
		for(int j=n;j>=v[i];j--)
		f[j]=max(f[j],f[j-v[i]]+w[i]);
	}
	
	cout<<f[n]<<endl;;
}

int main(){
	ios::sync_with_stdio(false);
	int t;
	cin>>t;
	while(t--)
	Test();
}

K、数位dp模板题 这里可以利用记忆化搜索来完成dp，直接上代码
代码：

#include<bits/stdc++.h>

using namespace std;
typedef long long ll;
ll dp[100][100];
int a[30];

ll dfs(int u,int per,bool flag){
	if(u==0) return 1;
	if(!flag&&dp[u][per]!=-1) return dp[u][per];
	int x;
	if(flag) x=a[u];
	else x=9;
	ll ans=0;
	for(int i=0;i<=x;i++){
		if(per==4&&i==9) continue;
		ans+=dfs(u-1,i,flag&&i==a[u]);
	}
	if(!flag) dp[u][per]=ans;
	return ans;
}


ll f(ll n){
	int cnt=0;
	while(n){
		a[++cnt]=n%10;
		n/=10;
	}
	return dfs(cnt,-1,1);
}

int main(){
	ios::sync_with_stdio(false);
	memset(dp,-1,sizeof dp);
	int t;
	cin>>t;
	while(t--){
		ll n;
		cin>>n;
		cout<<n-f(n)+1<<"\n";	
	}
	
}