[pe531]Chinese leftovers

题意：1e6~1e6+5000之间任意两个之间同余方程组的解。余数为欧拉函数。

解题关键：线性筛预处理，扩展中国剩余定理暴力求解。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<iostream>
#define MAXN 1005002
using namespace std;
typedef long long ll;
ll x,y,r[MAXN],m[MAXN],n;

//将求素数和欧拉函数值都线性解出
int prime[MAXN];//保存素数
bool vis[MAXN];//初始化
int phi[MAXN];//欧拉函数
void Lphisieve(int n){
    int cnt=0;
    for(int i=2;i<n;i++){
        if(!vis[i]){
            prime[cnt++]=i;
            phi[i]=i-1;// if p is prime,then phi[i]=i-1
        }
        for(int j=0;j<cnt&&i*prime[j]<n;j++){
            int k=i*prime[j];
            vis[k]=true;
            if(i%prime[j]==0){
                phi[k]=phi[i]*prime[j];
                break;
            }
            else   phi[k]=phi[i]*(prime[j]-1);
        }
    }
}

ll extgcd(ll a,ll b,ll &x,ll &y){
    ll d=a;
    if(b)   d=extgcd(b,a%b,y,x),y-=a/b*x;
    else  x=1,y=0;
    return d;
}
ll excrt(int n,ll *m,ll *r){
    ll M=m[0],pre=r[0],d;//a是模数 
    for(int i=1;i<n;i++){
        d=extgcd(M,m[i],x,y);
        if((pre-r[i])%d!=0) return -1;
        x=(pre-r[i])/d*x%m[i];
        pre-=x*M;
        M=M/d*m[i];//lcm 
        pre%=M;
    }
    return (pre%M+M)%M;
}
int main(){
    Lphisieve(1005001);
    ll sum=0;
    for(ll i=1000000;i<1005000;i++){
        m[0]=i,r[0]=phi[i];
        for(ll j=i+1;j<1005000;j++){
            m[1]=j,r[1]=phi[j];
            ll t=excrt(2,m,r);
            if(t==-1) t=0;
            sum+=t;    
        }
    }
    printf("%lld
",sum);
}