[LOJ2288][THUWC2017]大葱的神力：搜索+背包DP+费用流+随机化
分析
测试点1、2：搜索+剪枝。
测试点3：只有一个抽屉，直接01背包。
测试点4、5：每个物品体积相同，说明每个抽屉能放下的物品个数固定，建图跑费用流。
测试点6：每个物品体积相近，经过验证发现每个抽屉能放下的物品个数仍然固定，费用流。
测试点7：除了第一个物品，其他物品体积相同。所以可以枚举第一个物品放在哪个抽屉，然后费用流。
测试点8、9、10：随机化贪心，如果不是人品爆发基本不可能过（博主没过）。
代码
测试点3
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cctype>
#include <algorithm>
#define rin(i,a,b) for(int i=(a);i<=(b);i++)
#define rec(i,a,b) for(int i=(a);i>=(b);i--)
#define trav(i,a) for(int i=head[(a)];i;i=e[i].nxt)
typedef long long LL;
using std::cin;
using std::cout;
using std::endl;
inline int read(){
int x;scanf("%d",&x);return x;
}
const int MAXN=2005;
const int MAXS=10005;
int n,m,s,v[MAXN],w[MAXN];
int f[MAXN][MAXS],pre[MAXN][MAXS];
int ans[MAXN];
int main(){
freopen("drawer3.in","r",stdin);
freopen("drawer3.out","w",stdout);
n=read(),m=read();
rin(i,1,n) v[i]=read();
s=read();
rin(i,1,n) w[i]=read();
rin(i,1,n) rec(j,s,v[i]){
f[i][j]=f[i-1][j];
pre[i][j]=j;
if(f[i-1][j-v[i]]+w[i]>f[i][j]){
f[i][j]=f[i-1][j-v[i]]+w[i];
pre[i][j]=j-v[i];
}
}
int now=s,pos=n;
while(pos){
if(pre[pos][now]<now) ans[pos]=1;
now=pre[pos][now];
pos--;
}
rin(i,1,n) printf("%d ",ans[i]);
printf("\n");
return 0;
}
测试点4、5、6
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cctype>
#include <algorithm>
#include <queue>
#define rin(i,a,b) for(int i=(a);i<=(b);i++)
#define rec(i,a,b) for(int i=(a);i>=(b);i--)
#define trav(i,a) for(int i=head[(a)];i;i=e[i].nxt)
typedef long long LL;
using std::cin;
using std::cout;
using std::endl;
inline int read(){
int x;scanf("%d",&x);return x;
}
const int MAXN=1005;
const int MAXM=205;
int n,m,S,T,mincost,v,c[MAXM],w[MAXN][MAXM],ecnt=1,head[MAXN+MAXM];
int dis[MAXN+MAXM],pre[MAXN+MAXM];
bool book[MAXN+MAXM];
std::queue<int> q;
struct Edge{
int to,nxt,cap,ost;
}e[(MAXN+MAXM+MAXN*MAXM)<<1];
inline void add_edge(int bg,int ed,int cap,int ost){
ecnt++;
e[ecnt].to=ed;
e[ecnt].nxt=head[bg];
e[ecnt].cap=cap;
e[ecnt].ost=ost;
head[bg]=ecnt;
}
inline bool spfa(){
memset(dis,0x3f,sizeof dis);
while(!q.empty()) q.pop();
dis[S]=0;
q.push(S);
book[S]=1;
while(!q.empty()){
int x=q.front();q.pop();
trav(i,x){
int ver=e[i].to;
if(dis[ver]>dis[x]+e[i].ost&&e[i].cap){
dis[ver]=dis[x]+e[i].ost;
pre[ver]=i;
if(!book[ver]){
q.push(ver);
book[ver]=1;
}
}
}
book[x]=0;
}
return dis[T]<1e9;
}
inline void upd(){
int now=T,flow=1e9;
while(now!=S){
flow=std::min(flow,e[pre[now]].cap);
now=e[pre[now]^1].to;
}
now=T;
while(now!=S){
e[pre[now]].cap-=flow;
e[pre[now]^1].cap+=flow;
now=e[pre[now]^1].to;
}
mincost+=flow*dis[T];
}
inline void ek(){
while(spfa()) upd();
}
int main(){
freopen("drawer6.in","r",stdin);
freopen("drawer6.out","w",stdout);
n=read(),m=read();
S=n+m+1,T=S+1;
rin(i,1,n){
v=read();
add_edge(S,i,1,0);
add_edge(i,S,0,0);
}
rin(i,1,m){
c[i]=read()/v;
add_edge(n+i,T,c[i],0);
add_edge(T,n+i,0,0);
}
rin(i,1,n) rin(j,1,m){
w[i][j]=read();
add_edge(i,n+j,1,-w[i][j]);
add_edge(n+j,i,0,w[i][j]);
}
ek();
rin(i,1,n){
bool flag=0;
trav(j,i){
int ver=e[j].to;
if(ver==S) continue;
if(!e[j].cap){
printf("%d ",ver-n);
flag=1;break;
}
}
if(!flag) printf("0 ");
}
printf("\n");
return 0;
}
测试点7
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cctype>
#include <algorithm>
#include <queue>
#define rin(i,a,b) for(int i=(a);i<=(b);i++)
#define rec(i,a,b) for(int i=(a);i>=(b);i--)
#define trav(i,a) for(int i=head[(a)];i;i=e[i].nxt)
typedef long long LL;
using std::cin;
using std::cout;
using std::endl;
inline int read(){
int x;scanf("%d",&x);return x;
}
const int MAXN=1005;
const int MAXM=205;
int n,m,S,T,mincost,v,c[MAXM],w[MAXN][MAXM],ecnt=1,head[MAXN+MAXM];
int dis[MAXN+MAXM],pre[MAXN+MAXM];
bool book[MAXN+MAXM];
std::queue<int> q;
struct Edge{
int to,nxt,cap,ost;
}e[(MAXN+MAXM+MAXN*MAXM)<<1];
inline void add_edge(int bg,int ed,int cap,int ost){
ecnt++;
e[ecnt].to=ed;
e[ecnt].nxt=head[bg];
e[ecnt].cap=cap;
e[ecnt].ost=ost;
head[bg]=ecnt;
}
inline bool spfa(){
memset(dis,0x3f,sizeof dis);
while(!q.empty()) q.pop();
dis[S]=0;
q.push(S);
book[S]=1;
while(!q.empty()){
int x=q.front();q.pop();
trav(i,x){
int ver=e[i].to;
if(dis[ver]>dis[x]+e[i].ost&&e[i].cap){
dis[ver]=dis[x]+e[i].ost;
pre[ver]=i;
if(!book[ver]){
q.push(ver);
book[ver]=1;
}
}
}
book[x]=0;
}
return dis[T]<1e9;
}
inline void upd(){
int now=T,flow=1e9;
while(now!=S){
flow=std::min(flow,e[pre[now]].cap);
now=e[pre[now]^1].to;
}
now=T;
while(now!=S){
e[pre[now]].cap-=flow;
e[pre[now]^1].cap+=flow;
now=e[pre[now]^1].to;
}
mincost+=flow*dis[T];
}
inline void ek(){
while(spfa()) upd();
}
int main(){
freopen("drawer7.in","r",stdin);
freopen("drawer7.out","w",stdout);
n=read(),m=read();
S=n+m+1,T=S+1;
int v1=read();
rin(i,2,n){
v=read();
}
rin(i,1,m){
c[i]=read();
}
rin(i,1,n) rin(j,1,m){
w[i][j]=read();
}
int minpos=-1,umincost=0;
rin(k,1,m){
ecnt=1,mincost=0;
memset(head,0,sizeof head);
c[k]-=v1;
mincost-=w[1][k];
rin(i,2,n){
add_edge(S,i,1,0);
add_edge(i,S,0,0);
}
rin(i,1,m){
add_edge(n+i,T,c[i]/v,0);
add_edge(T,n+i,0,0);
}
rin(i,2,n) rin(j,1,m){
add_edge(i,n+j,1,-w[i][j]);
add_edge(n+j,i,0,w[i][j]);
}
ek();
if(mincost<umincost){
umincost=mincost;
minpos=k;
}
c[k]+=v1;
}
ecnt=1,mincost=0;
memset(head,0,sizeof head);
rin(i,2,n){
add_edge(S,i,1,0);
add_edge(i,S,0,0);
}
rin(i,1,m){
add_edge(n+i,T,c[i]/v,0);
add_edge(T,n+i,0,0);
}
rin(i,2,n) rin(j,1,m){
add_edge(i,n+j,1,-w[i][j]);
add_edge(n+j,i,0,w[i][j]);
}
ek();
if(mincost<umincost){
umincost=mincost;
minpos=0;
}
else{
int k=minpos;
ecnt=1,mincost=0;
memset(head,0,sizeof head);
c[k]-=v1;
mincost-=w[1][k];
rin(i,2,n){
add_edge(S,i,1,0);
add_edge(i,S,0,0);
}
rin(i,1,m){
add_edge(n+i,T,c[i]/v,0);
add_edge(T,n+i,0,0);
}
rin(i,2,n) rin(j,1,m){
add_edge(i,n+j,1,-w[i][j]);
add_edge(n+j,i,0,w[i][j]);
}
ek();
}
printf("%d ",minpos);
rin(i,2,n){
bool flag=0;
trav(j,i){
int ver=e[j].to;
if(ver==S) continue;
if(!e[j].cap){
printf("%d ",ver-n);
flag=1;break;
}
}
if(!flag) printf("0 ");
}
printf("\n");
return 0;
}
测试点8、9、10
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cctype>
#include <algorithm>
#include <ctime>
#include <vector>
#define rin(i,a,b) for(int i=(a);i<=(b);i++)
#define rec(i,a,b) for(int i=(a);i>=(b);i--)
#define trav(i,a) for(int i=head[(a)];i;i=e[i].nxt)
typedef long long LL;
using std::cin;
using std::cout;
using std::endl;
inline int read(){
int x;scanf("%d",&x);return x;
}
const int MAXN=1005;
const int MAXM=305;
int n,m,v[MAXN],c[MAXM],blg[MAXN];
bool vis[MAXN];
std::vector<int> vec[MAXM],maxvec[MAXM];
struct Pair{
int id,w;
inline friend bool operator < (Pair A,Pair B){
return 1.0*A.w/v[A.id]>1.0*B.w/v[B.id];
}
}a[MAXM][MAXN];
inline int work(){
memset(vis,0,sizeof vis);
int ret=0;
rin(i,1,m){
vec[i].clear();
int rem=c[i];
rin(j,1,n){
if(vis[a[i][j].id]||rem<v[a[i][j].id]) continue;
if((rand()%(n*2))>n+vec[i].size()) continue;
vec[i].push_back(a[i][j].id);
vis[a[i][j].id]=1;
rem-=v[a[i][j].id];
ret+=a[i][j].w;
}
}
return ret;
}
int main(){
freopen("drawer10.in","r",stdin);
freopen("drawer10.out","w",stdout);
srand(time(0));
n=read(),m=read();
rin(i,1,n) v[i]=read();
rin(i,1,m) c[i]=read();
rin(i,1,n) rin(j,1,m) a[j][i]=(Pair){i,read()};
rin(i,1,m) std::sort(a[i]+1,a[i]+n+1);
int maxans=0;
rin(k,1,10000){
int temp=work();
if(temp>maxans){
memset(blg,0,sizeof blg);
rin(i,1,m) rin(j,0,(int)vec[i].size()-1)
blg[vec[i][j]]=i;
maxans=temp;
}
}
rin(i,1,n) printf("%d ",blg[i]);
printf("\n");
return 0;
}
[LOJ2288][THUWC2017]大葱的神力：搜索+背包DP+费用流+随机化的更多相关文章
- 配置魔药 [NOIP模拟] [DP] [费用流]
问题描述在<Harry Potter and the Chamber of Secrets>中,Ron 的魔杖因为坐他老爸的 Flying Car 撞到了打人柳,不幸被打断了,从此之后,他 ...
- AIM Tech Round 3 (Div. 1) (构造,树形dp,费用流,概率dp)
B. Recover the String 大意: 求构造01字符串使得子序列00,01,10,11的个数恰好为$a_{00},a_{01},a_{10},a_{11}$ 挺简单的构造, 注意到可以通 ...
- 背包dp整理
01背包 动态规划是一种高效的算法.在数学和计算机科学中,是一种将复杂问题的分成多个简单的小问题思想 ---- 分而治之.因此我们使用动态规划的时候,原问题必须是重叠的子问题.运用动态规划设计的算法比 ...
- LOJ 2288「THUWC 2017」大葱的神力
LOJ 2288「THUWC 2017」大葱的神力 Link Solution 比较水的提交答案题了吧 第一个点爆搜 第二个点爆搜+剪枝,我的剪枝就是先算出 \(mx[i]\) 表示选取第 \(i \ ...
- Codeforces 946 课程表背包DP 数位DFS构造
A B 给你A,B 两个数 1.a=0 OR b=0 break 2.a>=2b a=a-2b 3.b>=2a b=b-2a 如果只是单纯模拟肯定会超时 ...
- @loj - 2288@「THUWC 2017」大葱的神力
目录 @description@ @solution@ @data - 1@ @data - 2@ @data - 3@ @data - 4@ @data - 5@ @data - 6@ @data ...
- hdu 5534 Partial Tree 背包DP
Partial Tree Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid= ...
- HDU 5501 The Highest Mark 背包dp
The Highest Mark Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?p ...
- Codeforces Codeforces Round #319 (Div. 2) B. Modulo Sum 背包dp
B. Modulo Sum Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://codeforces.com/contest/577/problem/ ...
随机推荐
- centos7 VM虚拟机在主机关机重启后，无法ping通
解决办法 1.虚拟机的某些网络相关的服务没有启动,打开电脑的服务启动相关服务 2.打开虚拟机的虚拟网络设置,恢复默认设置即可 3.判定虚拟网卡的网关和centos的网关是否一致,如果不一致,改成一致, ...
- python 并发编程 多进程 JoinableQueue
JoinableQueue和Queue 使用一样 这就像是一个Queue对象,但队列允许项目的使用者通知生成者项目已经被成功处理.通知进程是使用共享的信号和条件变量来实现的. JoinableQueu ...
- linux 系统目录权限
小结: 目录的读.写.执行权限 1. 可读r:表示具有浏览目录 下面文件及子目录的权限 ls dir 1)如果没有x权限,则不能进到目录,既无法执行cd dir 2)如果没有x权限,ls列表时可以看到 ...
- [转帖]CGI与ISAPI的区别(转)
CGI与ISAPI的区别(转) 不知道原始网站是哪个 博客园里面也是转帖的 https://www.cnblogs.com/eret9616/p/8515095.html 不过我还是不了解CGI和IS ...
- debian上安装mysql server
1 将mysql添加到apt的repository中 第一步,下载mysql提供的ppa文件 wget https://dev.mysql.com/get/mysql-apt-config_0.8.1 ...
- SELECT 与 SET给标量赋值
本文转自 :https://blog.csdn.net/perddy/article/details/4033406 SQL Server推荐使用 SET 而不是 SELECT 对变量进行赋值.当表达 ...
- SQL如何通过当前日期获取上周一日期【转】
--当前时间 select getdate() --当前时间周的起始日期(以周一为例) ,) --上周起始: ,,)) --上上周起始: ,,)) --上上上周起始:s elect ,,))
- 使用git命令删除branch
使用git删除server上的一个branch注意事项: 1. 使用命令:git push origin –delete 分支名 (分支名称需要去掉origin,如果有), git branch ...
- selenium 教程
selenium 本身是一套web自动化测试工具,但其经常被用于爬虫,解决一些复杂爬虫的问题. selenium 用于爬虫时,相当于模拟人操作浏览器. 浏览器驱动 使用 selenium 需要先安装 ...
- makemap - 为sendmail创建数据库映像表
SYNOPSIS(总览) [-N ] [-d ] [-f ] [-o ] [-r ] [-s ] [-v ] maptype mapname DESCRIPTION(描述) 创建 sendmail(8 ...