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

随机推荐

【HDOJ】3451 Beat drop
BFS.当水滴破裂飞溅后,直到碰到水滴才会停止(观察case1).同时,考虑当水滴飞溅到点(x,y)并且该点同一时间破裂的情况,该水滴算作吸收. /* 3451 */ #include <ios ...
HDOJ(HDU) 2309 ICPC Score Totalizer Software(求平均值)
Problem Description The International Clown and Pierrot Competition (ICPC), is one of the most disti ...
android关于图片缩放
网上有许多关于图片缩放的demo,本人都感觉不怎么好用,最近在github看到了一个简单的支持多指缩放图片的Android View类 gesture-imageview (地址:https://g ...
postgresql使用文档之一初始化数据存储区
17.2. 创建一个数据库集群(Database Cluster) 在你能做任何事情之前,你必须在磁盘上初始化一块存储空间.我们称这为一个数据库集群(database cluster). 一个Data ...
selenium webdriver启动IE浏览器失败的解决办法
通过selenium webdriver启动IE浏览器失败,报错:selenium.common.exceptions.WebDriverException: Message: Unexpected ...
单件模式Singleton来控制窗体被重复或多次打开
本文转载:http://blog.csdn.net/a0700746/article/details/4473796 一般在百度搜一下,会出来一下内容,看来很好用.Singleton很方便的一个用处就 ...
Swift-Dictionary
1.字典写法 Dictionary<KeyType,ValueType>,KeyType是你想要储存的键,ValueType是你想要储存的值. 唯一的限制就是KeyType必须是可哈希的, ...
Redis配置文件分析
#Redis演示示例配置文件 # 注意单位问题:当须要设置内存大小的时候,能够使用类似1k.5GB.4M这种常见格式: # # 1k=> 1000 bytes #1kb => 1024 b ...
Ext入门的第一个程序（1）
1.Ext是什么? extjs是集UI和ajax框架与一身的,界面又好看,又有很强的ajax交互功能,适合不会做漂亮页面的程序员用的,缺点就是太大了,要导入近800KB左右的js和css文件,这对于w ...
[转] 消息系统该Push/Pull模式分析
信息推拉技术简介 “智能信息推拉(IIPP)技术”是在网上信息获取技术中加入了智能成份,从而有助于用户在海量信息中高效.及时地获取最新信息,提高了信息系统主动信息服务的能力.如果引入基于IIPP的主 ...

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

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

随机推荐

热门专题