TwoSat

有N个布尔变量Xi,另外有M个需要满足的条件,每个条件形式为“Xi为真/假 则 Xj必为真/假”。是否能给所有的Xi赋值并且满足所有的条件。
struct TwoSat {
    vector<int> nxt[MAXN<<1];
    int mark[MAXN<<1];
    int stk[MAXN<<1], top;
    int n;
    
    void init(int tn) {
        n = tn;
        for(int i = 0; i < n; i++) {
            nxt[i<<1].clear();
            nxt[i<<1|1].clear();
            mark[i<<1] = mark[i<<1|1] = 0;
        }
    }
    
    void add(int x, int vx, int y, int vy) { //最多其中一个 
        int tx = (x<<1) + vx;
        int ty = (y<<1) + vy;
        nxt[tx].push_back(ty^1);
        nxt[ty].push_back(tx^1);
    }
    
    int dfs(int u) {
        if(mark[u^1]) return 0;
        if(mark[u]) return 1;
        mark[u] = 1;
        stk[top++] = u;
        int i, sz = nxt[u].size();
        for(i = 0; i < sz; i++) {
            if(!dfs(nxt[u][i])) return 0;
        }
        return 1;
    }
    
    int judge() {
        for(int i = 0; i < n; i++) {
            if(!mark[i<<1] && !mark[i<<1|1]) {
                top = 0;
                if(!dfs(i<<1)) {
                    while(top) mark[ stk[--top]] = 0;
                    if(!dfs(i<<1|1)) return 0;
                }
            }
        }
        return 1;
    }
}twosat;


#include <bits/stdc++.h> // #include "atcoder/convolution" // #include "atcoder/dsu" // #include "atcoder/fenwicktree" // #include "atcoder/lazysegtree" // #include "atcoder/math" // #include "atcoder/maxflow" // #include "atcoder/mincostflow" // #include "atcoder/modint" // #include "atcoder/scc" // #include "atcoder/segtree" // #include "atcoder/string" // #include "atcoder/twosat" // using namespace __gnu_cxx; // using namespace __gnu_pbds; using namespace std; auto rng = mt19937(random_device()()); auto rngl = mt19937_64(random_device()()); // Let's set a bit and flow! // I came, I divided, I conquered! int main() { ios_base::sync_with_stdio(false); cin.tie(0); cout.tie(0); int t; cin >> t; while (t --) { int n, q; cin >> n >> q; vector<long long> nums(n); for (auto &x: nums) cin >> x; vector<long long> vis = nums, cur = nums; vector<pair<int, int>> updates(q); for (auto &[idx, val]: updates) { cin >> idx >> val; idx --; cur[idx] += val; vis.emplace_back(cur[idx]); } sort(vis.begin(), vis.end()); vis.erase(unique(vis.begin(), vis.end()), vis.end()); int k = vis.size(); vector<int> fen(k + 1, 0); auto getidx = [&] (long long val) -> int { return lower_bound(vis.begin(), vis.end(), val) - vis.begin(); }; auto add = [&] (int idx, int val) -> void { idx ++; while (idx <= k) { fen[idx] += val; idx += idx & -idx; } }; auto not_exceed = [&] (int val) -> int { int res = 0, idx = 0; for (int i = 20; i >= 0; i --) { if ((idx | (1 << i)) <= k && res + fen[idx | (1 << i)] <= val) { idx |= 1 << i; res += fen[idx]; } } return res; }; for (
07-18
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值