51Nod 1014 X^2 Mod P

X*X mod P = A，其中P为质数。给出P和A，求<=P的所有X。

 

Input

两个数P A，中间用空格隔开。(1 <= A < P <= 1000000, P为质数)

Output

输出符合条件的X，且0 <= X <= P，如果有多个，按照升序排列，中间用空格隔开。
如果没有符合条件的X，输出：No Solution

Input示例

13 3

Output示例

4 9

#include <iostream>
#include <cstring>
#include <string>
#include <algorithm>
#include <stdio.h>
using namespace std;
#define ll long long
ll quickm(ll m,ll n,ll p)
{
    ll b=1;
    while(n){
        if(n&1) b=(b*m)%p;
        n=n>>1;
        m=(m*m)%p;
    }
    return b;
}
int main()
{
    ll p,a;
    scanf("%d%d",&p,&a);
    ll t=0,i;
    for(i=1;i<=p;i++){
        if(quickm(i,2,p)==a){
            if(t==0) printf("%d",i);
            else printf(" %d",i);
            t++;
        }
    }
    if(t==0) printf("No Solution
");
    else printf("
");
    printf("
");
    return 0;
}