【美团杯2020】平行四边形(原根)

【美团杯2020】平行四边形

蒜斜非常喜欢下围棋。自从AlphaOg面世以来，他就立志一定要研究出AlphaOg的破绽。 终于，他发现当AlphaOg遇到一种特殊局面后，它的神经网络会自动输出“投降”！

随着进一步的研究，蒜斜发现这种局面有着更一般的特性，不仅仅局限于固定大小棋盘。 具体来说，当棋盘大小是 n（n+1 是一个质数）且棋盘上恰好有 n 个棋子的时候，如果这些棋子的位置满足下列条件，那么 AlphaOg 就会直接投降。假设第 i 个棋子的位置是点 Pi，处在第 xi 行第 yi 列，那么这些坐标需要满足：

x1 至 xn 是 1−n 的排列。
y1 至 yn 是 1−n 的排列。
这些点之间不构成平行四边形（包括退化）。即对于任何两个不完全相同的棋子对 (Pa,Pb),(Pc,Pd)（允许它们之间共用至多一个棋子），线段 PaPb 与 PcPd 要么长度不同，要么所在的直线不平行且不重合。
凭借这项发现，蒜斜荣获了“北大算协吉祥物”的称号。 如果你也能找出一种合法方案，蒜斜的称号就是你的了！

输入格式
输入第一行包含一个整数 t(1≤t≤10)，表示数据组数。

对于每组数据，输入第一行包含一个整数 n(4≤n≤1000)，保证 n+1 是一个质数。

输出格式
对于每组数据，如果无解输出一行一个整数 -1。否则输出 n 行，每行两个整数 (xi,yi)(1≤xi≤n,1≤yi≤n)，表示第 i 个棋子的坐标。如果坐标方案不唯一，你只需要输出任意一种。

样例一
input
1
4

output
1 1
3 2
4 3
2 4

限制与约定
Small Task: n≤12。

Large Task: n≤1000。

时间限制：1s
空间限制：512MB

一个n*n棋盘(n+1为奇质数)每行和每列有且只有一个棋子，且没有四个棋子可以构成一个平行四边形，求任意一组解

• 原根简单介绍：
  • 如果m为奇质数一定存在模m的原根.
  • 原根的性质:假设一个数g是m的原根，那么gⁱ mod m的结果两两不同，且有 1
  • 原根的求法：如果g是m的原根，gⁱ = 1 (mod m)当且仅当i=m-1的时候成立
• 先求出模n+1的原根g，答案就是(1,g mod(n+1)),(2,g² mod(n+1)),(3,g³ mod(n+1))…(n,gⁿ mod(n+1))
• 正确性证明：
  • 假设四个可以构成平行四边形的点为（a,g^a）(b,g^b),(c,g^c),(d,g^d)
  • 四个点构成平行四边形，当且仅当满足以下两个条件
  • 1、a-b=c-d
  • 2、g^a-g^b=g^c-g^d <==> g^a(1-g^b-a) = g^c(1-g^d-c) <==> g^a=g^c
  • 根据原根的性质，当a!=c时g^a!=g^c,结论和假设矛盾
  • 所以就证明了不会构成平行四边形

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e3+9;
long long qui_pow(long long a,long long n,long long mod)
{
     
    long long ans = 1;
    while(n)
    {
     
        if(n&1)
            ans = ans*a%mod;
        a=a*a%mod;
        n>>=1;
    }
    return ans;
}
bool judge(long long a,long long n)
{
     
    for(long long i = 1; i<n; i++)
    {
     
        long long num = qui_pow(a,i,n)%n;
        if(i!=n-1&&num==1)
            return 0;
        if(i==n-1&&num!=1)
            return 0;
    }
    return 1;
}

int main()
{
     
    ios::sync_with_stdio(false);
    long long i,j,n,m,t,num;
    cin>>t;
    while(t --)
    {
     
        cin>>n;
        num = -1;
        for(i = 1; i<n; i++)
        {
     
            if(judge(i,n+1))
            {
     
                num = i;
                break;
            }
        }
        for(i = 1;i<=n;i++)
        {
     
            printf("%lld %lld\n",i,qui_pow(num,i,n+1)%(n+1));
        }
    }

    return 0;
}