#525. 【美团杯2020】平行四边形【原根】

题目链接

---

  既然x和y都是排列的话，我们不妨让x先升序从1～n如此输出，这样只用管y了。如果在不要求退化平行四边形的时候，我们可以用1, n, 2, n-1, 3,……这样的不断的两边互取的方式来完成，但是本题却要求求一个同时还要不满足退化的平行四边形的。

  于是这里引入了原根的思想，什么是原根，就是一个可以保证、、、……、他们（%p）之后，是1～p-1的排列，这p-1个元素互不相同，这题由于n+1是个质数，所以我们只需要对n+1取求一个原根就可以了。

  快速求原根的方式：从2到p-1进行枚举，找到一个数x满足：

对于，n有质因子、、……，满足对于这k个都有。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <bitset>
//#include <unordered_map>
//#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 1e3 + 7;
bool f[maxN];
int Prime[maxN], p_num;
inline void get_Prime()
{
    for(int i=2; i<maxN; i++)
    {
        if(!f[i])
        {
            f[i] = true;
            Prime[++p_num] = i;
        }
        for(int j=1; j<=p_num && i * Prime[j] < maxN; j++)
        {
            f[i * Prime[j]] = true;
            if(i % Prime[j] == 0) break;
        }
    }
}
inline ll fast_mi(ll a, ll b, ll mod)
{
    ll ans = 1;
    while(b)
    {
        if(b & 1) ans = ans * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return ans;
}
int N, P, elem[maxN], ele_cnt;
bool judege(int x)
{
    int k = N;
    ele_cnt = 0;
    for(int i=1; i<=p_num; i++)
    {
        bool flag = false;
        while(!(k % Prime[i]))
        {
            flag = true;
            k /= Prime[i];
        }
        if(flag)
        {
            elem[++ele_cnt] = Prime[i];
        }
        if(k == 1) break;
    }
    bool ok = true;
    for(int i=1; i<=ele_cnt; i++)
    {
        if(fast_mi(x, (P - 1) / elem[i], P) == 1)
        {
            ok = false;
            break;
        }
    }
    return ok;
}
signed main()
{
    get_Prime();
    int T; scanf("%d", &T);
    while(T--)
    {
        scanf("%d", &N); P = N + 1;
        int x;
        for(x = 2; ; x++)
        {
            if(judege(x))
            {
                break;
            }
        }
        for(int i=1, tmp = 1; i<=N; i++)
        {
            tmp = tmp * x % P;
            printf("%d %d\n", i, tmp);
        }
    }
    return 0;
}