最小费用最大流可解最优解。至于dif如何解,可以把w扩大100倍,如果mission编号和排列P相等则对w+1,然后建立网络流。
对结果取模100可以得到没有改变mission的company数目,用company数目减之可以得到dif.

 /* 2853 */
#include <iostream>
#include <string>
#include <map>
#include <queue>
#include <set>
#include <stack>
#include <vector>
#include <deque>
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <cstring>
#include <climits>
#include <cctype>
#include <cassert>
#include <functional>
#include <iterator>
#include <iomanip>
using namespace std;
//#pragma comment(linker,"/STACK:102400000,1024000") #define sti set<int>
#define stpii set<pair<int, int> >
#define mpii map<int,int>
#define vi vector<int>
#define pii pair<int,int>
#define vpii vector<pair<int,int> >
#define rep(i, a, n) for (int i=a;i<n;++i)
#define per(i, a, n) for (int i=n-1;i>=a;--i)
#define clr clear
#define pb push_back
#define mp make_pair
#define fir first
#define sec second
#define all(x) (x).begin(),(x).end()
#define SZ(x) ((int)(x).size())
#define lson l, mid, rt<<1
#define rson mid+1, r, rt<<1|1 const int INF = 0x1f1f1f1f;
const int maxn = ;
const int maxv = maxn * ;
const int maxe = maxv * maxv * ;
int V[maxe], F[maxe], W[maxe], nxt[maxe];
int head[maxv], dis[maxv], pre[maxv], ID[maxv];
bool visit[maxv];
int M[maxn][maxn], P[maxn];
int s, t, m; void addEdge(int u, int v, int f, int w) {
V[m] = v;
F[m] = f;
W[m] = w;
nxt[m] = head[u];
head[u] = m++; V[m] = u;
F[m] = ;
W[m] = -w;
nxt[m] = head[v];
head[v] = m++;
} bool bfs() {
queue<int> Q;
int u, v, k; memset(dis, INF, sizeof(dis));
memset(visit, false, sizeof(visit));
Q.push(s);
dis[s] = ; while (!Q.empty()) {
u = Q.front();
Q.pop();
visit[u] = false;
for (k=head[u]; k!=-; k=nxt[k]) {
v = V[k];
if (F[k] && dis[v]>dis[u]+W[k]) {
dis[v] = dis[u] + W[k];
ID[v] = k;
pre[v] = u;
if (!visit[v]) {
visit[v] = true;
Q.push(v);
}
}
}
} return dis[t]==INF;
} int MCMF() {
int ret = , tmp;
int u, v, k; while () {
if (bfs())
break; tmp = INF;
for (v=t, u=pre[v]; v!=s; v=u, u=pre[v]) {
k = ID[v];
tmp = min(F[k], tmp);
} for (v=t, u=pre[v]; v!=s; v=u, u=pre[v]) {
k = ID[v];
F[k] -= tmp;
F[k^] += tmp;
} ret += dis[t] * tmp;
} return ret;
} int main() {
ios::sync_with_stdio(false);
#ifndef ONLINE_JUDGE
freopen("data.in", "r", stdin);
freopen("data.out", "w", stdout);
#endif int r, c;
int ans, tot, tmp, dif; while (scanf("%d %d", &r, &c) != EOF) { s = m = ;
t = r + c + ;
memset(head, -, sizeof(head)); rep(i, , r+)
rep(j, , c+)
scanf("%d", &M[i][j]);
rep(i, , r+)
scanf("%d", &P[i]); rep(i, , r+)
addEdge(s, i, , ); rep(j, , c+)
addEdge(j+r, t, , ); tot = ;
rep(i, , r+) {
rep(j, , c+) {
if (j == P[i]) {
tmp = M[i][j] * + ;
} else {
tmp = M[i][j] * ;
}
addEdge(i, j+r, , -tmp);
}
tot += M[i][P[i]];
} tmp = -MCMF();
ans = tmp/ - tot;
#ifndef ONLINE_JUDGE
printf("tot = %d, MCMF = %d\n", tot, tmp);
#endif
dif = r - tmp%;
printf("%d %d\n", dif, ans);
} #ifndef ONLINE_JUDGE
printf("time = %d.\n", (int)clock());
#endif return ;
}

最新文章

  1. 第五篇 基于.net搭建热插拔式web框架(拦截器---请求管道)
  2. (一)安卓小app开发之基础环境搭建
  3. 修改更新源sources.list,提高软件下载安装速度(提供Kali 2.0 更新源)
  4. 启动hadoop报192.168.1.151: Address 192.168.1.151 maps to node1, but this does not map back to the address - POSSIBLE BREAK-IN ATTEMPT!
  5. ubuntu 折腾之路
  6. 15_会话技术_Cookie
  7. Delphi 客户端调用Webservice 的TClientdataset 报出“http://www.borland.com/namespaces/Types-IAppServerSOAP”
  8. Python sql数据的增删改查简单操作
  9. .NET中公共变量与属性的区别
  10. Netty(二)——TCP粘包/拆包
  11. MAVEN自动发布更新本地和远程仓库
  12. 知识点:Mysql 索引优化实战(3)
  13. Promise及Async/Await
  14. Xtion pro live OpenNI2.2 Nite 2.2 安装配置1.0
  15. 使ipconfig命令结果更整洁
  16. 深度学习demo
  17. 服务注册和发现(Consul)
  18. Word在转PDF的过程中如何创建标签快速方便阅读(图文详解)
  19. xargs 原理&amp;使用
  20. 初探Angular_03 组件中模板数据绑定

热门文章

  1. Android Activity的生命周期详解
  2. Angularjs中使用$location获取url参数时,遇到的坑~~~
  3. 记:mysql 连接超时解决办法
  4. 销毁session
  5. 04_过滤器Filter_01_入门简述
  6. 一次ora-1113 记录
  7. imagecreatefromjpeg(): gd-jpeg, libjpeg: recoverable error: Corrupt JPEG data: 1 extraneous bytes be
  8. void void*
  9. ci 多个文件同时上传
  10. 已经安装php后,再增加扩展模块(不重新编辑php)