状压dp（3）

vijos严厉的班长

木姑娘在班级里面是班长。虽然是副班长，却有着比正班长更高的威信，并深受小朋友们的爱戴。
每天眼保健操时间，木姑娘都要监督所有小朋友认真做眼保健操。整个过程被描述为n个时间段，第i个时间段内对于揉捏穴位的速度被要求为a[i]。
如果某一位小朋友在n个时间段内，速度分别为b[1],b[2],…,b[n]，则他的最终得分为 score = |a[1]-b[1]| + |a[2]-b[2]| + … + |a[n]-b[n]|。
分数越大，则说明表现越差，如果b[1]到b[n]中，又满足任意两个数字均互素，则会受到木姑娘的严厉批评。
对于爱慕木姑娘的孩子来说，能和木姑娘有更多时间待在一起，才是最关键的。

两个数字不能有相同质因子。30-1是29 1和所有互素，怎么都可以选，，，那59以外的选不到，每个质因子在所有选的数质因子分解之后只能出现一次，那么dp[i][s]表示选了i个b了，s：59以内质因子就16个，所以0/1表示选没选，状压，

当前i选什么数字，就要1~59跑一遍，看看谁可以（预处理所有的质因数分解开的状态），然后对每个状态，取答案的最小值

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
int prime[20] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};
const int m = 16;
const int maxn = 110;
int dp[maxn][(1<<m)+10];
int n,a[maxn],bi[60];

void print(int tmp){
    cout <<tmp <<endl;
    while(tmp){
        printf("%d ",tmp%2);
        tmp /= 2;
    }
    printf("\n");
}

void init(){
    scanf("%d",&n);
    for(int i = 1; i <= n ; i++)
        scanf("%d",&a[i]);

    for(int i = 2; i < 59 ; i++){
        for(int j = 0 ; j < 16 && i >= prime[j] ; j++){
            if(i % prime[j] == 0)
                bi[i] |= (1 << j);
        }
       // printf("bi[%d] = ",i);
        //print(bi[i]);
    }

}

void sov(){
    for(int i = 0; i <= n ; i++)
        for(int j =  0 ; j < (1<<m) ;j++)
            dp[i][j] = (1<<30);
    dp[0][0] = 0;
    int ans = (1<<30);
    for(int i = 1; i <= n ; i++){
        for(int pre = 0 ; pre < (1<<m) ; pre++){
            //1:
            dp[i][pre] = min(dp[i][pre],dp[i-1][pre]+a[i]-1);
          //  cout << "dp [ "<<i <<" ] [ "<<pre<<" ] = "<<dp[i][pre]<<endl;
            if(i == n)  ans = min(ans,dp[i][pre]);
            if(a[i] == 1)  continue;

            for(int now = 2 ; now <= 59; now++){
               // if(now == 5 && pre == 0 && i == 2)  printf("dp[i-1][pre] = %d\n",pre&bi[now]);
                if((pre&bi[now])==0){
                    dp[i][pre|bi[now]] = min(dp[i][pre|bi[now]] ,dp[i-1][pre]+max(now-a[i],a[i]-now));
                    if(i == n)
                        ans = min(ans,dp[i][pre|bi[now]]);
                }
            //qif(now == 5 && pre == 0 && i == 2)  printf("dp[i-1][pre] = %d\n",dp[i][pre|bi[now]]);
            }
        }

    }
    printf("%d\n",ans);
}

int main(){
    init();
    sov();
}

ural1326. Bottle Taps

http://acm.timus.ru/problem.aspx?space=1&num=1326
给你一个n代表要从n个物品中买东西，然后告诉你这n个东西的单价，在给你m个集合的情况，就是每个结合中有x件物品，他们合起来买的价格是k。这x件物品依次是：p1……px。之后给你一个kk，表示你要买的物品的编号。让你求出来如何花费最少的钱买到要求的序列。
状压：s表示每个买了还是没买

#include <cstdio>

#include <iostream>

#include <cstring>

using namespace std;

int k ,p,t,s[(1<<20)+10];

const int maxn = 125;

int dp[(1<<20)+10],m,n,cost[maxn];

void init(){

    for(int i = 1; i <= n ; i++){

        scanf("%d",&cost[i]);

        s[i] = (1<<(i-1));

    }

    scanf("%d",&m);

    for(int i = 1; i <= m ; i++){

        scanf("%d",&cost[i+n]);

        scanf("%d",&k);

        for(int j = 1; j <= k ; j++){

            scanf("%d",&p);

            s[i+n] |= (1<<(p-1));

        }//cout << " i = "<<i <<"  s  = "<<s[i+n]<<endl;

    }

    scanf("%d",&k);

    for(int i = 1; i <= k ; i++){

        scanf("%d",&p);

        t |= (1<<(p-1));

       // cout <<" t =  "<<t<<endl;

    }

}



void sov(){

    for(int i = 0 ; i < (1<<n) ; i++)

            dp[i] = (1<<30);

    dp[0] = 0;

    int ans = (1<<30);

    for(int i = 0; i < (1<<n); i++){

       // cout <<"dp  "<< i <<endl;

        for(int j = 1; j <= n+m; j++){

            //printf("i = %d  s[%d] = %d  i|s[j] = %d dp[i|s[j]] = %d  dp[i] = %d  cost  = %d\n",i,j,s[j],i|s[j],dp[i|s[j]], dp[i], cost[j]);

            dp[i|s[j]] = min(dp[i|s[j]],dp[i]+cost[j]);

           // cout << "dp  ["<<(i|s[j])<<" ] = "<<dp[i|s[j]]<<endl;

           // if(j == 7)  cout <<"i|s[j] = "<<(i|s[j])<<endl;

            if( ((i|s[j]) & t) == t){

               // cout <<" i = "<<i<<endl;

                ans = min(ans,dp[i|s[j]]);

            }

        }

    }

    printf("%d\n",ans);

}



int main()`


    scanf("%d",&n);

    init();

    sov();

}

ural1152. False Mirrors

http://acm.timus.ru/problem.aspx?space=1&num=1152
有n个阳台，首尾相接形成一个环，阳台里有怪物，伤害为a[i]，你每次可以打掉连续的三个阳台，那么同时，没被打掉的那些怪物会对你造成响应a[i]的伤害。经过几次战斗，你打掉所有阳台的怪物，问受到的伤害最小是多少。
觉得区间dp也可以写啊？

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
const int m = 20;
int dp[(1<<m)+10];
int n,a[m];
int sum(int s){
    int sum = 0;
    for(int i = 0 ; i < n ; i++){
        if(s & (1 << i))
            sum += a[i];
    }
    return sum;
}
int dfs(int s){
    if(dp[s] != -1)   return dp[s];
    if(s == 0)   return dp[s] = 0;

    int cnt = (1<<30);
    for(int i = 0; i < n ; i++){

        if(s&(1<<i)){
            int next = s & ( ~(1<<((i-1+n)%n))) ;
            next &= ~(1<<i);
            next &= ~(1<<((i+1)%n));

            cnt = min(cnt,dfs(next)+sum(next));
        }
    }
    return dp[s] = cnt;
}

void init(){
    for(int i = 0 ; i < n ; i++){
        scanf("%d",&a[i]);
    }
    cout << dfs((1<<n)-1) <<endl;
}

int main(){
    while(~scanf("%d",&n)){
        memset(dp,-1,sizeof(dp));
        init();
    }
}