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. 设计模式之观察者模式（Observer Pattern）

    一.什么是观察者模式? 把现实世界中的报纸与订阅者的关系抽象出来就是观察者模式,一种报纸对应多个订阅者,订阅者可以随时解除订阅,未订阅的读者也可以随时开始订阅.一旦有新报纸发布,所有的订阅者都会收到新 ...

  2. C++ STL的各种实现版本

    ANSI/ISO的C++ STL规范版本正式通过以后,各个C++编译器厂商就可以依照标准所描述的原型去实现C++ STL泛型库,于是出现多种符合标准接口,但具体实现代码不同的泛型库,主要有: HP S ...

  3. 更改Visual Studio 2010/2012/2008的主题设置

    一.更改主题: 主题网站:http://studiostyl.es/ Visual Studio 2010发布也已经有一段时间了,不过安装后默认的白底的主题长时间看代码可能会感觉眼睛酸痛,况且时间长了 ...

  4. 实战weblogic集群之创建节点和集群

    一.启动weblogic,访问控制台 weblogic的domain创建完成后,接下来就可以启动它,步骤如下: $ cd /app/sinova/domains/base_domain/bin $ . ...

  5. Manacher Ural 1297 Palindrome

    1297. Palindrome Time limit: 1.0 secondMemory limit: 64 MB The “U.S. Robots” HQ has just received a ...

  6. UTF8与GBK、GB2312等其他字符编码的相互转换

    utf8与其他字符编码的转换是国际化必然遇到的问题. windows的api函数对: WideCharToMultiByte和MultiByteToWideChar 个人喜欢的跨平台库: 单纯做字符编 ...

  7. hdu 2795 线段树

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2795 #include <cstdio> #include <cmath> # ...

  8. Postman 基本操作学习

    History 所有使用postman发送的request都会保存在这里.点击之后会在当前Tab打开. 参考: Requests History Environments 这里用来设定当前reques ...

  9. windows路由命令route

    route print 查看当前的路由信息 route add 10.0.0.0 mask 255.0.0.0 10.1.1.1 增加一条到10.0.0.0/8网络的路由,网关是10.1.1.1 ro ...

  10. Java 中的泛型详解-Java编程思想

    Java中的泛型参考了C++的模板,Java的界限是Java泛型的局限. 2.简单泛型 促成泛型出现最引人注目的一个原因就是为了创造容器类. 首先看一个只能持有单个对象的类,这个类可以明确指定其持有的 ...