poj1722 SUBTRACT

有一种非常直观且显然的dp做法，就是用$f(i,j,k)$维护区间$[i,j]$是否可以合并出$k$。但是这样的复杂度为$O(20000n^2)$，而且转移的复杂度为$O(20000^2)$，当然是不能用的。
既然不能带括号dp，那么直接把括号去掉，剩下的只是一个包含$n-1$个+或-的序列且第一个符号肯定是-。而把任何一个这样的序列还原成带括号的样子是很简单的，只需要遍历序列两次，第一次对所有+号进行操作，第二次对所有-号进行操作即可。
那么就可以直接对这个+和-的序列进行dp，$f(i,k)$表示前$i$个数能否组合成$k$的值，转移复杂度为$O(1)$。用$O(20000n)$的复杂度求出序列，还原括号时用树状数组维护一下消除了那些位置的数即可。

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
#define clr(a) memset(a,0,sizeof(a))
#define md 10000
typedef long long ll;
int _t,ar[100],ts[100],n;
bool dp[100][20010];
int te[110];
inline int _lw(int x){
    return x&(-x);
};
void upd(int d){
    for(;d<=n;d+=_lw(d))
        te[d]++;
};
int qry(int d){
    int r=0;
    for(;d;d-=_lw(d))
        r+=te[d];
    return r;
};
bool cl(){
    int i,j,k;
    if(scanf("%d %d",&n,&_t)==-1)return 0;
    for(i=0;i<n;scanf("%d",&ar[i++]));
    if(n==1)return 1;
    else if(n==2)
        printf("1\n");
    else{
        clr(dp);
        dp[1][ar[0]-ar[1]+md]=1;
        for(i=2;i<n;++i){
            for(j=0;j<=20000;++j){
                if(j+ar[i]<=20000)
                    dp[i][j]|=dp[i-1][j+ar[i]];
                if(j-ar[i]>=0)
                    dp[i][j]|=dp[i-1][j-ar[i]];
            }
        }
        for(i=n-1,j=_t+md,ts[0]=0;i>1;--i){
            if(j+ar[i]<=20000&&dp[i-1][j+ar[i]])
                ts[i-1]=0,j+=ar[i];
            else{
                ts[i-1]=1,j-=ar[i];
            }
        }
        clr(te);
        for(i=0;i<n-1;++i)if(ts[i])
            printf("%d\n",i+1-qry(i+1)),upd(i+2);
        for(i=0;i<n-1;++i)if(!ts[i])
            printf("%d\n",i+1-qry(i+1)),upd(i+2);
    }
    return 1;
};
int main(){
#ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
#endif
    while(cl());
    return 0;
};