庞果网--不可表示的数

题目：

    给定表达式[x/2] + y + x * y, 其中x,y都是正整数。其中的中括号表示下取整，例如[3/2] = 1 ， [5/2]  = 2。 有些正整数可以用上述表达式表达出来，例如正整数2，当取x = y = 1时，可以把2表达出来 （ 解释下：当x=y=1时， [x / 2] + y + x * y = [1 / 2] + 1 + 1 * 1 = 0+1+1 = 2 ）； 有些数可以有多种方式表达，例如13可以由 x = 2 y = 4 以及x = 3 y = 3来表示； 有些数无法用这个表达式表达出来，比如3。 从1开始第n个不能用这个表达式表示出来的数，我们叫做an，例如a1=1 a2=3,给定n，求an。 输入：n值 1<=n<=40 输出：an % 1000000007的结果（因为结果较大，输出an %1000000007的结果）。 函数头部 C/C++： int givean(int n); java： public class Main {       public static int givean(int n) {     } } 

代码：

int prim_givean[]={
2,3,5,7,13,17,19,31,61,89,
107,127,521,607,1279,2203,2281,3217,4253,4423,
9689,9941,11213,19937,21701,23209,44497,86243,110503,132049,
216091,756839,859433,1257787,1398269,2976221,3021377,6972593,13466917,20996011
//,24036583,25964951,30402457,32582657,37156667,42643801,43112609,57885161
};//梅森素数
int givean(int n)
{
	int base=1;
	int bit=1;
	int prim_num=prim_givean[n-1]-1;
	while(prim_num>bit)
	{
		base=(base<<bit)%1000000007;
		prim_num-=bit;
	}
	base=(base<<prim_num)%1000000007;
	//printf("%d:%d,%d\n",n,prim_givean[n-1],((base)-1));
		


	return base-1;
}

又是一次不理想的解答～～