POJ3658Matrix( 双重二分+负数+死循环)

　　双重二分,wa了一下午,实在不太明白为啥一写二分就会进入死循环.

　　INF要设的大一些,本题设0x3f3f3f3f会wa.

　　本题有负数, 二分时(l+r)/2与(l+r)>>1的结果有所不同;

　　如 l=0,r=-1,则 (l+r)/2=0,而(l+r)>>1=-1,而我们需要的正确答案是-1,所以第一种二分必须写成(l+r)>>1

　　以下两种二分可过:

　　(l+r)>>1

/*

* FillName:     二重二分

* Created:     2016年04月02日 14时48分41秒 星期六

* Author:      Akrusher

*

*/

#include <cstdio>

#include <cstdlib>

#include <cstring>

#include <cmath>

#include <ctime>

#include <iostream>

#include <algorithm>

#include <string>

#include <vector>

#include <deque>

#include <list>

#include <set>

#include <map>

#include <stack>

#include <queue>

#include <numeric>

#include <iomanip>

#include <bitset>

#include <sstream>

#include <fstream>

using namespace std;

#define rep(i,a,n) for (int i=a;i<n;i++)

#define per(i,a,n) for (int i=n-1;i>=a;i--)

#define in(n) scanf("%d",&(n))

#define in2(x1,x2) scanf("%d%d",&(x1),&(x2))

#define inll(n) scanf("%I64d",&(n))

#define inll2(x1,x2) scanf("%I64d%I64d",&(x1),&(x2))

#define inlld(n) scanf("%lld",&(n))

#define inlld2(x1,x2) scanf("%lld%lld",&(x1),&(x2))

#define inf(n) scanf("%f",&(n))

#define inf2(x1,x2) scanf("%f%f",&(x1),&(x2))

#define inlf(n) scanf("%lf",&(n))

#define inlf2(x1,x2) scanf("%lf%lf",&(x1),&(x2))

#define inc(str) scanf("%c",&(str))

#define ins(str) scanf("%s",(str))

#define out(x) printf("%d\n",(x))

#define out2(x1,x2) printf("%d %d\n",(x1),(x2))

#define outf(x) printf("%f\n",(x))

#define outlf(x) printf("%lf\n",(x))

#define outlf2(x1,x2) printf("%lf %lf\n",(x1),(x2));

#define outll(x) printf("%I64d\n",(x))

#define outlld(x) printf("%lld\n",(x))

#define outc(str) printf("%c\n",(str))

#define pb push_back

#define mp make_pair

#define fi first

#define se second

#define SZ(x) ((int)(x).size())

#define mem(X,Y) memset(X,Y,sizeof(X));

typedef vector<int> vec;

typedef long long ll;

typedef pair<int,int> P;

const int dx[]={,,-,},dy[]={,,,-};

const ll INF=1e12; //INF要设的大一些,本题设0x3f3f3f3f会wa

const ll mod=1e9+;

ll powmod(ll a,ll b) {ll res=;a%=mod;for(;b;b>>=){if(b&)res=res*a%mod;a=a*a%mod;}return res;}

const bool AC=true;

int t,n;

ll m;

ll cal(ll i,ll j){ //返回矩阵a[i][j]的值,数组会超内存

    return     i*i+j*j+*i-*j+i*j;      //int可能会爆

}

ll C(ll x){ //此处改为int会莫名奇妙的wa,返回值也得是ll

    int l,r,mid2;

    ll cnt; //此处一定要是ll

    cnt=;

    rep(j,,n+){

    l=,r=n+;

    while(r>l){

     mid2=(r+l)/;

     if(cal(mid2,j)>=x) r=mid2; //本题要求的是第m大的值

     else l=mid2+;

     }

     cnt+=r-;

     }

     return cnt;

}

int main()

{

    ll ub,lb,mid;

    in(t);

    while(t--){

        in(n);

        inlld(m);

        lb=-INF,ub=INF;            //INF要设的大一些

        while(ub>lb){

            mid=(ub+lb)>>;        //有负数,不能写成(l+r)/2

            if(C(mid)>=m) ub=mid; //类似中位数的那个题

            else lb=mid+;

            }

        outlld(ub-);//最后要减1

        }

    return ;

}

下面这种也可以过

(l+r)/2

/*

* FillName:     二重二分

* Created:     2016年04月02日 14时48分41秒 星期六

* Author:      Akrusher

*

*/

#include <cstdio>

#include <cstdlib>

#include <cstring>

#include <cmath>

#include <ctime>

#include <iostream>

#include <algorithm>

#include <string>

#include <vector>

#include <deque>

#include <list>

#include <set>

#include <map>

#include <stack>

#include <queue>

#include <numeric>

#include <iomanip>

#include <bitset>

#include <sstream>

#include <fstream>

using namespace std;

#define rep(i,a,n) for (int i=a;i<n;i++)

#define per(i,a,n) for (int i=n-1;i>=a;i--)

#define in(n) scanf("%d",&(n))

#define in2(x1,x2) scanf("%d%d",&(x1),&(x2))

#define inll(n) scanf("%I64d",&(n))

#define inll2(x1,x2) scanf("%I64d%I64d",&(x1),&(x2))

#define inlld(n) scanf("%lld",&(n))

#define inlld2(x1,x2) scanf("%lld%lld",&(x1),&(x2))

#define inf(n) scanf("%f",&(n))

#define inf2(x1,x2) scanf("%f%f",&(x1),&(x2))

#define inlf(n) scanf("%lf",&(n))

#define inlf2(x1,x2) scanf("%lf%lf",&(x1),&(x2))

#define inc(str) scanf("%c",&(str))

#define ins(str) scanf("%s",(str))

#define out(x) printf("%d\n",(x))

#define out2(x1,x2) printf("%d %d\n",(x1),(x2))

#define outf(x) printf("%f\n",(x))

#define outlf(x) printf("%lf\n",(x))

#define outlf2(x1,x2) printf("%lf %lf\n",(x1),(x2));

#define outll(x) printf("%I64d\n",(x))

#define outlld(x) printf("%lld\n",(x))

#define outc(str) printf("%c\n",(str))

#define pb push_back

#define mp make_pair

#define fi first

#define se second

#define SZ(x) ((int)(x).size())

#define mem(X,Y) memset(X,Y,sizeof(X));

typedef vector<int> vec;

typedef long long ll;

typedef pair<int,int> P;

const int dx[]={,,-,},dy[]={,,,-};

const ll INF=1e12; //INF要设的大一些,本题设0x3f3f3f3f会wa

const ll mod=1e9+;

ll powmod(ll a,ll b) {ll res=;a%=mod;for(;b;b>>=){if(b&)res=res*a%mod;a=a*a%mod;}return res;}

const bool AC=true;

int t,n;

ll m;

ll cal(ll i,ll j){ //返回矩阵a[i][j]的值,数组会超内存

    return     i*i+j*j+*i-*j+i*j;      //int可能会爆

}

ll C(ll x){ //此处改为int会莫名奇妙的wa,返回值也得是ll

    int l,r,mid2;

    ll cnt; //此处一定要是ll

    cnt=;

    rep(j,,n+){

    l=,r=n+;

    while(r>l){

     mid2=(r+l)/;

     if(cal(mid2,j)>=x) r=mid2; //本题要求的是第m大的值

     else l=mid2+;

     }

     cnt+=r-;

     }

     return cnt;

}

int main()

{

    ll ub,lb,mid;

    in(t);

    while(t--){

        in(n);

        inlld(m);

        lb=-INF,ub=INF;            //INF要设的大一些

        while(ub-lb>){

            mid=(ub+lb)/;

            if(C(mid)>=m) ub=mid;

            else lb=mid;

            }

        outlld(lb);

        }

    return ;

}

下面这种做标记的二分也挺有意思

/*

* FillName:     二重二分

* Created:     2016年04月02日 14时48分41秒 星期六

* Author:      Akrusher

*

*/

#include <cstdio>

#include <cstdlib>

#include <cstring>

#include <cmath>

#include <ctime>

#include <iostream>

#include <algorithm>

#include <string>

#include <vector>

#include <deque>

#include <list>

#include <set>

#include <map>

#include <stack>

#include <queue>

#include <numeric>

#include <iomanip>

#include <bitset>

#include <sstream>

#include <fstream>

using namespace std;

#define rep(i,a,n) for (int i=a;i<n;i++)

#define per(i,a,n) for (int i=n-1;i>=a;i--)

#define in(n) scanf("%d",&(n))

#define in2(x1,x2) scanf("%d%d",&(x1),&(x2))

#define inll(n) scanf("%I64d",&(n))

#define inll2(x1,x2) scanf("%I64d%I64d",&(x1),&(x2))

#define inlld(n) scanf("%lld",&(n))

#define inlld2(x1,x2) scanf("%lld%lld",&(x1),&(x2))

#define inf(n) scanf("%f",&(n))

#define inf2(x1,x2) scanf("%f%f",&(x1),&(x2))

#define inlf(n) scanf("%lf",&(n))

#define inlf2(x1,x2) scanf("%lf%lf",&(x1),&(x2))

#define inc(str) scanf("%c",&(str))

#define ins(str) scanf("%s",(str))

#define out(x) printf("%d\n",(x))

#define out2(x1,x2) printf("%d %d\n",(x1),(x2))

#define outf(x) printf("%f\n",(x))

#define outlf(x) printf("%lf\n",(x))

#define outlf2(x1,x2) printf("%lf %lf\n",(x1),(x2));

#define outll(x) printf("%I64d\n",(x))

#define outlld(x) printf("%lld\n",(x))

#define outc(str) printf("%c\n",(str))

#define pb push_back

#define mp make_pair

#define fi first

#define se second

#define SZ(x) ((int)(x).size())

#define mem(X,Y) memset(X,Y,sizeof(X));

typedef vector<int> vec;

typedef long long ll;

typedef pair<int,int> P;

const int dx[]={,,-,},dy[]={,,,-};

const ll INF=1e12; //INF要设的大一些,本题设0x3f3f3f3f会wa

const ll mod=1e9+;

ll powmod(ll a,ll b) {ll res=;a%=mod;for(;b;b>>=){if(b&)res=res*a%mod;a=a*a%mod;}return res;}

const bool AC=true;

int t,n;

ll m;

ll cal(ll i,ll j){ //返回矩阵a[i][j]的值,数组会超内存

    return     i*i+j*j+*i-*j+i*j;      //int可能会爆

}

ll C(ll x){ //此处改为int会莫名奇妙的wa,返回值也得是ll

    ll l,r,mid2;

    ll cnt,res; //此处一定要是ll

    cnt=;

    rep(j,,n+){

    l=,r=n;

    res=;

    while(r>=l){

     mid2=(r+l)>>;//起标记的作用

     if(cal(mid2,j)<=x) {

     res=mid2;

     l=mid2+;

    }

     else r=mid2-;

     }

     cnt+=res;

     }

     return cnt;

}

int main()

{

    ll ub,lb,mid,ans;

    in(t);

    while(t--){

        in(n);

        inll(m);

        lb=-INF,ub=INF;        //INF要设的大一些

        ans=;

        while(ub>=lb){

            mid=(ub+lb)>>;        //有负数,不能写成(l+r)/2

            if(C(mid)>=m) {

                ans=mid;

                ub=mid-; //带标记的二分

            }

            else lb=mid+;

            }

        outll(ans);

        }

    return ;

}

POJ3658Matrix( 双重二分+负数+死循环)的更多相关文章

51nod 1105 第K大的数【双重二分/二分套二分/两数组任意乘积后第K大数】
1105 第K大的数基准时间限制:1 秒空间限制:131072 KB 分值: 40 难度:4级算法题收藏关注数组A和数组B,里面都有n个整数.数组C共有n^2个整数,分别是A[0] * ...
二分套二分 hrbeu.acm.1211Kth Largest
Kth Largest TimeLimit: 1 Second MemoryLimit: 32 Megabyte Description There are two sequences A and ...
LA 3971 (二分) Assemble
题意: 你有b块钱想要组装一台电脑.给出n个配件的种类,品质和价格,要求每个种类的配件各买一个总价格不超过b且“品质最差配件”的品质因子应尽量大. 这种情况下STL的map的确很好用,学习学习这种最 ...
Gym 101064 D Black Hills golden jewels 【二分套二分/给定一个序列，从序列中任意取两个数形成一个和，两个数不可相同，要求求出第k小的组合】
D. Black Hills golden jewels time limit per test 2 seconds memory limit per test 256 megabytes input ...
poj 3685 矩阵问题查找第K小的值
题意:N阶矩阵Aij= i2 + 100000 × i + j2 – 100000 × j + i × j,求第M小的元素. 思路:双重二分考虑到,aij是跟着i递增的,所以i可以作为一个二分搜索 ...
CDQ分治&整体二分学习个人小结
目录小结 CDQ分治二维LIS 第一道裸题 bzoj1176 Mokia bzoj3262 陌上花开 bzoj 1790 矩形藏宝地 hdu5126四维偏序 P3157 [CQOI2011]动态逆 ...
整体二分QAQ
POJ 2104 K-th Number 时空隧道题意: 给出一个序列,每次查询区间第k小分析: 整体二分入门题? 代码: #include<algorithm> #include&l ...
JAVA源码走读（二）二分查找与Arrays类
给数组赋值:通过fill方法. 对数组排序:通过sort方法,按升序.比较数组:通过equals方法比较数组中元素值是否相等.查找数组元素:通过binarySearch方法能对排序好的数组进行二分查找 ...
用c语言编写二分查找法
二分法的适用范围为有序数列,这方面很有局限性. #include<stdio.h> //二分查找法 void binary_search(int a[],int start,int mid ...

随机推荐

HDU 5396 Expression(DP+组合数)（详解）
题目大意: 给你一个n然后是n个数. 然后是n-1个操作符,操作符是插入在两个数字之间的. 由于你不同的运算顺序,会产生不同的结果. 比如: 1 + 1 * 2 有两种 (1+1)*2 或者 ...
HDOJ 2055 An easy problem
Problem Description we define f(A) = 1, f(a) = -1, f(B) = 2, f(b) = -2, - f(Z) = 26, f(z) = -26; Giv ...
高效算法——E - 贪心-- 区间覆盖
E - 贪心-- 区间覆盖题目链接:http://acm.hust.edu.cn/vjudge/contest/view.action?cid=85904#problem/E 解题思路: 贪心思想, ...
VirtualBox扩展Ubuntu磁盘空间
有时候我们在使用virtualBox虚拟机时,创建虚拟机时并没有考虑到所占用硬盘大小,后来可能磁盘空间不够用了. 以下方法可以帮你扩展虚拟机的存储空间(以下Ubuntu为例). 1. 在宿主机器进入V ...
[JAVA关键字] synchronized
synchronized, Example: public synchronized void XXX() {} 参考 http://wenku.baidu.com/link?url=ecb1Zivf ...
组合数学：容斥原理(HDU1976)
●容斥原理所研究的问题是与若干有限集的交.并或差有关的计数. ●在实际中, 有时要计算具有某种性质的元素个数. 例: 某单位举办一个外语培训班, 开设英语, 法语两门课.设U为该单位所有人集合, A, ...
SRM 390(1-250pt)
DIV1 250pt 题意:给定整数n和k,问最少需要多少个n连接在一起形成的新整数t,使得t是k的倍数.如果不能形成倍数,输出-1.k <= 10^5,n <= 10^9. 解法:不断增 ...
Hive 创建和生成Rcfile 和SequenceFile格式的表
rcfile格式表需要从原始的textfile 文件格式表导出数据并导入到新建好的rcfile格式表里转: HIVE RCFile高效存储结构创建rcfile 格式的表: CREATE TABLE ...
[置顶] C语言单元测试框架
unitest.h /****************************************************************************** * * * This ...
NYOJ 491 幸运三角形
幸运三角形 import java.util.Scanner; public class Main{//搜索 static int sum=0,n; public static void main(S ...

POJ3658Matrix( 双重二分+负数+死循环)

POJ3658Matrix( 双重二分+负数+死循环)的更多相关文章

随机推荐

热门专题