E - Product Development
题意
有一个产品有
K
K
K 个参数,初始都为
0
0
0,公司计划把所有参数都调到
≥
P
\geq P
≥P
有
N
N
N 个计划,每个计划都会对第
j
j
j 个参数有
A
i
,
j
A_{i,j}
Ai,j 的增加,并且这个计划花费是
C
i
C_i
Ci
求出满足条件的最小花费
思路
注意到 K K K 最大是 5 5 5,所以每个参数最多就是 0 − 5 , 6 0-5 , 6 0−5,6 个状态,那么 5 5 5 个参数就是 6 5 = 7776 6^5 = 7776 65=7776 种状态 。我们可以利用背包 和 状态压缩,从前面已经有的状态,推出新的状态。
可以使用一个 5 5 5 位的十进制数来表示状态,然后每一位模 10 10 10 就是当前参数的状态了。至于所有 ≥ P \geq P ≥P 的参数,全部表示成 P P P 就可以了。
时间复杂度: O ( N K ( P + 1 ) 5 ) O(NK(P+1)^5) O(NK(P+1)5)
// Problem: E - Product Development
// Contest: AtCoder - AtCoder Beginner Contest 322
// URL: https://atcoder.jp/contests/abc322/tasks/abc322_e
// Memory Limit: 1024 MB
// Time Limit: 4000 ms
//
// Powered by CP Editor (https://cpeditor.org)
#include<bits/stdc++.h>
#define fore(i,l,r) for(int i=(int)(l);i<(int)(r);++i)
#define fi first
#define se second
#define endl '\n'
const int INF=0x3f3f3f3f;
const long long INFLL=0x3f3f3f3f3f3f3f3fLL;
typedef long long ll;
const int N=105;
const int K=6;
ll c[N];
int A[N][K];
ll dp[N][100005];
int main(){
std::ios::sync_with_stdio(false);
std::cin.tie(nullptr);
std::cout.tie(nullptr);
int n,k,p;
std::cin>>n>>k>>p;
fore(i,1,n+1){
std::cin>>c[i];
fore(j,1,k+1) std::cin>>A[i][j];
}
fore(i,0,n+1)
fore(j,0,100005)
dp[i][j] = INFLL;
std::set<int> mp;
mp.insert(0);
dp[0][0] = 0;
fore(i,1,n+1){
fore(j,0,100005) dp[i][j] = dp[i-1][j]; //这里可以压缩成一位dp数组
std::vector<int> temp;
for(auto x : mp){
std::vector<int> v(k+1,0); //k parameters
int idx = 0;
int tmp = x;
while(x){ //old status
v[idx++] = x%10;
x/=10;
}
x = tmp;
fore(j,1,k+1){ //new status
v[j-1] += A[i][j];
if(v[j-1] > p) v[j-1] = p;
}
int status = 0;
fore(j,0,k) status += v[j] * pow(10,j);
dp[i][status] = std::min(dp[i][status],dp[i-1][x] + c[i]);
temp.push_back(status); //记录一下新增加的状态
}
for(auto x : temp) mp.insert(x); //插入新状态
}
int res = 0;
fore(j,0,k) res += p * pow(10,j);
std::cout<<(dp[n][res] == INFLL?-1:dp[n][res]);
return 0;
}
F - Vacation Query
题意
给定一个长度为 N N N 的 01 01 01 串,有 q q q 次操作,每次操作的形式是: ( c , L , R ) (c,L,R) (c,L,R)
- 如果 c = 1 c = 1 c=1:将 [ L , R ] [L,R] [L,R] 区间的 S i S_i Si 所有都取反
- 如果 c = 2 c = 2 c=2:输出 [ L , R ] [L,R] [L,R] 区间的最长连续 1 1 1 的长度
思路
由于涉及到 区间修改 和 区间查询,考虑使用 懒标记线段树 来实现
首先只考虑查询,不难发现:区间
[
L
,
R
]
[L,R]
[L,R] 的最长连续
1
1
1,一定来自三种途径:
- 来自区间 [ L , m i d ] [L,mid] [L,mid]
- 来自区间 [ m i d + 1 , R ] [mid+1,R] [mid+1,R]
- 来自左儿子的最长后缀 和 右儿子的最长前缀 合并起来,也就是贯穿中间的位置
因此我们就需要维护一个区间的 : 前缀 1 1 1 和 后缀 1 1 1。
再考虑修改,可以发现其实区间取反,就等价于把 前缀 1 1 1 和 前缀 0 0 0 交换一下,后缀也是一样的。因此我们还要维护一个区间的:前缀 0 0 0 和 后缀 0 0 0。
时间复杂度: O ( N l o g N ) O(NlogN) O(NlogN)
// Problem: F - Vacation Query
// Contest: AtCoder - AtCoder Beginner Contest 322
// URL: https://atcoder.jp/contests/abc322/tasks/abc322_f
// Memory Limit: 1024 MB
// Time Limit: 2000 ms
//
// Powered by CP Editor (https://cpeditor.org)
#include<bits/stdc++.h>
#define int long long
#define fore(i,l,r) for(int i=(int)(l);i<(int)(r);++i)
#define fi first
#define se second
#define endl '\n'
const int INF=0x3f3f3f3f;
const long long INFLL=0x3f3f3f3f3f3f3f3fLL;
typedef long long ll;
const int N=500005;
struct node{
int maxlen1;
int maxlen0;
int tag;
int pref0;
int pref1;
int suff0;
int suff1;
int l;
int r;
}tree[N<<2];
std::string s;
void build(int p,int l,int r){
tree[p] = {0,0,0,0,0,0,0,l,r};
if(l==r){
tree[p].maxlen1 = (s[l] == '1');
tree[p].maxlen0 = (s[l] == '0');
tree[p].pref0 = tree[p].suff0 = (s[l] == '0');
tree[p].pref1 = tree[p].suff1 = (s[l] == '1');
return;
}
int mid = l + r >> 1;
build(p<<1,l,mid);
build(p<<1|1,mid+1,r);
tree[p].pref0 = tree[p<<1].pref0;
if(mid - l + 1 == tree[p<<1].pref0) tree[p].pref0 += tree[p<<1|1].pref0;
tree[p].suff0 = tree[p<<1|1].suff0;
if(r - mid == tree[p<<1|1].suff0) tree[p].suff0 += tree[p<<1].suff0;
tree[p].pref1 = tree[p<<1].pref1;
if(mid - l + 1 == tree[p<<1].pref1) tree[p].pref1 += tree[p<<1|1].pref1;
tree[p].suff1 = tree[p<<1|1].suff1;
if(r - mid == tree[p<<1|1].suff1) tree[p].suff1 += tree[p<<1].suff1;
tree[p].maxlen1 = std::max({tree[p<<1].maxlen1,tree[p<<1|1].maxlen1,tree[p<<1].suff1 + tree[p<<1|1].pref1});
tree[p].maxlen0 = std::max({tree[p<<1].maxlen0,tree[p<<1|1].maxlen0,tree[p<<1].suff0 + tree[p<<1|1].pref0});
}
inline void push_down(int p){
tree[p].tag = 0;
tree[p<<1].tag ^= 1;
tree[p<<1|1].tag ^= 1;
std::swap(tree[p<<1].pref0,tree[p<<1].pref1);
std::swap(tree[p<<1].suff0,tree[p<<1].suff1);
std::swap(tree[p<<1].maxlen0,tree[p<<1].maxlen1);
std::swap(tree[p<<1|1].pref0,tree[p<<1|1].pref1);
std::swap(tree[p<<1|1].suff0,tree[p<<1|1].suff1);
std::swap(tree[p<<1|1].maxlen0,tree[p<<1|1].maxlen1);
}
void update(int p,int l,int r,int L,int R){
if(r < L || R < l) return;
if(l <= L && R <= r){
tree[p].tag ^= 1;
std::swap(tree[p].pref0,tree[p].pref1);
std::swap(tree[p].suff0,tree[p].suff1);
std::swap(tree[p].maxlen0,tree[p].maxlen1);
return;
}
int mid = L + R >> 1;
if(tree[p].tag) push_down(p);
if(l <= mid) update(p<<1,l,r,L,mid);
if(r > mid) update(p<<1|1,l,r,mid+1,R);
/* push_up */
tree[p].pref0 = tree[p<<1].pref0;
if(mid - L + 1 == tree[p<<1].pref0) tree[p].pref0 += tree[p<<1|1].pref0;
tree[p].suff0 = tree[p<<1|1].suff0;
if(R - mid == tree[p<<1|1].suff0) tree[p].suff0 += tree[p<<1].suff0;
tree[p].pref1 = tree[p<<1].pref1;
if(mid - L + 1 == tree[p<<1].pref1) tree[p].pref1 += tree[p<<1|1].pref1;
tree[p].suff1 = tree[p<<1|1].suff1;
if(R - mid == tree[p<<1|1].suff1) tree[p].suff1 += tree[p<<1].suff1;
tree[p].maxlen1 = std::max({tree[p<<1].maxlen1,tree[p<<1|1].maxlen1,tree[p<<1].suff1 + tree[p<<1|1].pref1});
tree[p].maxlen0 = std::max({tree[p<<1].maxlen0,tree[p<<1|1].maxlen0,tree[p<<1].suff0 + tree[p<<1|1].pref0});
}
int query(int p,int l,int r,int L,int R){
if(r < L || R < l) return 0;
if(l <= L && R <= r) return tree[p].maxlen1;
int mid = L + R >> 1;
if(tree[p].tag) push_down(p);
int res = 0;
if(l <= mid) res = query(p<<1,l,r,L,mid);
if(r > mid) res = std::max(res,query(p<<1|1,l,r,mid+1,R));
/* 可能 r 并没有包含右儿子的全部前缀1 */
if(l <= mid && r > mid) res = std::max(res,std::min(tree[p<<1].suff1,mid-l+1) + std::min(tree[p<<1|1].pref1,r-mid));
return res;
}
signed main(){
std::ios::sync_with_stdio(false);
std::cin.tie(nullptr);
std::cout.tie(nullptr);
int n,q;
std::cin>>n>>q;
std::cin>>s;
s = '0' + s;
build(1,1,n);
int opt,l,r;
while(q--){
std::cin>>opt>>l>>r;
if(opt==1) update(1,l,r,1,n);
else std::cout<<query(1,l,r,1,n)<<endl;
}
return 0;
}