二分答案全攻略

当直接求解答案比较困难时，可以枚举答案，再判断枚举得到的这个答案是否合法。
如果这个答案还具有单调性，就可以用二分加速枚举过程。

Kth Largest(http://acm.hrbeu.edu.cn/index.php?act=problem&id=1005&cid=14)(二次二分)
题意:
 A,B两个数组，求出两个数组互相乘以后生成的n*n个数中的第k大的数。

解法:
  对A，B从小到大排序。
  二分第k大的数的值ans。
  判定过程:sum置0，枚举指向A的i,二分B中的j，找到A[i]*B[j]>=ans的最小的j，将n-j+1累加到sum中，sum==k这ans合法否则不合法。
  最后算出ans再进行微调。

代码:

 #include <stdio.h> #include <stdlib.h> int A[10000+5],B[10000+5]; int n,k; int cmp(const void *a,const void *b) { return (*((int *)a))-(*((int *)b)); } int find(int total,int anther) { int left,right; int mid; left=0;right=n-1; while(left<right) { mid=(left+right)/2; if(anther*B[mid]<total) left=mid+1; else right=mid; } if(anther*B[left]<total&&left==n-1) return n; return left; } int getx() { int i,j; int right,left; int mid; int sum; left=A[0]*B[0];right=A[n-1]*B[n-1]; while(left<right) { mid=(right+left+1)/2; sum=0; for(i=0;i<n;i++) sum+=n-find(mid,A[i]); if(sum>=k) left=mid; else right=mid-1; } return left; } int min(int a,int b) { return a<b?a:b; } int main() { int cases; int i,j; int x; int ans; int temp; scanf("%d",&cases); while(cases--) { scanf("%d%d",&n,&k); for(i=0;i<n;i++) scanf("%d",&A[i]); for(i=0;i<n;i++) scanf("%d",&B[i]); qsort(A,n,sizeof(int),cmp); qsort(B,n,sizeof(int),cmp); x=getx(); ans=A[n-1]*B[n-1]; for(i=0;i<n;i++) { temp=find(x,A[i]); if(temp!=n) { ans=min(ans,B[temp]*A[i]); } } printf("%d/n",ans); } return 0; }
A Task Proces（http://acm.hdu.edu.cn/showproblem.php?pid=3433）（二分加DP）
题意:
  有x个任务A,y个任务B，有n个人做A的时间和做B的时间不一样。问最少的做完所有x个A任务和y个B任务的时间。

解法:
  二分做完的时间，对每个人计算：在该时间内，他在做了u个A任务的情况下还能做多少个B任务，记为f[i][u].然后进行一下DP，DP[i][u]表示前i个人在时间限制内做了u个任务的同时还能做多少个B

任务。DP[i][u]=max{dp[i-1][u-x]+F[i][x]},最后判断DP[n][x]>=y是否成立即可。

代码:

#include <stdio.h> #include <algorithm> #include <memory.h> #include <iostream> using namespace std; int f[50+5][200+5]; int a[200+5],b[200+5]; int dp[50][200+5]; int n,x,y; int max(int x,int y){return x>y?x:y;} bool go(int t) { int i,j,k; for(i=0;i<=n;i++) for(j=0;j<=x;j++) { f[i][j]=-10000000; dp[i][j]=-10000000; } for(i=1;i<=n;i++) { for(j=0;j<=x&&j*a[i]<=t;j++) { f[i][j]=(t-j*a[i])/b[i]; } } dp[0][0]=0; for(i=1;i<=n;i++) { for(j=0;j<=x;j++) { for(k=0;k<=j;k++) { //if((j-k)*b[i]>=t)continue; dp[i][j]=max(dp[i-1][k]+f[i][j-k],dp[i][j]); } } } if(dp[n][x]>=y)return true; else return false; } int main() { int t,cas=0,i; scanf("%d",&t); while(t--) { scanf("%d %d %d",&n,&x,&y); for(i=1;i<=n;i++) { scanf("%d %d",&a[i],&b[i]); } int left=0,right=400000; while(left<right) { int mid=(left+right)/2; if(go(mid)) { right=mid; } else left=mid+1; } printf("Case %d: %d/n",++cas,left); } return 0; }
organ（http://acm.nbu.edu.cn/JudgeOnline/showproblem?problem_id=1820）(二分加贪心)
题意:
 将一棵树分成m份,使每份节点的最小权值最小,求这个最小值。
解法:
先二分答案转换成判定性问题，判定的时候可以用贪心，从叶子节点往根节点枚举，如果当前子树能力和大于等于二分的那个平均值，就直接把这个子树独立为一组。到最后如果独立出来的组数大于等

于m，即当前二分的解是一个可行的解。
代码:
#include<iostream> using namespace std; #define MAXN 10005 int n,m,i,j; int a[MAXN],father[MAXN],as[MAXN],ras[MAXN]; bool calc(int v) { int i,k; k=0; for(i=1;i<=n;i++) ras[i]=a[i]; for(i=n;(i>=1)&&(k<m);i--) if (ras[i]>=v) { k++; } else ras[father[i]]+=ras[i]; if (k==m) return true; else return false; } int BinarySearch(int l,int r) { int mid; while(l<r) { mid=(l+r+1)>>1; if (!calc(mid)) r=mid-1; else l=mid; } return l; } int main() { int sum; while(scanf("%d%d%d",&n,&m,&a[1])!=EOF) { sum=a[1]; memset(as,0,sizeof(as)); for(i=2;i<=n;i++) { scanf("%d %d",&father[i],&a[i]); sum+=a[i]; } for(i=n;i>=1;i--) { as[i]+=a[i]; as[father[i]]+=as[i]; } printf("%d/n",BinarySearch(1,sum)); } }

Defend the Bases(http://acm.cs.ecnu.edu.cn/problem.php?problemid=2595)(二分加匹配)
题意：
  n个部队要开赴n个地方，每个部队到每个地方有个距离,求一个部队开赴指派方案，使最少的时间内所有地方都有部队到达。

解法:
 二分到达时间，判断是否能得到完美匹配，能就可行，否则不可行。
代码:
#include <iostream> #include<cmath> #include<algorithm> using namespace std; #define N 510 #define MAXN 500 #define M 510 #define INF 2147483647 int n, m; double g[N][M]; //记录边的权值 int match[M]; //记录匹配的点 int map[N][M]; bool chk[M]; int a[300][4],b[300][4]; double times[10005]; int dfs(int p) { int t,i,j; for(i=n+1;i<=n+m;i++) if (map[p][i]) if (!chk[i]) { chk[i]=1; t=match[i]; match[i]=p; if ((t==0)||(dfs(t))) return 1; match[i]=t; } return 0; } int xyl() { int i,res; res=0; memset(match,0,sizeof(match)); for(i=1;i<=n;i++) { memset(chk,0,sizeof(chk)); if (dfs(i)) res++; } return res; } void ReBuildGraph(int mid) { int i,j; memset(map,0,sizeof(map)); for(i=1;i<=n;i++) for(j=1;j<=m;j++) { if ( g[i][j]<=times[mid] ) map[i][n+j]=1; } } int BinarySearch(int l,int h) { if (l==h) return l; int mid=(l+h)/2; ReBuildGraph(mid); if (xyl()<n) return BinarySearch(mid+1,h); else return BinarySearch(l,mid); } int main() { int i,j,l; while(cin>>m>>n) { l=0; for(i=1;i<=m;i++) cin>>a[i][1]>>a[i][2]>>a[i][3]; for(i=1;i<=n;i++) cin>>b[i][1]>>b[i][2]; for(i=1;i<=n;i++) for(j=1;j<=m;j++) { g[i][j]=sqrt((double(b[i][1])-a[j][1])*(b[i][1]-a[j][1])+(b[i][2]-a[j][2])*(b[i][2]-a[j][2]))/a[j][3]; times[++l]=g[i][j]; } sort(times+1,times+1+n*m); printf("%.8lf/n",times[BinarySearch(1,n*m)]); } }