POJ 3658 Matrix

　　双重二分,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( 双重二分+负数+死循环)的更多相关文章

  1. 51nod 1105 第K大的数 【双重二分/二分套二分/两数组任意乘积后第K大数】

    1105 第K大的数  基准时间限制:1 秒 空间限制:131072 KB 分值: 40 难度:4级算法题  收藏  关注 数组A和数组B,里面都有n个整数.数组C共有n^2个整数,分别是A[0] * ...

  2. 二分套二分 hrbeu.acm.1211Kth Largest

    Kth Largest TimeLimit: 1 Second   MemoryLimit: 32 Megabyte Description There are two sequences A and ...

  3. LA 3971 (二分) Assemble

    题意: 你有b块钱想要组装一台电脑.给出n个配件的种类,品质和价格,要求每个种类的配件各买一个总价格不超过b且“品质最差配件”的品质因子应尽量大. 这种情况下STL的map的确很好用,学习学习 这种最 ...

  4. 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 ...

  5. poj 3685 矩阵问题 查找第K小的值

    题意:N阶矩阵Aij= i2 + 100000 × i + j2 – 100000 × j + i × j,求第M小的元素. 思路:双重二分 考虑到,aij是跟着i递增的,所以i可以作为一个二分搜索 ...

  6. CDQ分治&整体二分学习个人小结

    目录 小结 CDQ分治 二维LIS 第一道裸题 bzoj1176 Mokia bzoj3262 陌上花开 bzoj 1790 矩形藏宝地 hdu5126四维偏序 P3157 [CQOI2011]动态逆 ...

  7. 整体二分QAQ

    POJ 2104 K-th Number 时空隧道 题意: 给出一个序列,每次查询区间第k小 分析: 整体二分入门题? 代码: #include<algorithm> #include&l ...

  8. JAVA源码走读（二）二分查找与Arrays类

    给数组赋值:通过fill方法. 对数组排序:通过sort方法,按升序.比较数组:通过equals方法比较数组中元素值是否相等.查找数组元素:通过binarySearch方法能对排序好的数组进行二分查找 ...

  9. 用c语言编写二分查找法

    二分法的适用范围为有序数列,这方面很有局限性. #include<stdio.h> //二分查找法 void binary_search(int a[],int start,int mid ...

随机推荐

  1. 【转】eclipse -- the project was not built due to a resource exists with a different case...

    原文网址:http://blog.csdn.net/mylinx/article/details/44280563 进行编码时,工程前面莫名有个红X,正当百思不得其解时,发现在[problems]下有 ...

  2. 数据结构（树套树）：ZJOI 2013 K大数查询

    有几个点卡常数…… 发现若第一维为位置,第二维为大小,那么修改时第一维修改区间,查询时第一维查询区间,必须挂标记.而这种情况下标记很抽象,而且Push_down不是O(1)的,并不可行. 那要怎么做呢 ...

  3. 软工UML学习札记

    UML模型由:事物.关系和图组成 (1)类(class)── 类用带有类名.属性和操作的矩形框来表示. (2)主动类(active class)── 主动类的实例应具有一个或多个进程或线程,能够启动控 ...

  4. Codeforces Round #335 (Div. 1)--C. Freelancer's Dreams 线性规划对偶问题+三分

    题意:p, q,都是整数. sigma(Ai * ki)>= p, sigma(Bi * ki) >= q; ans = sigma(ki).输出ans的最小值 约束条件2个,但是变量k有 ...

  5. linux内网机器访问外网代理设置squid

    公司一般出于安全考虑, 在同一局域网中只有一台机器可以访问外网,运维进行了整体的限制, 但是在后面的工作中,需要在机器上安装一些软件,及命令,所以其他的机器需要访问外网来简化工作, 但又不能打乱原有运 ...

  6. hive 模块

  7. JSP执行过程详解

    复习JSP的概念 JSP是Java Server Page的缩写,在传统的HTML页面中加入JSP标签和java的程序片段就构成了JSP. JSP的基本语法:两种注释类型.三个脚本元素.三个元素指令. ...

  8. [分享]WPF 虚拟键盘

    场景:用WPF做触屏的软件难免会需要用户输入的问题,至少是简单的数字,这个时候就免不了需要键盘输入. 思路:既然是虚拟键盘,那么我的目的就是模拟键盘输入即可. 1.模拟键盘输入 模拟键盘输入?那么肯定 ...

  9. 利用NIO建立Socket服务器

    传统的Java 的IO,利用Socket建立服务器,接收客户端连接,一般都是为每一个连接建立一个线程,如果连接数巨大,那么服务器开销也将巨大..NIO的原理,可以参照图:http://new.51ct ...

  10. padding与margin的差别

    之前一直没有搞懂android:padding和android:layout_margin的差别,事实上概念非常easy,padding是站在父view的角度描写叙述问题,它规定它里面的内容必须与这个 ...