Unsolvable CF 225E

题目链接：CodeForces - 225E

题目大意：

给出一个等式，z=[x/2]+y+x⋅y，求另这个等式无解的正整数z，从小到大排序，访问第i个。

思路：

先把取整去掉。

因为若x为偶数是不需要取整的，只有x为奇数时才需要取整，而且是向小的取整。

（1）对于奇数情况：

令x=2k+1;

z=k+y+(2k+1)y;（已经对x/2取过整）

z=k+2y(k+1);

z+1=(k+1)(2y+1);

因为k+1代表了所有的正整数，而2y+1代表了所有的奇数，所以两者的乘积一定含有质因子，所以无解的条件就是z+1=2^t，即z=2^t-1。

（2）对于偶数情况：

令x=2k;

z=k+y+2ky;

z=k+y(2k+1);

2z=2k+2y(2k+1);

2z+1=(2k+1)(2y+1);

这里的2k+1与2y+1都是奇数，乘积代表了所有的奇数。对于这里的z，一定满足x是奇数时候的情况，所以2z+1=2^(t+1)-1，发现这样的结果一定也是一个奇数，所以2^(t+1)-1只能是一个质数。

 

无解的z为2^(t+1)-1，其中，长成2^n-1这样的质数叫梅森质数（第一次听说！），目前有49个梅森素数，打表所有可能构成梅森素数的n(就是t+1，别忘了z=2^t-1,所以算结果的时候不是算梅森素数，而是算“类梅森素数”其中的n换成了n-1），然后用快速幂求解就可以啦~

代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<map>
#include<vector>
#include<set>
#include<string>
#include<cmath>
#define ll long long
using namespace std;
int n,m=1e9+7;
int a[50]={0,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};
ll mi(ll a,ll b)
{
    int ans=1;
    while(b)
    {
        if(b&1)
            ans=ans*a%m;
        a=a*a%m;
        b>>=1;
    }
    return ans;
}
int main()
{
    cin>>n;
    cout<<(mi(2,a[n]-1)-1+m)%m<<endl;//a[n]还要减1，因为n=t+1，所以t=n-1.+m%m是为了防止取负模。
    return 0;
}