动态规划---51nod1371

来源：互联网发布：红外相机软件编辑：程序博客网时间：2024/06/03 03:20

这一题考察的是完全状态的获取能力:--->好题,回头一定要再复习一遍.主要是复习思想.
DP状态的获得是由小状态通过找出后效型的问题关键然后再一步一步增加限制的.所以一开始的dp[i]很重要.往往先不出来的时候可以开绿一下是不是dp[i]的状态是否正确.注意细节-----C(n , 2) ;

#include <stdio.h>#include <string>#include <string.h>#include <queue>#include <stack>#include <map>#include <iostream>#include <stdlib.h>#include <math.h>#include <set>#include <algorithm>#define inf 0x3f3f3f3f#define mem0(x , y)  memset(x , y , sizeof(x))#define ll long long#define rep(x , y)   for(int x=0;x<y;x++)#define per(x , y)   for(int x=y-1;x>=0;x--)#define lowbit(x) (x & (-x))#define read(x) scanf("%d",&x)#define print(x) printf("%d\n" , x) ;#define middle int mid = (l + r)/2#define lson rt<<1 , l , mid#define rson rt<<1|1 , mid+1, r#define CASE int T ; int ca = 0 ; read(T) ; while(T--)#define RT tree[rt]#define FP freopen("1" , "r" , stdin)#define mem0(x , y) memset(x , y , sizeof(x))#define BIN(x) ll bin ; bin = (ll)1 << x ; ///2的x次方#define E Edge[i]const int MAXN = 500000 ;const int MAXV = 200000 ;const int MAXE = 200000 ;const int mod  = 1e8+7  ;using namespace std ;/****************前向星****************struct edge{    int s , e , w , next ;    bool operator < (const edge & e) const{        return w < e.w ;    }}Edge[MAXE] ;int head[MAXV] , hct = 0 ;void add_edge(int s,int e,int w){    Edge[hct].s = s ; Edge[hct].e =e ;Edge[hct].w = w ; Edge[hct].next = head[s] ;    head[s] = hct ++ ;}/*************************************/#define MOD dp[i][j][k] %= mod#define DP  dp[i][j][k]ll dp[210][210][210] ;ll ans = 0 ;map <ll , int> maps ;int main(){    maps [1] = 10 ;    for(int i=0;i<10000;i++){        for(int j=0;j<10000;j++){            maps[i] = 10 ;        }    }    return 0 ;    int n ; read(n) ;    dp[0][0][0] = 1  ;    for(int i=1;i<=n;i++){        for(int j=0;j<=i;j++){            for(int k=0;k+j<=i;k++){                if(k >= 1) DP += dp[i-1][j][k-1] * (i-j-k+1); MOD ;                if(j >= 1) DP += dp[i-1][j-1][k] * (i-j-k+1); MOD ;                if(k >= 1) DP += dp[i-1][j+1][k-1] * (j+1 ) ; MOD ;                if(j >= 2) DP += dp[i-1][j-2][k] * (i-j-k+2)*(i-j-k+1)/2 ; MOD ;                if(k >= 1) DP += dp[i-1][j][k-1]*(i-j-k+1)*j; MOD ;                if(k >= 2) DP += dp[i-1][j+2][k-2] * (j+2)*(j+1)/2  ; MOD ;                DP += dp[i-1][j][k] ; MOD ;                if(i == n) {                    ans += DP  ;                    ans %= mod ;                }            }        }    }    cout << ans << endl ;}

0 0