HRBUST 1846 方格涂色 （DP+递推）

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))

#define pii pair<ll,ll>
#define mk(x,y) make_pair(x,y)
#define fi first
#define se second

const int  maxn =1e5+5;
const int mod=1e9+7;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
主体思路还是dp+递推。
正难则反，我们用dp(i)表示前i个
中连续个数不超过m个的方案数，
那么我们只要枚举当前位置连续段长度就可以得到递推式。
这样维护个dp的前缀和就方便推的下一位的dp值。
复杂度：O(n).
*/

ll dp[maxn],n,m,pw[maxn];

int main(){
    pw[0]=1LL;rep(i,1,maxn) pw[i]=2*pw[i-1]%mod;
    while(scanf("%lld%lld",&n,&m)!=EOF){
        rep(i,1,m){
            dp[i]=powmod(2LL,i);
            dp[i]=(dp[i]+dp[i-1])%mod;
        }
        dp[m]=(2*dp[m-1]%mod+1)%mod;///再次初始化
        rep(i,m+1,n+1){
            dp[i]=(dp[i-1]-dp[i-m-1]+mod)%mod;
            dp[i]=(dp[i]+dp[i-1])%mod;
        }
        ll ans=(dp[n]-dp[n-1]+mod)%mod;
        ans=(pw[n]-ans+mod)%mod;
        printf("%lld\n",ans);
    }
    return 0;
}