[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=;
     for(int i=;i<n;i++){
         if(!vis[i]){
             prime[cnt++]=i;
             phi[i]=i-;// if p is prime,then phi[i]=i-1
         }
         for(int j=;j<cnt&&i*prime[j]<n;j++){
             int k=i*prime[j];
             vis[k]=true;
             if(i%prime[j]==){
                 phi[k]=phi[i]*prime[j];
                 break;
             }
             else   phi[k]=phi[i]*(prime[j]-);
         }
     }
 }

 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=,y=;
     return d;
 }
 ll excrt(int n,ll *m,ll *r){
     ll M=m[],pre=r[],d;//a是模数
     for(int i=;i<n;i++){
         d=extgcd(M,m[i],x,y);
         if((pre-r[i])%d!=) return -;
         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();
     ll sum=;
     for(ll i=;i<;i++){
         m[]=i,r[]=phi[i];
         for(ll j=i+;j<;j++){
             m[]=j,r[]=phi[j];
             ll t=excrt(,m,r);
             if(t==-) t=;
             sum+=t;
         }
     }
     printf("%lld\n",sum);
 }