【BZOJ 3661】 Hungry Rabbit （贪心、优先队列）
3661: Hungry Rabbit
Time Limit: 100 Sec Memory Limit: 512 MBSec Special Judge
Submit: 67 Solved: 47Description
可怕的洪水在夏天不期而至，兔子王国遭遇了前所未有的饥荒，它们不得不去外面的森林里寻找食物。
为了简化起见，我们假设兔子王国中有n只兔子，编号为1n。在救济粮到来之前的m天中，每天恰好有k只兔子需要去森林里寻找粮食。森林里居住着可怕的大灰狼，所幸兔子已经摸清了大灰狼捕食习惯，即狼们在每一天只会捕食特定编号的兔子。为了安全起见，兔子们需要保证每次出去觅食的k只兔子都不会被狼捕食。由于每天出去捕食的兔子都不尽相同，它们为每一天定义了一个生疏度pi，即第i天出来寻找食物，但是第i1天却没有出来觅食的兔子个数。规定第1天的生疏度为0.现在兔子们希望在保证安全的前提下，每天的生疏度不能超过L，请为兔子们构造一个合法的方案。Input
第一行包括四个整数n,m,k和L.
接下来n行，每行一个长度为m的01串。其中第i行第j个字符若为0，则表示狼在第j天会捕食编号为i的兔子，为1则表示不捕食。Output
m行，每行k个1-n之间互不相同的整数，代表这一天出去寻找食物的兔子编号。如果没有合法方案，则输出一行1即可。
Sample Input
5 4 3 1
1001
1101
1111
1110
0111
Sample Output
2 3 4
2 3 4
3 4 5
2 3 5
HINT
对于 100% 的测试数据，1 <= n;m <= 800; 1 <= k <= n; 1 <= l <= k
Source
【分析】
唉、、、我打了网络流，70分。。。【其实数据弱恰好我数组开小本来90分的。。
说说我的网络流打法吧，毕竟想了很久。
就是1100011的兔子可以看成两只1100000、0000011的兔子。
然后就是说每个兔子表示一个区间，表示它可以从l~r这几天工作。
就是要找若干个区间覆盖全区间k次。
就类似这样子的建图：
拆点那里的点容限制了<=L次交换。反向INF表示区间相交求并集。然后判断是否满流就是否有解。
然后有意义的表示区间的边就是答案，输出即可。【好像输出答案有点困难啊？？
然后正解是贪心！【表示考场上想过不知道为什么脑抽觉得n^2log过不了？？？【黑人问号？？？
就是每次贪心找尽量长的区间。。。
有优先队列维护，每次用最优的更新最差的看看是否成立即可。
1 #include<cstdio> 2 #include<cstdlib> 3 #include<cstring> 4 #include<iostream> 5 #include<algorithm> 6 #include<queue> 7 using namespace std; 8 #define Maxn 810 9 #define INF 0xfffffff 10 11 struct node 12 { 13 int l,r,id; 14 }t[Maxn*Maxn]; 15 struct cmp{bool operator()(node x,node y){return x.r<y.r;}}; 16 struct cmp2{bool operator()(node x,node y){return x.r>y.r;}}; 17 // bool cmp(node x,node y) {return x.r>y.r;} 18 // bool cmp2(node x,node y) {return x.r<y.r;} 19 bool cmp3(node x,node y) {return x.l<y.l;} 20 priority_queue<node,vector<node>,cmp> q1; 21 priority_queue<node,vector<node>,cmp2> q2; 22 23 char s[Maxn]; 24 25 bool vis[Maxn]; 26 int op[Maxn][Maxn]; 27 int main() 28 { 29 // int T; 30 // scanf("%d",&T); 31 // while(T--) 32 { 33 int n,m,k,l; 34 scanf("%d%d%d%d",&n,&m,&k,&l); 35 36 int cnt=0,id; 37 for(int i=1;i<=n;i++) 38 { 39 scanf("%s",s+1); 40 for(int j=1;j<=m;j++) 41 { 42 if((j==1||s[j-1]=='0')&&s[j]=='1') id=j; 43 else if(j!=1&&s[j-1]=='1'&&s[j]=='0') t[++cnt].l=id,t[cnt].r=j-1,t[cnt].id=i; 44 if(s[j]=='1'&&j==m) t[++cnt].l=id,t[cnt].r=m,t[cnt].id=i; 45 } 46 } 47 sort(t+1,t+1+cnt,cmp3); 48 // int qz=0; 49 while(!q1.empty()) q1.pop(); 50 while(!q2.empty()) q2.pop(); 51 int nw=0; 52 bool ok=1; 53 memset(vis,0,sizeof(vis)); 54 for(int i=1;i<=m;i++) 55 { 56 int hh=0; 57 while(t[nw+1].l<=i&&nw<cnt) q1.push(t[++nw]); 58 while(!q2.empty()&&q2.top().r<i) {vis[q2.top().id]=0;q2.pop();} 59 while(q2.size()<k) 60 { 61 if(q1.empty()||q1.top().r<i) {ok=0;break;} 62 vis[q1.top().id]=1; 63 q2.push(q1.top());q1.pop(); 64 hh++; 65 if(hh>l&&i!=1) break; 66 } 67 if((hh>l&&i!=1)||!ok) {ok=0;break;} 68 while(hh<l&&!q1.empty()&&q1.top().r>q2.top().r) 69 { 70 hh++; 71 vis[q2.top().id]=0; 72 q2.pop(); 73 vis[q1.top().id]=1; 74 q2.push(q1.top());q1.pop(); 75 } 76 op[i][0]=0; 77 for(int j=1;j<=n;j++) if(vis[j]) op[i][++op[i][0]]=j; 78 } 79 if(!ok) printf("1\n"); 80 else 81 { 82 for(int i=1;i<=m;i++) 83 { 84 for(int j=1;j<=op[i][0];j++) printf("%d ",op[i][j]); 85 printf("\n"); 86 } 87 } 88 } 89 return 0; 90 }
2017-04-24 20:03:42
