hihocoder 1044 状态压缩dp

最新推荐文章于 2018-02-01 21:40:00 发布

点羽染墨

最新推荐文章于 2018-02-01 21:40:00 发布

阅读量716

点赞数

CC 4.0 BY-SA版权

文章标签：状态压缩dp

本文链接：https://blog.youkuaiyun.com/qq_24667639/article/details/45417551

本文介绍了一种使用C++实现的动态规划算法，通过定义状态转移方程解决特定问题。该算法涉及位运算和状态压缩技巧，适用于解决具有重叠子问题和最优子结构特征的问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>
#include <climits>
#include <cstring>
#include <cmath>
#include <map>
#include <set>
#define INF 100000000

using namespace std;
int n,m,q;
int dp[1050][1100];
int num[1100];
int a[1050];

int f(int i){
    int ans = 0;
    while(i){
        ans += i&1;
    //  cout <<"a " <<  i << endl;
    //  cout << ( i^0 )<< endl;
        i >>= 1;
    }
    return ans;
}

void fun(){
    for(int i = 0;i <= (1<<10);i++){
        num[i] = f(i); 
    }
}

int main(){
    fun();
    while(cin >> n >> m >> q){

        for(int i = 0;i < n;i++){
            scanf("%d",&a[i]);  
        }
        memset(dp[0],0,sizeof(dp[0]));

        for(int i = 0;i < n;i++){
            for(int s = 0; s < (1 << m);s++){
                if(num[s] <= q){
                    if(s&1){
                        dp[i+1][s] = a[i];
                        dp[i+1][s] += max(dp[i][(s>>1)^(1<<(m-1))],dp[i][(s>>1)]);                  
                    }
                    else{
                        if(num[s] == q){
                            dp[i+1][s] = dp[i][s>>1];
                        }
                        else{
                            dp[i+1][s] = max(dp[i][(s>>1)^(1<<(m-1))],dp[i][(s>>1)]);
                        }
                    }
                }
                else{
                    dp[i][s] = 0;
                }
            }
        }
        int max = 0;
        for(int i = 0;i < (1<<m);i++){
            if(num[i] <= q && dp[n][i] > max){
                max = dp[n][i];
            }
        }
        printf("%d\n",max);
    }
    return 0;
}