蓝桥杯 【dp?】.cpp

题意:

  给出一个2*n的方格,当刷完某一个方格的漆后可以且只可以走到相邻的任何一格,即上 下 左 右 左上 左下 右上 右下。可以从任意一个格子开始刷墙,问有多少种刷法,因为随着n的增大方案数会变多,因此输出方案数mod 1000000007.

 

思路:

  dp[][2], dp[i][0]表示2*i的格子从第i列开始刷,最后回到该格子下面

         dp[i][1]表示2*i的格子从第i列开始刷,最后无法回到该格子下面

  状态转移方程是:dp[i][0] = 2*dp[i-1][0] 即 A1->B1->...->B2->A2, A1->B2->...->B1->A2,  A2也一样

          dp[i][1] = 2*(dp[i-1][0]+dp[i-1][1]) + 4*(dp[i-2][0]+dp[i-2][1])

    因为不能回到那一端,所以 一种情况是 1- >2  然后是dp1[i-1]+dp2[i-1],或者2->1  dp1[i-1]+dp2[i-1],所以乘以2。

或者可以1->3->2->4,1->4->2->3,2->3->1->4,2->4->1->3,所以乘以4。

    初始状态:dp[0][0] = 1, dp[0][1] = 0, dp[1][0] = 2, dp[1][1] = 0;

       没有格子的情况下每个子涂所以回到该列只有一种情况,而回不到该列完全不可能。只有一列的情况下回到该列的方案数为2,即A1->A2, 或者A2->A1

  最后罗列出从第一列到最后一列的总共刷墙方案数。

  ans += 2*((dp[i-1][0]+dp[i-1][1])*dp[n-i][0])) + 2*((dp[n-i)[0]+dp[n-i][1])*dp[i-1][0])

        即从第i列开始刷,先刷右边的格子,有dp[n-i][0]种方案,因为要刷回来才可以刷左边的格子,然后再刷左边的格子,总共有(dp[i-1][0]+dp[i-1][1])种方案。同理可以先刷左边的格子再刷右边的格子

 

Tips:

  nothing?因为公式的推理比较重要..

 

Code:

 1 #include <stdio.h>
 2 #include <cstring>
 3 #include <algorithm>
 4 using namespace std;
 5 
 6 const int MAXN = 1010;
 7 const int mod = 1000000007;
 8 
 9 long long dp[MAXN][2];
10 long long ans[MAXN] = {0, 2};
11 
12 void init()
13 {
14     dp[0][0] = 1, dp[0][1] = 0;
15     dp[1][0] = 2, dp[1][1] = 0;
16     for (int i = 2; i < MAXN; ++i) {
17         dp[i][0] = 2*dp[i-1][0]%mod;
18         dp[i][1] = (2*(dp[i-1][0]%mod+dp[i-1][1]%mod)%mod + 4*(dp[i-2][0]%mod+dp[i-2][1]%mod)%mod)%mod;
19     }
20 
21     for (int i = 2; i < MAXN; ++i) {
22         ans[i] = (dp[i][0]%mod+dp[i][1]%mod)%mod;
23         ans[i] = (ans[i]*2)%mod;
24         for (int j = 2; j < i; ++j) {
25             ans[i] += (2*((dp[j-1][0]+dp[j-1][1])%mod*dp[i-j][0]%mod)%mod)%mod;
26             ans[i] += (2*((dp[i-j][0]+dp[i-j][1])%mod*dp[j-1][0]%mod)%mod)%mod;
27         }
28         ans[i] %= mod;
29     }
30 
31 }
32 
33 int main()
34 {
35     int n;
36     init();
37     while (~scanf("%d", &n)) {
38         printf("%lld\n", ans[n]);
39     }
40     return 0;
41 }
View Code

 

链接:http://acm.nbut.cn/Problem/view.xhtml?id=1476

posted @ 2013-08-11 14:05  Griselda.  阅读(342)  评论(0编辑  收藏  举报