2017 湘潭邀请赛&JSCPC G&J

训练的时候对G想了一个假算法。。也有很大可能是写错了。。

下来一看别人的G 看起来很奇妙。。

开始把所有的左括号翻成右括号，然后cost*=-1 这样在优先队列中就是最优的 然后for每一段 如果前缀和小于0就从优先队列中取右括号翻转

最后的结果一定是一个可行的括号序列

1 如果所有的原左括号都被重新反转，那么我们的选择仍然是最优，因为实际上没有在这个上面消耗

2 如果原左括号没有都被反转回去 那么我们的实际消耗也是最少的

J题一看就能想出来nmmm的想法 但是很远。。并没有想到nmm的 没想到可以nmmlogm过 就比较神奇了

画图可以发现 数字的大小其实是从右向左不断扩散的

思考nmmm的做法 对于第i个数字 枚举他是第k个a中的数字，那么就需要找到 第j个数字让它在a中做第k+1位来承接状态

但是这样状态不全 第j个数字是有限制的 它需要比第i个数字大或者小 而又有范围

因为数字的大小有规律 所以把范围也记录下来 第i个数字一定是范围的一边 于是只需要记录另一边就可以了

对于“另一边” 要么是第j个数字要么是之前传下来的 总之我们需要for一下这个范围来让他们加上当前状态

考虑区间加直接把状态加上去就可以把m变成logm了。。虽然复杂度仍然感觉不科学

G

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<vector>
#include<queue>
#include<map>
#include<algorithm>
#include<stack>
using namespace std;
#define L long long
 
#define rep(i,l,r) for(L i = l ; i <= r ; i ++ )
 
const L mod = 1000000000 + 7 ;
 
struct node {
    L cost ;
    L num ;
    friend bool operator < (node A , node B) {
        return A.cost > B.cost ;
    }
};
 
L a[100050] ;
L b[100050] ;
 
int main () {
    L n ;
    while(scanf("%lld" , &n) != EOF) {
        L ans = 0 ;
        rep(i,1,n) {
            L l,d ; char s[20] ;
            scanf("%lld%s%lld" , &l,s,&d) ;
            a[i] = l ;
            b[i] = d ;
            if(s[0] == '(') {
                ans += a[i] * b[i] ;
                b[i] *= -1 ;
            }
        }
        priority_queue<node> q ;
        L sum = 0 ;
        rep(i,1,n) {
            node c ;
            c.num = a[i] ;
            c.cost = b[i] ;
            q.push(c) ;
            sum -= c.num ;
            if(sum >= 0) continue ;
            L ned = sum * (-1) ;
            ned ++ ; ned /= 2 ;
            sum += ned * 2;
            while(ned > 0) {
                node c = q.top() ; q.pop() ;
                if(c.num >= ned) {
                    c.num -= ned ;
                    ans += ned * c.cost ;
                    q.push(c) ;
                    ned = 0 ;
                }
                else {
                    ned -= c.num ;
                    ans += c.num * c.cost ;
                }
            }
        }
        printf("%lld\n" , ans) ;
    }
}

J

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<vector>
#include<queue>
#include<map>
#include<algorithm>
#include<stack>
using namespace std;
#define L long long
 
#define rep(i,l,r) for(L i = l ; i <= r ; i ++ )
 
const L mod = 1000000000 + 7 ;
 
L n , m ;
 
L a[22] ;
L b[505] ;
 
L dp[505][505][22] ;
L s[505][505][22] ;
 
L lowbit(L x) {
    return (x & (-x)) ;
}
 
void add(L x,L j,L k,L val) {
    while(x <= m) {
        s[x][j][k] += val ;
        while(s[x][j][k] < 0) s[x][j][k] += mod ;
        s[x][j][k] %= mod ;
        x += lowbit(x) ;
    }
}
L sum(L x,L j,L k) {
    L sum = 0 ;
    while(x) {
        sum += s[x][j][k] ;
        sum %= mod ;
        x -= lowbit(x) ;
    }
    return sum ;
}
void upda(L l,L r,L val,L j,L k) {
    add(l,j,k,val) ;
    add(r+1,j,k,-val) ;
}
 
int main () {
    while(cin >> n >> m) {
        rep(i,1,n) cin >> a[i] ;
        rep(i,1,m) cin >> b[i] ;
        memset(dp, 0, sizeof(dp)) ;
        memset(s, 0, sizeof(s)) ;
        L ans = 0 ;
        for(L i = 1 ; i <= m ; i ++ ) {
            for(L j = 1 ; j <= m ; j ++ ) {
                for(L k = 1 ; k <= n ; k ++ ) {
                    if(k == 1) {
                        L bord ;
                        if(a[1] == 0) bord = m ;
                        else bord = 1 ;
                        if (bord != j) continue ;
                        dp[i][bord][1] = 1 ;
                        if (k == n){
                            ans += dp[i][j][k];
                            ans %= mod ;
                        }
                    } else {
                        dp[i][j][k] = sum(b[i], j, k);
                        if (k == n){
                            ans += dp[i][j][k];
                            ans %= mod ;
                        }
                    }
                }
            }
            for(L j = 1 ; j <= m ; j ++ ) {
                for(L k = 1 ; k <= n ; k ++ ) {
                    if(k == 1) {
                        L bord ;
                        if(a[1] == 0) bord = m ;
                        else bord = 1 ;
                        if (bord != j) continue ;
                        L x = min(bord,b[i]) ;
                        L y = max(bord,b[i]) ;
                        if(a[2] == 0) {
                            upda(x,y,1,y,2) ;
                        } else {
                            upda(x,y,1,x,2) ;
                        }
                    } else {
                        L x = min(j, b[i]);
                        L y = max(j, b[i]);
                        if (dp[i][j][k] == 0) continue;
                        if (a[k + 1] == 0) {
                            upda(x, y, dp[i][j][k], y, k + 1);
                        } else {
                            upda(x, y, dp[i][j][k], x, k + 1);
                        }
                    }
                }
            }
 
        }
        printf("%lld\n" , ans) ;
    }
}

J线段树会T。。。