专栏首页数据结构与算法loj#6073. 「2017 山东一轮集训 Day5」距离(费用流)

loj#6073. 「2017 山东一轮集训 Day5」距离(费用流)

题意

题目链接

Sol

我们可以把图行列拆开,同时对于行/列拆成很多个联通块,然后考虑每个点所在的行联通块/列联通块的贡献。

可以这样建边

从S向每个行联通块连联通块大小条边,每条边的容量为1,费用为\(i\)(i表示这是第几条边)。

从每个点所在的行联通块向列联通块连边,容量为1,费用为0

从每个列联通块向T连联通块大小条边,每条边的容量为1,费用为\(i\)(i表示这是第几条边)。

这样跑最小费用最大流,每增光一次的费用就是答案。预处理后O(1)回答即可

#include<bits/stdc++.h> 
#define Pair pair<int, int>
#define MP(x, y) make_pair(x, y)
#define fi first
#define se second
//#define int long long 
#define LL long long 
#define ull unsigned long long 
#define Fin(x) {freopen(#x".in","r",stdin);}
#define Fout(x) {freopen(#x".out","w",stdout);}
using namespace std;
const int MAXN = 5001, mod = 1e9 + 7, INF = 1e9 + 10;
const double eps = 1e-9;
template <typename A, typename B> inline bool chmin(A &a, B b){if(a > b) {a = b; return 1;} return 0;}
template <typename A, typename B> inline bool chmax(A &a, B b){if(a < b) {a = b; return 1;} return 0;}
template <typename A, typename B> inline LL add(A x, B y) {if(x + y < 0) return x + y + mod; return x + y >= mod ? x + y - mod : x + y;}
template <typename A, typename B> inline void add2(A &x, B y) {if(x + y < 0) x = x + y + mod; else x = (x + y >= mod ? x + y - mod : x + y);}
template <typename A, typename B> inline LL mul(A x, B y) {return 1ll * x * y % mod;}
template <typename A, typename B> inline void mul2(A &x, B y) {x = (1ll * x * y % mod + mod) % mod;}
template <typename A> inline void debug(A a){cout << a << '\n';}
template <typename A> inline LL sqr(A x){return 1ll * x * x;}
template <typename A, typename B> inline LL fp(A a, B p, int md = mod) {int b = 1;while(p) {if(p & 1) b = mul(b, a);a = mul(a, a); p >>= 1;}return b;}
template <typename A> A inv(A x) {return fp(x, mod - 2);}
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}

int N, S, T , TT;
char s[51][51];
int id[51][51][2], c1 = 1, c2 = 1, ans[MAXN * MAXN], tot1[20 * MAXN], tot2[20 * MAXN ], num1, num2;
struct Edge {
    int u, v, w, f, nxt;
}E[2 * MAXN * MAXN];
int head[MAXN  * 20 + 1], num;
void add_edge(int x, int y, int w, int f) {
    E[num] = (Edge){x, y, w, f, head[x]};
    head[x] = num++;
}
void AddEdge(int x, int y, int w, int f) {
    //printf("%d %d %d %d\n", x, y, w, f);
    add_edge(x, y, w, f);
    add_edge(y, x, -w, 0);
}
int dis[MAXN * 10], vis[MAXN * 10], pre[MAXN * 10];
int SPFA() {
    memset(dis, 0x3f, sizeof(dis));
    memset(vis, 0, sizeof(vis));
    queue<int> q; q.push(S); dis[S] = 0;
    while(!q.empty()) {
        int p = q.front(); q.pop(); vis[p] = 0;
        for(int i = head[p]; ~i; i = E[i].nxt) {
            int to = E[i].v, w = E[i].w;
            if(dis[to] > dis[p] + w && E[i].f) {
                dis[to] = dis[p] + w; pre[to] = i;
                if(!vis[to]) vis[to] = 1, q.push(to);
            }
        }
    }
    return dis[TT];
}
int MCMF() {
    int val = SPFA(), dec = INF;
    for(int k = TT; k != S; k = E[pre[k]].u) chmin(dec, E[pre[k]].f);
    for(int k = TT; k != S; k = E[pre[k]].u) E[pre[k]].f -= dec, E[pre[k] ^ 1].f += dec;
    return dec * val;
}
signed main() {
    //freopen("a.in", "r", stdin);
    memset(head, -1, sizeof(head));
    N = read(); S = 0; T = N * N * 10, TT = T + 1; c2 = N * N * 3 + 1;
    for(int i = 1; i <= N; i++) scanf("%s", s[i] + 1);
    for(int i = 1; i <= N; i++) {
        for(int j = 1; j <= N; j++) {
            if(s[i][j] == '#') tot1[c1] = num1, num1 = 0, c1++; 
            else id[i][j][0] = c1, num1++;
            if(s[j][i] == '#') tot2[c2] = num2, num2 = 0, c2++;
            else id[j][i][1] = c2, num2++;
        }
        if(num1) tot1[c1++] = num1, num1 = 0;
        if(num2) tot2[c2++] = num2, num2 = 0;
    }
    for(int i = 1; i <= N; i++)
        for(int j = 1; j <= N; j++)
            if(id[i][j][0] && id[i][j][1])
                AddEdge(id[i][j][0], id[i][j][1], 0, 1);
    for(int i = 1; i <= c1; i++) 
        for(int j = 0; j < tot1[i]; j++)
            AddEdge(S, i, j, 1);
    for(int i = N * N * 3 + 1; i <= c2; i++)
        for(int j = 0; j < tot2[i]; j++)
            AddEdge(i, T, j, 1);
    for(int i = 1; i <= 2 * N * N; i++)
        AddEdge(T, TT, 0, 1);
    for(int i = 1; i <= N * N; i++)
        ans[i] = ans[i - 1] + MCMF();
    int Q = read();
    while(Q--) cout << ans[read()] << '\n';
    return 0;
}

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • loj#6073. 「2017 山东一轮集训 Day5」距离(树链剖分 主席树)

    首先对询问差分一下,我们就只需要统计\(u, v, lca(u, v), fa[lca(u, v)]\)到根的路径的贡献。

    attack
  • 互联网早知道

    1、拼多多第三季度总营收33.72亿元 同比增长697% 股价大涨 距赶超京东市值仅差26亿美元 2、酷派向深圳法院提起诉讼,称锤子科技“欠钱不还” 3、有偿删...

    程序员的酒和故事
  • 深度|猎聘大数据​:北京人才供需及流动全景分析

    2017年5月2日,国内领先的中高端人才职业发展平台猎聘发布了《2017年一季度北京人才供需及流动全景报告》。本次报告主要基于猎聘大数据 ,结果显示,在2017...

    灯塔大数据
  • 3.4 VR扫描:Google取消2020 I/O开发者大会;Facebook Horizon开测Alpha版本

    (VRPinea 3月4日电)今日重点新闻:Google取消2020 I/O开发者大会现场活动,或探索线上直播;3月底前购买国行任天堂Switch的用户,将免费...

    VRPinea
  • 冲击IPO:达达的负“重”上市之路

    放在如今的美股市场而言,敢于面对惨淡的市场逆势冲击上市,也称得上真的勇士。前有疫情的影响,后有瑞幸自爆财务造假引发中概股被众多做空机构轮番做空,美股市场到如今更...

    刘旷
  • Few-shot Learning进展调研

    小样本学习主要研究如何通过少量样本学习识别模型。目前学术界普遍研究的是N-way-K-shot问题,即进行N个类别的识别,每类有K个样本。训练过程以task为单...

    公众号机器学习与生成对抗网络
  • 使用Python对图像进行中值滤波

    首先解答上一篇文章Python使用标准库subprocess调用外部程序中的问题,该题答案为['1', '2', '3', '4'],在正则表达式中,问号(?)...

    Python小屋屋主
  • 地铁数据下的真实生活:通勤缩短三分钟,房价贵出一辆车

    根据中原集团研究中心分析,每距离办公区域近一站地铁,或者说缩短三分钟上班路程,上海的一套房均价就要抬高18.75万元,相当于一辆普通轿车。在11月30日的“20...

    DT数据侠
  • 小米造车再加码!收购深动科技自动驾驶公司| 镁客网每周硬科技领域投融资汇总(7.3-7.9)

    本周硬科技领域投融资事件一共31起,人工智能领域发生12起融资事件,占比39%;半导体领域发生10起融资事件,占比33%;生物医药领域发生7起融资事件,占比22...

    镁客网
  • 选好域名很重要 3.35亿美元投资花落大搜车

    阿里巴巴领投,华平投资、春华资本、招银国际跟投完成对大搜车3.35亿美元(22.115亿元)E轮融资。这次融资距离上次2016年11月蚂蚁金服领投1亿美金投资不...

    躲在树上的域小名
  • 数据告诉你丨北京上班族的“三宗罪”

    北京很大,大得可以容下所有人的野心和梦想;北京很小,小得容不下一个人的孤独和泪水。在这里,繁华与残酷同在,挑战与机遇并存。每天,无数怀揣梦想的人向这里汇聚,每天...

    钱塘数据
  • 耶鲁大学博士讲述:中美精英生活方式大不同

    我的生活圈子里,成功者较多。这当然和高学历的教育背景有关系。最近渐渐发现,这些朋友中的中国人,虽然大多和我一样是博士,但居住在美国的和居住在中国的在业余习好上越...

    用户1756920
  • 1.22 VR扫描:SNOW中国完成5000万美元A轮融资;Greenlight报告:2017年AR/VR融资达19亿美元

    VRPinea
  • 29岁苏宁上市副总裁,跟随王健林万达黄金十年,他为何选择扎根于乡村?

    当BAT们“下沉”到三四线城市寻求增长空间时,旅游行业也在迎来一次前所未有的下沉潮流:乡村旅游正在中国兴起,并且成为国家乡村振兴战略过程中的重要抓手。 日前,...

    罗超频道
  • 塔荐 |【重磅】国家旅游局数据中心中国电信联合实验室发布《2017年上半年全国自驾游(跨市)报告》

    摘要:8月30日上午,由国家旅游局数据中心、中国电信集团公司、银联商务股份有限公司共同举办的“旅游大数据联合实验室签约揭牌仪式暨国内旅游大数据发布会”在国家旅游...

    灯塔大数据
  • 骑行大数据:北京最热门的购物圣地,居然不是西单

    畅骑1公里 骑行者们有更多消费选择 随着城市范围的扩大和多中心化,多类商圈鼎立的场景屡见不鲜。想在这个商圈吃网红餐厅,又想在隔壁商圈“血拼”,如何安排好一天的行...

    BestSDK
  • 新技术重塑全球供应链格局

    在智能制造和万物互联分论坛上,嘉宾们正在分享观点。 周明阳摄 新一代信息技术与制造业深度融合,正在引发一场影响深远的产业变革,形成新的生产方式、产业形态、商业模...

    企鹅号小编
  • 一周AI看点 | 苹果Swift之父离职特斯拉加盟Google Brain,全球首颗AI芯片麒麟970下单4000万颗

    本期一周AI看点包括:行业热点、技术应用、投融资、大咖观点以及技术前沿。 行业 苹果 Swift 之父 Chris Lattner 离职特斯拉加盟 Google...

    AI科技大本营
  • 【分析】比特大陆招股说明书解读[附pdf下载]

    2018年9月26日晚间消息,今天晚上,比特大陆在港交所递交A1招股文件,正式启动了香港联交所主板上市的计划,中金公司为独家保荐人。

    辉哥

扫码关注云+社区

领取腾讯云代金券