poj3685Matrix（二分）

看到题完全没想到是二分（太弱，太渣了），别人告诉本弱这道题是二分依旧不会做。
思路：二分答案ans，然后枚举j，然后再二分i，求出矩阵中小于等于ans元素的个数……
注意二分时如果l或r出现负数，最好直接向下取整，否则如(-1+0)/2=0，本来是想让他等于-1的，这样就可能进入死循环。

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <stdlib.h>
#include <stack>
#include <vector>
#include <string.h>
#include <math.h>
#include <queue>
#define msc(X) memset(X,-1,sizeof(X))
#define ms(X) memset(X,0,sizeof(X))
typedef long long LL;
using namespace std;
const LL inf=1ll<<33;
#define f(X,Y) ((X)*(X)+100000*(X)+(Y)*(Y)-100000*(Y)+(X)*(Y))
const double eps=1e-5;
int n;
LL m;
bool judge(LL num)
{
    LL cnt=0;
    for(int j=1;j<=n;j++)
    {
        int l=1,r=n+1,mid;
        while(l<r){
            mid=(l+r)>>1;
            if(f((LL)mid,(LL)j)>num) r=mid;
            else if(f((LL)mid,(LL)j)<=num&&f((LL)(mid+1),(LL)j)>num)
                break;
            else l=mid+1;
        }
        if(r!=mid)
            cnt+=mid;//当前<=num有多少个
    }
    return cnt<m;
}
int main(int argc, char const *argv[])
{
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%d %lld",&n,&m);
        LL l=-inf,r=inf;
        while(l<r){
            LL mid=(LL)(floor((l+r)/2.0+eps));
            if(judge(mid)) l=mid+1;
            else r=mid;
        }
        printf("%lld\n",l );
    }
    return 0;
}