Bzoj的翻译出锅了所以来官方题面:

这个题应该是单向边而BZOJ说的是双向边,什么你WA了?谁叫你懒得看英文......

显然我们能正向反向两遍SPFA处理出每个点到总部的距离和总部到每个点的距离。
如果某个点所在的部门的大小为S,那么这个点需要送出S-1次消息并接收S-1次消息。
我们把每个点的两个距离求和并排序,显然在一个块中的是这个序列上的一个区间(脑补一下为什么不这样不优),我们做一下前缀和。
然后就开始DP了,f[i][j]表示前i个点分j个块,最小代价。f[i][j] = min( f[k][j-1] + ( i - k - 1 ) * ( sum[i] - sum[k] ) )。
这个DP是O(n^3)的,考虑优化。
显然更小的边权所在的块应该更大,所以我们能从区间[i-(i/j),i-1]枚举k,这样能优化到n^2logn。
然而有更优美的做法:显然随着j增大,对于每个i最优的k也是递增的,直接指针扫过去,复杂度O(n^2)。
两种做法都可以AC,反正我6代i5的机器上时间差异不大,不知道BZOJ能不能卡出来(我只交了第一种)。
(为什么我现在在刷这种水题?我也不知道啊!!!)

第一种代码:

 #include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
typedef long long int lli;
using namespace std;
const int maxn=5e3+1e2,maxe=5e4+1e2; lli su[maxn],f[][maxn];
int b,s,cur; struct Graph {
int s[maxn],t[maxe],nxt[maxe],l[maxe],dis[maxn],inq[maxn],cnt;
inline void addedge(int from,int to,int len) {
t[++cnt] = to , l[cnt] = len , nxt[cnt] = s[from] , s[from] = cnt;
}
inline void spfa(int st) {
memset(dis,0x3f,sizeof(dis)) , dis[st] = ;
queue<int> q; q.push(st) , inq[st] = ;
while( q.size() ) {
const int pos = q.front(); q.pop() , inq[pos] = ;
for(int at=s[pos];at;at=nxt[at])
if( dis[t[at]] > dis[pos] + l[at] ) {
dis[t[at]] = dis[pos] + l[at];
if( !inq[t[at]] ) q.push(t[at]);
}
}
}
}gra,rev; inline void dp() {
for(int i=;i<=b;i++) su[i] = (lli) gra.dis[i] + rev.dis[i];
sort(su+,su++b) , memset(f,0x3f,sizeof(f)) , **f = ;
for(int i=;i<=b;i++) su[i] += su[i-];
for(int j=;j<=s;j++) { // j is number of groups .
cur ^= , memset(f[cur],0x3f,sizeof(f[]));
for(int i=;i<=b;i++) // i is last node .
for(int lst=i/j;lst;lst--)
f[cur][i] = min( f[cur][i] , f[cur^][i-lst] + ( lst - ) * ( su[i] - su[i-lst] ) );
}
} int main() {
static int n,r;
scanf("%d%d%d%d",&n,&b,&s,&r);
for(int i=,a,b,l;i<=r;i++) scanf("%d%d%d",&a,&b,&l) , gra.addedge(a,b,l) , rev.addedge(b,a,l);
gra.spfa(b+) , rev.spfa(b+) , dp();
printf("%lld\n",f[cur][b]);
return ;
}

第二种代码:

 #include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
typedef long long int lli;
using namespace std;
const int maxn=5e3+1e2,maxe=5e4+1e2; int tp[maxn];
lli su[maxn],f[][maxn];
int b,s,cur; struct Graph {
int s[maxn],t[maxe],nxt[maxe],l[maxe],dis[maxn],inq[maxn],cnt;
inline void addedge(int from,int to,int len) {
t[++cnt] = to , l[cnt] = len , nxt[cnt] = s[from] , s[from] = cnt;
}
inline void spfa(int st) {
memset(dis,0x3f,sizeof(dis)) , dis[st] = ;
queue<int> q; q.push(st) , inq[st] = ;
while( q.size() ) {
const int pos = q.front(); q.pop() , inq[pos] = ;
for(int at=s[pos];at;at=nxt[at])
if( dis[t[at]] > dis[pos] + l[at] ) {
dis[t[at]] = dis[pos] + l[at];
if( !inq[t[at]] ) q.push(t[at]);
}
}
}
}gra,rev; inline void dp() {
for(int i=;i<=b;i++) su[i] = (lli) gra.dis[i] + rev.dis[i];
sort(su+,su++b) , memset(f,0x3f,sizeof(f)) , **f = ;
for(int i=;i<=b;i++) su[i] += su[i-];
for(int j=;j<=s;j++) { // j is number of groups .
cur ^= , memset(f[cur],0x3f,sizeof(f[]));
for(int i=;i<=b;i++) // i is last node .
for(int lst=tp[i];lst<i;lst++)
if( f[cur][i] >= f[cur^][lst] + ( i - lst - ) * ( su[i] - su[lst] ) ) f[cur][i] = f[cur^][lst] + ( i - lst - ) * ( su[i] - su[lst] ) , tp[i] = lst;
}
} int main() {
static int n,r;
scanf("%d%d%d%d",&n,&b,&s,&r);
for(int i=,a,b,l;i<=r;i++) scanf("%d%d%d",&a,&b,&l) , gra.addedge(a,b,l) , rev.addedge(b,a,l);
gra.spfa(b+) , rev.spfa(b+) , dp();
printf("%lld\n",f[cur][b]);
return ;
}

良心的我给的数据下载:
链接:https://pan.baidu.com/s/19EmgxmCYDASzpTaqr0alkA 密码:00qr

果てないこの闇の向こうにも
在无尽的这黑暗的对面
光があると信じてる
我相信也存在光明
ここから生まれ変わる世界だけ
从此只盯住这脱胎换骨的世界
见つめて离さないよ
目不转睛
広がるこの空を见上げると
仰望这广袤的天空
あの日の戦いが映る
那一日的决战映在脑海
いつかは全て消えてしまうのか
终有一天一切将会消失
栄光を取り戻せ
重新夺回这荣光

4609: [Wf2016]Branch Assignment 最短路 DP (阅读理解题)的更多相关文章

  1. Gym - 101611D Decoding of Varints(阅读理解题 )

    Decoding of Varints ​ 题意&思路: 首先根据红色边框部分的公式算出x,再有绿色部分得知,如果x是偶数则直接除以2,x是奇数则(x+1)/-2. PS:这题有数据会爆掉un ...

  2. 阅读关于DuReader:百度大规模的中文机器阅读理解数据集

    很久之前就得到了百度机器阅读理解关于数据集的这篇文章,今天才进行总结!.... 论文地址:https://arxiv.org/abs/1711.05073 自然语言处理是人工智能皇冠上的明珠,而机器阅 ...

  3. BZOJ_4609_[Wf2016]Branch Assignment_决策单调性+带权二分

    BZOJ_4609_[Wf2016]Branch Assignment_决策单调性+带权二分 Description 要完成一个由s个子项目组成的项目,给b(b>=s)个部门分配,从而把b个部门 ...

  4. hdu 4568 Hunter 最短路+dp

    Hunter Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Subm ...

  5. BZOJ_1003_[ZJOI2006]物流运输_最短路+dp

    BZOJ_1003_[ZJOI2006]物流运输_最短路+dp 题意:http://www.lydsy.com/JudgeOnline/problem.php?id=1003 分析: 这种一段一段的显 ...

  6. POJ 3635 Full Tank? 【分层图/最短路dp】

    任意门:http://poj.org/problem?id=3635 Full Tank? Time Limit: 1000MS   Memory Limit: 65536K Total Submis ...

  7. [USACO07NOV]牛继电器Cow Relays (最短路,DP)

    题目链接 Solution 非正解 似乎比较蛇啊,先个一个部分分做法,最短路+\(DP\). 在求最短路的堆或者队列中存储元素 \(dis_{i,j}\) 代表 \(i\) 这个节点,走了 \(j\) ...

  8. P1772 [ZJOI2006]物流运输 最短路+DP

    思路:最短路+DP 提交:1次 题解: $f[i]$表示到第$i$天的最小代价,我们可以预先处理出$i,j$两天之间(包括$i,j$)都可通行的最短路的代价记做$s[i][j]$,然后有$f[i]=m ...

  9. bzoj1003物流运输 最短路+DP

    bzoj1003物流运输 题目描述 物流公司要把一批货物从码头A运到码头B.由于货物量比较大,需要n天才能运完.货物运输过程中一般要转停好几个码头.物流公司通常会设计一条固定的运输路线,以便对整个运输 ...

随机推荐

  1. mysql系列十一、mysql优化笔记:表设计、sql优化、配置优化

    可以从这些方面进行优化: 数据库(表)设计合理 SQL语句优化 数据库配置优化 系统层.硬件层优化 数据库设计 关系数据库三范式 1NF:字段不可分; 2NF:有主键,非主键字段依赖主键; 3NF:非 ...

  2. Python3学习笔记06-字符串

    可以使用引号('或")来创建字符串. var1 = 'Hello World!' var2 = "Runoob" 在最新的Python 3版本中,字符串是以Unicode ...

  3. 【转】snprintf()函数使用方法

    众所周知,sprintf不能检查目标字符串的长度,可能造成众多安全问题,所以都会推荐使用snprintf. 注:sprintf()函数:int sprintf( char *buffer, const ...

  4. C:malloc/calloc/realloc/alloca内存分配函数

    原文地址:http://www.cnblogs.com/3me-linux/p/3962152.html calloc(), malloc(), realloc(), free(),alloca() ...

  5. windows系统中搭建Jenkins服务器

    1       须知 100.126.36.232等Jenkins服务器是通过设置代理访问外网,管理Jenkins和插件升级站点的,本地安装受黄区网络限制需要特殊配置,且有些插件无法下载. 前提条件: ...

  6. 转载:2.2.3 配置项的注释《深入理解Nginx》(陶辉)

    原文:https://book.2cto.com/201304/19628.html 如果有一个配置项暂时需要注释掉,那么可以加"#"注释掉这一行配置.例如: #pid       ...

  7. shell expect权威指南和实战

    一.概述 我们通过Shell可以实现简单的控制流功能,如:循环.判断等.但是对于需要交互的场合则必须通过人工来干预,有时候我们可能会需要实现和交互程序如telnet服务器等进行交互的功能.而expec ...

  8. Fragment的详细使用

    一直在用Fragment,但是没有系统的整理过,Google了一下相关文章,看到了几篇,将几篇还不错的文章重点整理了下,很多是直接Copy的,只为做个笔记,以后翻来看比较方便,建议大家看一下下面几篇, ...

  9. CGAffineTransform 缩放 / 旋转 / 平移

    CGAffineTransform此类是一个3*3矩阵的变换. - (void)transformImageView { CGAffineTransform t = CGAffineTransform ...

  10. 将日期或数据转换为char数据类型 TO_CHAR(x[[,c2],C3])

    TO_CHAR(x[[,c2],C3])[功能]将日期或数据转换为char数据类型[参数]x是一个date或number数据类型.c2为格式参数c3为NLS设置参数如果x为日期nlsparm=NLS_ ...