博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
Codeforces 348C Subset Sums 分块思想
阅读量:4478 次
发布时间:2019-06-08

本文共 2480 字,大约阅读时间需要 8 分钟。

把所有集合分成大于sqrt(n)和小于sqrt(n)的集合, 处理出每两个集合有交集是多大之后就可以在sqrt(n)的复杂内完成一次操作。

#include
#define LL long long#define LD long double#define ull unsigned long long#define fi first#define se second#define mk make_pair#define PLL pair
#define PLI pair
#define PII pair
#define SZ(x) ((int)x.size())#define ALL(x) (x).begin(), (x).end()#define fio ios::sync_with_stdio(false); cin.tie(0);using namespace std;const int N = 1e5 + 7;const int inf = 0x3f3f3f3f;const LL INF = 0x3f3f3f3f3f3f3f3f;const int mod = 998244353;const double eps = 1e-8;const double PI = acos(-1);template
inline void add(T& a, S b) {a += b; if(a >= mod) a -= mod;}template
inline void sub(T& a, S b) {a -= b; if(a < 0) a += mod;}template
inline bool chkmax(T& a, S b) { return a < b ? a = b, true : false;}template
inline bool chkmin(T& a, S b) { return a > b ? a = b, true : false;}const int B = 350;int n, m, q;LL a[N];int cnt[307][N];int cnt2[307][307];LL sum[N], addval[N];int to[N];int fto[N];vector
S[N];vector
big, sml;vector
vc[N];int main() { scanf("%d%d%d", &n, &m, &q); for(int i = 1; i <= n; i++) scanf("%lld", &a[i]); for(int i = 1; i <= m; i++) { int num; scanf("%d", &num); S[i].resize(num); for(int j = 0; j < num; j++) scanf("%d", &S[i][j]); if(num >= B) to[i]= SZ(big), fto[SZ(big)] = i, big.push_back(i); else to[i] = SZ(sml), fto[SZ(sml)] = i, sml.push_back(i); } for(int i = 0; i < SZ(big); i++) { for(auto& x : S[big[i]]) { vc[x].push_back(i); sum[fto[i]] += a[x]; } } for(int i = 0; i < SZ(sml); i++) { for(auto& x : S[sml[i]]) { for(auto& y : vc[x]) { cnt[y][i]++; } } } for(int i = 0; i < SZ(big); i++) { for(auto& x : S[big[i]]) { for(auto& y : vc[x]) { if(i < y) cnt2[i][y]++, cnt2[y][i]++; } } } while(q--) { char op[3]; scanf("%s", op); if(op[0] == '+') { int k, x; scanf("%d%d", &k, &x); if(SZ(S[k]) >= B) { sum[k] += 1LL * SZ(S[k]) * x; addval[k] += x; for(int i = 0; i < SZ(big); i++) { if(to[k] == i) continue; sum[fto[i]] += 1LL * x * cnt2[to[k]][i]; } } else { for(auto& id : S[k]) a[id] += x; for(int i = 0; i < SZ(big); i++) sum[fto[i]] += 1LL * cnt[i][to[k]] * x; } } else { int k; scanf("%d", &k); if(SZ(S[k]) >= B) { printf("%lld\n", sum[k]); } else { LL ans = 0; for(auto& x : S[k]) ans += a[x]; for(int i = 0; i < SZ(big); i++) ans += addval[fto[i]] * cnt[i][to[k]]; printf("%lld\n", ans); } } } return 0;}/**/

 

转载于:https://www.cnblogs.com/CJLHY/p/10863174.html

你可能感兴趣的文章
ios 查看UIView的层次继承关系工具
查看>>
C++二维数组讲解、二维数组的声明和初始化
查看>>
纹理映射和混合
查看>>
PHP获取域名、IP地址的方法
查看>>
php验证复选框的小例子
查看>>
Sql Server 判断表或数据库是否存在
查看>>
计算机网络
查看>>
iOS-浅谈runtime运行时机制
查看>>
redis运维的一些知识点
查看>>
【poj2459】 Feed Accounting
查看>>
js,jQuery获取html5的data-*属性
查看>>
优秀个人免费私有云OwnCloud 8.0终于发布 - 亮眼新功能初探简介
查看>>
Java 之外,是 Scala 还是 Groovy?【转载】
查看>>
冒泡排序
查看>>
软工之404 Note Found团队
查看>>
关于事件冒泡和浏览器默认行为
查看>>
自动测试框架PyUnit
查看>>
SVN使用—概念及生命周期
查看>>
Solr与Lucene的区别
查看>>
js设计模式--方法的链式调用及回调
查看>>