[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/ ...
随机推荐
- 20191128 Spring Boot官方文档学习（9.4-9.8）
9.4.Spring MVC Spring Boot有许多启动器包含Spring MVC.请注意,一些启动器包括对Spring MVC的依赖,而不是直接包含它. 9.4.1.编写JSON REST服务 ...
- JS中同步与异步
不讲过多定义,举两个例子说明下 例一: console.log(100); setTimeout(function(){ console.log(200); },1000); console.log( ...
- Spark出现java.lang.stackoverflowerror的解决方法
正在测试的程序需要多次迭代(400+次),每次迭代有复杂的运算 迭代到100多次的时候报java.lang.stackoverflowerror的错误 解决方法:先checkpoint()再count ...
- Paper Reading_Computer Architecture
最近(以及预感接下来的一年)会读很多很多的paper......不如开个帖子记录一下读paper心得 Computer Architecture Last level cache (llc) perf ...
- c# 杀死占用某个文件的进程
原文:c# 杀死占用某个文件的进程 需要使用微软提供的工具Handle.exe string fileName = @"H:\abc.dll";//要检查被那个进程占用的文件 Pr ...
- 小程序之如和使用view内部组件来进行页面的排版功能
这篇文章主要介绍了关于小程序之如和使用view内部组件来进行页面的排版功能,有着一定的参考价值,现在分享给大家,有需要的朋友可以参考一下 涉及知识点: 1.垂直排列,水平排列 2.居中对齐 示例: 1 ...
- 03-Django-models
# Models 模型 - ORM - ObjectRelationMap : 把面向对象思想转换成关系数据库思想.操作上把类等价于表格 - 类对应表格 - 类中的属性对应表中的字段 - 在应用中的m ...
- 梯度下降算法(Gradient descent)GD
1.我们之前已经定义了代价函数J,可以将代价函数J最小化的方法,梯度下降是最常用的算法,它不仅仅用在线性回归上,还被应用在机器学习的众多领域中,在后续的课程中,我们将使用梯度下降算法最小化其他函数,而 ...
- TCP/IP四层协议
1.数据链路层 数据链路层实现了网卡接口的网络驱动程序,处理数据在物理媒介(以太网,令牌环)上的传输,常用协议包含ARP(地址解析协议),RARP(逆地址解析协议)两个协议,他们实现了IP地址和物理地 ...
- Vue下简单分页及搜索功能
最近利用Vue和element ui仿写了个小页面,记一哈分页和搜索功能的简单实现. 首页 emmmm..... 搜索框输入..... 搜索完成 数据是直接写在这里面的: cardPhoto:[ ...