URAL——1204（数论之线性方程）

题目地址：http://acm.timus.ru/problem.aspx?space=1&num=1204

解析：n=pq，然后x（x-1）=pqM（M可以为任何正整数），因此就是解(x1)p-(y1)q=1,或者（x2）q-（y2）p=1；其中x1*y1=M，x2*y2=M，其实就是线性方程。

#include <iostream>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <cstdio>
using namespace std;
typedef long long ll;
#define INF 0xfffffff
#define MAX(a,b) a>b?a:b
#define MIN(a,b) a>b?b:a
#define N 101
#define M 100000
short int flag[M];
int prime[M/10];//只需要把sqrt（10e9）以内的筛选出来即可。
int q=0;//素数的下标 
void sieve_prime()
{
	memset(flag,0,sizeof(flag));
	int i; 
	for(i=2;i*i<M;i++) {    //边筛边存 
 		if(flag[i]) continue;
 		prime[q++]=i;
 		for(int j=i*i;j<M;j+=i)
 		  flag[j]=1;
 	}
 	for(i;i<M;i++)   //存另一部分素数 
 	  if(!flag[i]) prime[q++]=i;
} 
void extend_gcd(int a,int b,int &d,int &x,int &y){  //用扩展欧几里得求ax+by=c满足的一个解
	if(!b){
		d=a;x=1;y=0;
	}
	else {
		extend_gcd(b,a%b,d,y,x);
		y-=x*(a/b);
	}
}
int main()
{
	int i,j,k,t;
	int m,n;
	int a,b;
	sieve_prime();
	cin>>t;
	while(t--)
	 {	 	 
	         cin>>n;
		 cout<<"0 1 ";
	         int st;
	 	for(i=0;i<q;i++)
	 	  if(n%prime[i]==0) break;
 	    st=prime[i];
 	    int ed=n/st;
 	    int  x,y,d;
 	    extend_gcd(st,ed,d,x,y);   //本题一定有解且为4个。
 	    if(x<0)
		 {
 		        x+=ed;
			y-=st;	
 		 } 
 	    ll ans=max(x*st,abs(y)*ed);
 	    extend_gcd(ed,st,d,x,y);
 	    if(x<0)
		 {
 		    x+=st;
			y-=ed;	
 		 } 
 	    ll ans0=max(x*ed,abs(y)*st);
 	    cout<<min(ans,ans0)<<" "<<max(ans0,ans)<<endl;
 	    
 	 }	
	return 0;
}