PAT (Advanced Level) Practice 1010 Radix 二分法

一、概述

给出两个数和其中一个数的进制，求相等时另一个数的进制。

我做的时候就很奇怪，这进制也没给范围啊，以为默认是进制36，因为z代表35就按这个写了。结果错很多（这题的样例是真的多），看来坑点就在这个进制范围上。

要注意一点，对于不同的数字，进制的可能范围不同，进制的最小值为数字每位中的最大值加一，而最大值则可能是另外一个数的十进制值加一。因为如果进制超过了另外一个要比较的数加一，那就不可能选择它了。比如我十进制的10，你十二进制表示10必定是a，十一进制表示10也是啊；那么因为选择最小的进制，我肯定是选十一的。这样就确定了二分的上下限。

还有坑就是int全换为long long int，不然会溢出。

二、分析

这道题不能使用穷举做。实际上当我们理解了进制的范围时，我们就不可能选择穷举了。要用二分法做。使用string会很方便。

以N1的进制已知为例，先求N1的十进制值：

for(it=N1.rbegin();it!=N1.rend();it++)
{
	if(*it>='0'&&*it<='9')
    	Nindec+=(*it-'0')*pow(radix,bit);
	else
	    Nindec+=(*it-'a'+10)*pow(radix,bit);
	bit++;
}

使用倒序迭代器太方便了。

然后求N2各位的最大值：

for(it=N2.rbegin();it!=N2.rend();it++)
	{
		if(*it>='0'&&*it<='9')
		  {
		    if(*it-'0'>N2max)
		    N2max=*it-'0';
		  }
		else
		  {
		    if(*it-'a'+10>N2max)
		    N2max=*it-'a'+10;
		  }
	}

注意string上的每位是char类型，要减‘0’或‘a’+10。

之后二分的最大值最小值就可以确定了。

接下来是二分法：

while(min<=max)
	{
		long long int radix1=(min+max)/2;	
		int bit1=0;
		for(it1=N2.rbegin();it1!=N2.rend();it1++)
		{
			if(*it1>='0'&&*it1<='9')
				Otherdec+=(*it1-'0')*pow(radix1,bit1);
			else
				Otherdec+=(*it1-'a'+10)*pow(radix1,bit1);
			if(Otherdec>Nindec||Otherdec<0)
			break;
			bit1++;
		}
		if(Otherdec==Nindec)
		{
			equal=1;
			if(radix1<correct)
			{
				correct=radix1;
				max=radix1-1;
			}	
		}
		else if(Otherdec<0)
		  max=radix1-1;
		else if(Otherdec>Nindec)
			max=radix1-1;
		else
			min=radix1+1;
		Otherdec=0;
	}

有几点要注意：

当N2化为十进制期间，中间值大于N1或者小于0，那么都不用继续求了。直接往下在左子区间继续求就行。这个小于零溢出要注意。

循环条件是左小于等于右，而不是小于。这个要记住。

然后输出即可。

三、总结

一旦改成更大的类型，得分变多后，尽量把所有类型改大，不要偷懒。

二分法的端点要仔细思考。

PS：代码如下：

#include<stdio.h>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<math.h>
using namespace std;
int main()
{
	string N1,N2;
	long long int tag,radix;
	cin>>N1;
	cin>>N2;
	cin>>tag;
	cin>>radix;
	long long int Nindec=0;
	int equal=0;
	long long int correct;
	if(tag!=1)
	{
		string temp=N2;
		N2=N1;
		N1=temp;
	}
	string::reverse_iterator it;
	int bit=0;
	for(it=N1.rbegin();it!=N1.rend();it++)
	{
		if(*it>='0'&&*it<='9')
		Nindec+=(*it-'0')*pow(radix,bit);
		else
		Nindec+=(*it-'a'+10)*pow(radix,bit);
		bit++;
	}
	int N2max=0;
	for(it=N2.rbegin();it!=N2.rend();it++)
	{
		if(*it>='0'&&*it<='9')
		  {
		    if(*it-'0'>N2max)
		    N2max=*it-'0';
		  }
		else
		  {
		    if(*it-'a'+10>N2max)
		    N2max=*it-'a'+10;
		  }
	}
	string::reverse_iterator it1;
	long long int Otherdec=0;
	long long int min=N2max+1;
	long long int max;
	if(Nindec>min)
	max=Nindec+1;
	else
	max=min+1;
	while(min<=max)
	{
		long long int radix1=(min+max)/2;	
		int bit1=0;
		for(it1=N2.rbegin();it1!=N2.rend();it1++)
		{
			if(*it1>='0'&&*it1<='9')
				Otherdec+=(*it1-'0')*pow(radix1,bit1);
			else
				Otherdec+=(*it1-'a'+10)*pow(radix1,bit1);
			if(Otherdec>Nindec||Otherdec<0)
			break;
			bit1++;
		}
		if(Otherdec==Nindec)
		{
			equal=1;
			if(radix1<correct)
			{
				correct=radix1;
				max=radix1-1;
			}	
		}
		else if(Otherdec<0)
		  max=radix1-1;
		else if(Otherdec>Nindec)
			max=radix1-1;
		else
			min=radix1+1;
		Otherdec=0;
	}
	if(equal==1)
	{
		printf("%d",correct);
	}
	else
		printf("Impossible");
}