[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/ ...
随机推荐
- DB.JDBC_jar_下载
1.Download Microsoft JDBC Driver for SQL Server - SQL Server _ Microsoft Docs.html(https://docs.micr ...
- OSPF与ACL 综合应用
1.企业内网运行OSPF路由协议,区域规划如图所示:2.财务和研发所在的区域不受其他区域链路不稳定性影响:3.R1.R2.R3只允许被IT登录管理:4.YF和CW之间不能互通,但都可以与IT互通:5. ...
- 【监控实践】【4.2】perfmon监控性能计数器（基于typeperf命令）
关键词:typeperf typeperf 命令 使用示例: 案例1:#使用typeperf收集windows cpu.内存.硬盘性能 #使用typeperf收集windows cpu.内存.硬盘性能 ...
- Python 入门之 Python三大器 之 装饰器
Python 入门之 Python三大器 之 装饰器 1.开放封闭原则: (1)代码扩展进行开放 任何一个程序,不可能在设计之初就已经想好了所有的功能并且未来不做任何更新和修改.所以我们必须允许代 ...
- 数学: HDU Co-prime
Co-prime Time Limit : 2000/1000ms (Java/Other) Memory Limit : 32768/32768K (Java/Other) Total Subm ...
- vue-split-table【表格合并和编辑插件】
前言 vue-split-table应用的效果图 vue-split-table开源地址,欢迎star,现在已经开源和同步到npm上轻松搞定表格拆分或者合并,编辑,再也不怕被产品怼啦 1.核心源码分析 ...
- XADC
XADC实验 1.XADC概述 Xilinx7系列内部自带一个双通道12位分辨率的高速(1MSPS 1M sample per second)采样速率的模拟混合信号处理模块,双通道的ADC支持单极和差 ...
- IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactoryClass(); 时,报COMException
解决办法: 在Program.cs的Main函数中添加如下代码: ESRI.ArcGIS.RuntimeManager.Bind(ESRI.ArcGIS.ProductCode.Desktop);
- MySQL之Foreign_Key
MySQL之Foregin_Key 一\\一对多 一.员工表和部门表 dep emp 类似与我们将所有的代码都写在一个py文件内 确立标语表之间的关系 思路:一定要要换位思考问题(必须两方都考虑周全之 ...
- MongoDB与python 交互
一.安装pymongo 注意 :当同时安装了python2和python3,为区分两者的pip,分别取名为pip2和pip3. 推荐:https://www.cnblogs.com/thunderLL ...