码迷,mamicode.com
首页 > 其他好文 > 详细

【美团杯2020】平行四边形

时间:2020-05-27 23:19:49      阅读:130      评论:0      收藏:0      [点我收藏+]

标签:bit   amp   ref   lag   nbsp   cto   name   style   view   

题:http://uoj.ac/problem/525

分析:1 . 答案是(i,g^i%(n+1))后一项为原根的次方。采用反证法,假设有4个点:(a,g^a),(b,g^b),(c,g^c),(d,g^d);

   2 . 若形成平行四边形则要满足:(1) a-b==c-d. (2)g^a-g^b==g^c-g^d   (    化简为 g^b ( g^(a-b) - 1 ) == g^d ( g^(c-d) -1)   )

   3 . 联立(1)(2)得必须a==c&&b==d才能形成平行四边形,而我们第一项的 i 避免了这一点,所以直接求n+1的原根

技术图片
#include<bits/stdc++.h>
#define lowbit(i) i&(-i)
#define pb push_back
using  namespace std;
typedef long long ll;
ll ksm(ll a,ll b,ll mod){
    ll t=1;
    while(b){
        if(b&1)
            t=t*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return t;
}
///求原根
ll solve(ll p,ll m){///p==phi(m)
    ll x=p;
    vector<ll>tmp;
    tmp.clear();
    for(ll i=2;i*i<=x;i++){
        if(x%i==0){
            tmp.pb(i);
            while(x%i==0)
                x/=i;
        }
    }
    if(x>1)
        tmp.pb(x);
    for(ll i=2;;i++){
        ll flag=1;
        for(auto x:tmp){
            if(ksm(i,p/x,m)==1){
                flag=0;
                break;
            }
        }
        if(flag)
            return i;
    }
    return -1;
}
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        ll n;
        scanf("%lld",&n);
        ll g=solve(n,n+1);
        if(g==-1){
            printf("-1\n");
        }
        else{
            for(ll i=1;i<=n;i++)
                printf("%lld %lld\n",i,ksm(g,i,n+1));
        }
    }
    return 0;
}
View Code

 

【美团杯2020】平行四边形

标签:bit   amp   ref   lag   nbsp   cto   name   style   view   

原文地址:https://www.cnblogs.com/starve/p/12977294.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!