luogu1941 [NOIp2014]飞扬的小鸟 (dp)

设f[i][j]为到达(i,j)这个位置的最小操作数

就有$f[i][j]=min{f[i-1][j+Y[i-1]],f[i-1][j-X[i-1]*k]+k}$

然后考虑优化一下转移：

对于一系列模x[i-1]相同的高度，它们都可以转移到模x[i-1]相同的高度、而且在它们上边的点，所以只要从下往上做，不断地取一个最小值就可以了（会意 会意...）

要注意的是不管怎么操作，只要操作完高度>M都会变成=M，而且还可以从M点一下还维持在M

#include<bits/stdc++.h>
#define pa pair<int,int>
#define lowb(x) ((x)&(-(x)))
#define REP(i,n0,n) for(i=n0;i<=n;i++)
#define PER(i,n0,n) for(i=n;i>=n0;i--)
#define MAX(a,b) ((a>b)?a:b)
#define MIN(a,b) ((a<b)?a:b)
#define CLR(a,x) memset(a,x,sizeof(a))
#define rei register int
using namespace std;
typedef long long ll;
const int maxn=10010,maxm=1010;

inline ll rd(){
    ll x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*neg;
}

int N,M,K;
int x[maxn],y[maxn];
int f[2][maxm],l[maxn],h[maxn],ans=1,cnt;
bool istube[maxn];

int main(){
    // freopen("testdata.in","r",stdin);
    rei i,j,k;
    N=rd(),M=rd(),K=rd();
    for(i=0;i<N;i++) x[i]=rd(),y[i]=rd();
    for(i=0;i<=N;i++) l[i]=1,h[i]=M*2;
    for(i=1;i<=K;i++){
        int a=rd(),b=rd(),c=rd();
        istube[a]=1;
        l[a]=b+1;h[a]=c-1;
    }
    CLR(f,127);CLR(f[0],0);
    bool b=1;
    for(i=1;i<=N&&ans;i++){
        CLR(f[b],127);
        for(j=l[i];j<=min(M,h[i])&&j+y[i-1]<=min(M,h[i-1]);j++){
            if(j+y[i-1]>=l[i-1]) f[b][j]=f[b^1][j+y[i-1]];
        }
        for(j=l[i-1];j<l[i-1]+x[i-1]&&j<=min(M,h[i-1]);j++){
            int mm=f[b^1][j];
            // printf("!%d %d
",j,mm);
            for(k=1;j+x[i-1]*k<=h[i];k++){
                int jk=j+x[i-1]*k;
                // printf("!!!%d %d
",jk,mm);
                bool re=0;
                if(jk>M) jk=M,re=1;
                if(jk>=l[i]) f[b][jk]=min(f[b][jk],mm+1);
                ++mm;mm=min(mm,f[b^1][jk]);
                if(re) break;
            }
        }
        for(ans=0,j=l[i];j<=h[i];j++){
            // printf("%d %d %d
",i,j,f[b][j]);
            if(f[b][j]<=1e8) {ans=1;break;}
        }
        if(ans&&istube[i]) cnt++;
        b^=1;
    }
    printf("%d
",ans);
    if(ans){
        ans=1e8;
        for(i=l[N];i<=min(M,h[N]);i++) ans=min(ans,f[b^1][i]);
        printf("%d
",ans);
    }else{
        printf("%d
",cnt);
    }
    return 0;
}