• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            強烈推薦此題!
            一種巧妙的DP,稱為插頭DP。每個格子是一個階段,狀態是通過這個格子的右邊一個格子的左邊和上邊的S形折線(共由m+1條線段組成)上的3進制數。

            這是題目的第一個sample,我把行、列都從1開始編號,原始圖設為g,紅線表示狀態f[2][2][(100221)3]。
            狀態設計好了,狀態轉移也就不是難事了。根據g[i][j + 1]的值來確定狀態轉移的方式。
            我是用順推的轉移,不大好表述,不懂的話看看代碼,再看看圖,想想就能明白。

            /*************************************************************************
            Author: WHU_GCC
            Created Time: 2007-9-24 19:35:22
            File Name: pku3133.cpp
            Description: 
            ***********************************************************************
            */

            #include 
            <iostream>
            using namespace std;

            #define out(x) (cout << #x << ": " << x << endl)
            typedef 
            long long int64;
            const int maxint = 0x7FFFFFFF;
            const int64 maxint64 = 0x7FFFFFFFFFFFFFFFLL;
            template 
            <class T> void show(T a, int n) for (int i = 0; i < n; ++i) cout << a[i] << ' '; cout << endl; }
            template 
            <class T> void show(T a, int r, int l) for (int i = 0; i < r; ++i) show(a[i], l); cout << endl; }

            const int maxn = 10;
            const int mask[11= {1392781243729218765611968359049};

            int g[maxn][maxn];
            char f[maxn + 1][maxn][59049];
            int n, m;
            int pow3m;

            inline 
            int get_bit(int x, int k)
            {
                
            return x % mask[k + 1/ mask[k];
            }


            inline 
            int set_bit(int x, int k)
            {
                
            return x - get_bit(x, k) * mask[k] - get_bit(x, k + 1* mask[k + 1];
            }


            int dp()
            {
                pow3m 
            = mask[m + 1];
                memset(f, 
            0x7Fsizeof(f));
                f[
            1][0][0= 0;
                
            for (int i = 1; i <= n; i++)
                
            {
                    
            for (int j = 0; j < m; j++)
                        
            for (int k = 0; k < pow3m; k++if (f[i][j][k] != 0x7F)
                        
            {
                            
            int bit_j = get_bit(k, j);
                            
            int bit_j1 = get_bit(k , j + 1);
                            
            int pre_bit = set_bit(k, j);
                            
            if (g[i][j + 1== 1)
                            
            {
                                
            if (bit_j == 0 && bit_j1 == 0)
                                    f[i][j 
            + 1][k] <?= f[i][j][k];
                            }

                            
            else if (g[i][j + 1== 2)
                            
            {
                                
            if (bit_j + bit_j1 == 1)
                                    f[i][j 
            + 1][pre_bit] <?= f[i][j][k] + 1;
                                
            if (bit_j + bit_j1 == 0)
                                
            {
                                    f[i][j 
            + 1][pre_bit + mask[j]] <?= f[i][j][k] + 1;
                                    f[i][j 
            + 1][pre_bit + mask[j + 1]] <?= f[i][j][k] + 1;
                                }

                            }

                            
            else if (g[i][j + 1== 3)
                            
            {
                                
            if (bit_j == 0 && bit_j1 == 2 || bit_j == 2 && bit_j1 == 0)
                                    f[i][j 
            + 1][pre_bit] <?= f[i][j][k] + 1;
                                
            if (bit_j + bit_j1 == 0)
                                
            {
                                    f[i][j 
            + 1][pre_bit + 2 * mask[j]] <?= f[i][j][k] + 1;
                                    f[i][j 
            + 1][pre_bit + 2 * mask[j + 1]] <?= f[i][j][k] + 1;
                                }

                            }

                            
            else
                            
            {
                                
            if (bit_j == 1 && bit_j1 == 0 || bit_j == 0 && bit_j1 == 1)
                                
            {
                                    f[i][j 
            + 1][pre_bit + mask[j]] <?= f[i][j][k] + 1;
                                    f[i][j 
            + 1][pre_bit + mask[j + 1]] <?= f[i][j][k] + 1;
                                }

                                
            if (bit_j == 0 && bit_j1 == 2 || bit_j == 2 && bit_j1 == 0)
                                
            {
                                    f[i][j 
            + 1][pre_bit + 2 * mask[j]] <?= f[i][j][k] + 1;
                                    f[i][j 
            + 1][pre_bit + 2 * mask[j + 1]] <?= f[i][j][k] + 1;
                                }

                                
            if (bit_j == 1 && bit_j1 == 1 || bit_j + bit_j1 == 4)
                                    f[i][j 
            + 1][pre_bit] <?= f[i][j][k] + 1;
                                
            if (bit_j + bit_j1 == 0)
                                
            {
                                    f[i][j 
            + 1][pre_bit] <?= f[i][j][k];
                                    f[i][j 
            + 1][pre_bit + mask[j] + mask[j + 1]] <?= f[i][j][k] + 1;
                                    f[i][j 
            + 1][pre_bit + 2 * mask[j] + 2 * mask[j + 1]] <?= f[i][j][k] + 1;
                                }

                            }

                        }

                    
            for (int k = 0; k < pow3m / 3; k++)
                        f[i 
            + 1][0][k * 3= f[i][m][k];
                }

                
            if (f[n][m][0== 0x7F)
                    
            return 0;
                
            else
                    
            return f[n][m][0- 2;
            }


            int main()
            {
                
            while (scanf("%d%d"&n, &m), n + m != 0)
                
            {
                    
            for (int i = 1; i <= n; i++)
                        
            for (int j = 1; j <= m; j++)
                            scanf(
            "%d"&g[i][j]);
                    printf(
            "%d\n", dp());
                }

                
            return 0;
            }
            posted on 2007-09-24 21:12 Felicia 閱讀(1162) 評論(4)  編輯 收藏 引用 所屬分類: 動態規劃
            Comments
            • # re: [動態規劃]pku3133
              Run&Run
              Posted @ 2007-11-24 09:19
              0,1,2 分別表示什么狀態啊?  回復  更多評論   
            • # re: [動態規劃]pku3133
              Felicia
              Posted @ 2007-11-24 15:26
              看圖
              0表示沒有線穿過,1表示灰線穿過,2表示藍線穿過  回復  更多評論   
            • # re: [動態規劃]pku3133
              longshen
              Posted @ 2008-11-11 09:15
              看后對連通性狀態壓縮DP有點開竅了...

              謝謝!  回復  更多評論   
            • # re: [動態規劃]pku3133
              chuliuxiang
              Posted @ 2009-10-07 09:26
              還是有些不懂!  回復  更多評論   
             
            久久久久久av无码免费看大片| 国产高清美女一级a毛片久久w| 九九精品99久久久香蕉| 欧美精品福利视频一区二区三区久久久精品 | 久久久久久国产精品无码下载| 波多野结衣久久| 99精品久久久久久久婷婷 | 久久综合综合久久狠狠狠97色88 | 潮喷大喷水系列无码久久精品 | 99久久这里只精品国产免费 | 青青国产成人久久91网| 亚洲国产精品久久久久婷婷软件 | 久久成人影院精品777| 久久无码一区二区三区少妇| 久久久久波多野结衣高潮| 91精品国产综合久久久久久 | 蜜臀av性久久久久蜜臀aⅴ麻豆 | 亚洲va久久久久| 久久中文娱乐网| 狠狠色丁香婷婷久久综合| 97久久精品无码一区二区天美| 国产成人精品久久亚洲高清不卡| 日本五月天婷久久网站| 亚洲国产精品久久66| 亚洲精品白浆高清久久久久久| 日本亚洲色大成网站WWW久久 | 一本大道久久a久久精品综合| 亚洲综合久久夜AV | 久久最新免费视频| 久久er99热精品一区二区| 人妻精品久久无码区| 久久综合九色综合欧美就去吻| 99久久精品日本一区二区免费 | 精品无码久久久久久国产| 国产精品久久久久天天影视| 欧美粉嫩小泬久久久久久久| 99久久精品免费观看国产| 精品久久久久久成人AV| 亚洲va久久久噜噜噜久久男同| 亚州日韩精品专区久久久| 久久久久人妻精品一区三寸蜜桃|