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

随机推荐

c#,if 分支语句,条件运算符
//输入整数a和b, //若a²+b²大于100,则输出a²+b²百位以上数字, //否则输出两数之和 /*Console.Write("请输入整数a:"); int a = in ...
ASP.NET页面不被缓存
设置ASP.NET页面不被缓存分类:ASP.NET (3531) (1) /// <summary> /// 设置页面不被缓存 /// </summa ...
ZOJ-2587-Unique Attack(最小割的唯一性)
题意: 求无向图最小割是否唯一分析: 1.我们先对原图求一次最大流 2.对残留网络,我们从S开始,找到所有所有S能到达的点:再从T开始,找出所有能到达T的点. 3.判断原网络中是否还有没有访问到的点 ...
Android 开发无线调试
在进行android开发的时候有时候总是连接不上设备,因为设备的USB进行频繁插拔之后会导致一定的损坏.所以在设备已经无法利用数据线进行连接开发工具的时候就需要有那么一个功能,那就是无线调试,是不是感 ...
excel时会弹出向程序发送命令时出现问题的提示框
出现此问题需要做两个操作来解决: 1.在开始所有程序中找到Microsoft Excel 2007的运行程序,右键选择属性,在兼容性标签将“以管理员身份运行此程序”的勾去掉. 2.在打开的Excel程 ...
hdoj 1513 Palindrome【LCS+滚动数组】
Palindrome Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total ...
python面向对象【进阶篇】
静态方法通过@staticmethod装饰器即可把其装饰的方法变为一个静态方法,什么是静态方法呢?其实不难理解,普通的方法,可以在实例化后直接调用,并且在方法里可以通过self.调用实例变量或类变量 ...
Lucene简介
1 lucene简介1.1 什么是lucene Lucene是一个全文搜索框架,而不是应用产品.因此它并不像www.baidu.com 或者google Desktop那么拿来就能用,它只是提供 ...
VS 制作安装包小窥
难得忙里偷闲,看到有关VS制作安装包,按下文小试一把,还行,比不上Installshield. 首先在打开 VS2010 > 文件 >新建项目创建一个安装项目 XXX 在“目 ...
使用easy_install安装numpy、pandas、matplotlib及各种第三方模块
倒腾了一晚上最终把题目中的环境配好了.以下简要说明.留作资料.并共享. 1.安装python. 在cmd中能进入python环境,通过把python路径加入到系统路径中就可以实现. 2.安装easy- ...

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

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

随机推荐

热门专题