TOJ 1684 Corn Fields -- 状态压缩DP

9 篇文章 0 订阅

题目链接:http://acm.tzc.edu.cn/acmhome/problemdetail.do?&method=showdetail&id=1684

题目大意:给定一个0-1矩阵,要求从矩阵中选出一些位置满足以下条件:

1.只能在值为1的位置上选

2.选出的位置不能相邻(选出位置(i, j)后,不能选(i, j)上下左右4个位置)

问有多少种选法。

分析:矩阵维度很小,直接把一行的状态压缩然后对列动态规划。记下来是因为一开始错了两遍:我用的是滚动数组,在求新值之前忘了把数组清空。顺便提一下对状态的check,如何简洁高效地判断状态是否满足上述两个条件。

对于条件1,把一行的值视为一个二进制数,先按位取反得到a,那么状态t满足条件1只需要(a & 1 == 0)。即相同位上状态与矩阵中的值不同时为1。

对于条件2,判断相邻行的两个状态t1, t2时,只需要(t1 & t2 == 0),判断同一行时,即判断状态t是否满足条件2,只需要(t & (t << 1) == 0)

我事先把每一行的合法状态存下来了,这样可以避免对无用状态的判断。

#include 
   
   
    
    
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
         #include 
         #include 
         
           #include 
          
            #include 
           
             #include 
            
              #include 
             
               #include 
              
                #include 
               
                 #include 
                
                  #include 
                 
                   #include 
                  
                    #define mp make_pair #define X first #define Y second #define MEMSET(a, b) memset(a, b, sizeof(a)) using namespace std; typedef unsigned int ui; typedef long long ll; typedef unsigned long long ull; typedef pair 
                   
                     pii; typedef vector 
                    
                      vi; typedef vi::iterator vi_it; typedef map 
                     
                       mii; typedef priority_queue 
                      
                        pqi; typedef priority_queue 
                       
                         , greater 
                        
                          > rpqi; typedef priority_queue 
                         
                           pqp; typedef priority_queue 
                          
                            , greater 
                           
                             > rpqp; const int MAX_SIZE = 12 + 2; const int MOD = 100000000; int a[MAX_SIZE]; int coll[MAX_SIZE][1 << MAX_SIZE]; int dp[2][1 << MAX_SIZE]; void compute_coll(int k, int m) { for (int i = 0; i < (1 << m); ++i) { if (!(i & a[k]) && !(i & (i << 1))) { coll[k][++coll[k][0]] = i; } } } int main(int argc, char *argv[]) { // freopen("D:\\in.txt", "r", stdin); int n, m, i, j, k; cin >> n >> m; for (i = 0; i < n; ++i) { for (j = 0; j < m; ++j) { int x; cin >> x; a[i] = a[i] + a[i] + (x ^ 1); } compute_coll(i, m); } for (i = 1; i <= coll[0][0]; ++i) { dp[0][i] = 1; } for (i = 1; i < n; ++i) { int id1 = i & 1, id2 = (i & 1) ^ 1; MEMSET(dp[id1], 0); for (j = 1; j <= coll[i][0]; ++j) { for (k = 1; k <= coll[i - 1][0]; ++k) { if (!(coll[i][j] & coll[i - 1][k])) { dp[id1][j] = (dp[id1][j] + dp[id2][k]) % MOD; } } } } int ans = 0, id = (n - 1) & 1; for (i = 1; i <= coll[n - 1][0]; ++i) { ans = (ans + dp[id][i]) % MOD; } cout << ans << endl; return 0; } 
                            
                           
                          
                         
                        
                       
                      
                     
                    
                   
                  
                 
                
               
              
             
            
           
          
       
      
      
     
     
    
    
   
   

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值