【？】高精版同余方程

就是一个蒟蒻无聊搞的乱七八糟的东西 想法很美好 结果真正打的时候发现有一堆要考虑的细节 几乎崩溃

主要就是判断一下为负时的各种情况 其实都差不多 应该有更好的方法 但蒟蒻暂时只能想到这种每种情况讨论一下的方法

然后死于在mzoj提交时我的数据格式 导致无数次的编译错误

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define rg register
const int N=1000,base=10000,power=4,K=100,P=9901;
char a[N],b[N],E[1]={'1'};
template <class t>void rd(t &x)
{
    x=0;int w=0;char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
    x=w?-x:x;
}

struct num{
    int a[N];
    bool f=0;
    num(){memset(a,0,sizeof(a));}
    num(char *s)
    {
        memset(a,0,sizeof(a));
        int len=strlen(s);
        a[0]=(len+power-1)/power;
        for(rg int i=0,t=0,w;i<len;++i,w*=10)
        {
            if(i%power==0) w=1,++t;
            a[t]+=w*(s[i]-'0');
        }
    }
    void add(int k) {if(k||a[0]) a[++a[0]]=k;}
    void re() {reverse(a+1,a+1+a[0]);}
    void print()
    {
        printf("%d",a[a[0]]);
        for(rg int i=a[0]-1;i>0;--i) printf("%0*d",power,a[i]);
    }
}p,q,x,y,e;

bool operator <(const num &p,const num &q)
{
    if(p.f&&!q.f) return 1;
    if(!p.f&&q.f) return 0;
    if(!p.f&&!q.f)
    {
        if(p.a[0]<q.a[0]) return 1;
        if(p.a[0]>q.a[0]) return 0;
        for(rg int i=p.a[0];i>0;--i)
        if(p.a[i]!=q.a[i]) return p.a[i]<q.a[i];
        return 0;
    }
    if(p.f&&q.f)
    {
        if(p.a[0]<q.a[0]) return 0;
        if(p.a[0]>q.a[0]) return 1;
        for(rg int i=p.a[0];i>0;--i)
        if(p.a[i]!=q.a[i]) return p.a[i]>q.a[i];
        return 0;
    }
}


num operator +(const num &p,const num &q)//高精+
{
    num c;
    c.a[0]=max(p.a[0],q.a[0]);
    for(rg int i=1;i<=c.a[0];++i)
    {
        c.a[i]+=p.a[i]+q.a[i];
        c.a[i+1]+=c.a[i]/base,c.a[i]%=base;
    }
    if(c.a[c.a[0]+1]) ++c.a[0];
    return c;
}

num operator -(const num &p,const num &q)
{
    num c=p;
    for(int i=1;i<=q.a[0];++i)
    {
        c.a[i]-=q.a[i];
        if(c.a[i]<0) c.a[i]+=base,--c.a[i+1];
    }
    while(c.a[0]>0&&!c.a[c.a[0]]) --c.a[0];
    return c;
}

num operator %(const num &p,const num &q)
{
    num x,y;
    for(int i=p.a[0];i>=1;--i)
    {
        y.add(p.a[i]),y.re();
        while(!(y<q)) y=y-q;
        y.re();
    }
    y.re(),y.f=p.f;
    return y;
}

num operator *(const num &p,const num &q)//高精*
{
    num c;
    c.a[0]=p.a[0]+q.a[0]-1;
    for(rg int i=1;i<=p.a[0];++i)
    for(rg int j=1;j<=q.a[0];++j)
    {
        c.a[i+j-1]+=p.a[i]*q.a[j];
        c.a[i+j]+=c.a[i+j-1]/base,c.a[i+j-1]%=base;
    }
    //while(c.a[0]>1&&!c.a[c.a[0]]) --c.a[0];
    while(c.a[c.a[0]+1]) ++c.a[0];
    return c;
}

num operator / (const num &p, const num &q)
{  
    num x, y;  
    for (int i = p.a[0];i >= 1;--i)                       //从最高位开始取数  
    {  
        y.add(p.a[i]);  //把数添到末尾（最低位），这时候是高位在前，低位在后  
        y.re();           //把数反过来，变为统一的存储方式：低位在前，高位在后  
        while ( !(y < q) )         //大于等于除数的时候，如果小于的话，其实答案上的该位就是初始的“0”  
            y = y - q, ++x.a[i];   //看能减几个除数，减几次，答案上该位就加几次。  
        y.re();                    //将数反过来，为下一次添数做准备  
    }  
    x.a[0] = p.a[0];
    while (x.a[0] > 0 && !x.a[x.a[0]]) --x.a[0]; 
    return x;  
}

void exgcd(num p,num q,num &x,num &y)
{
    if(!q.a[0]) {x=e,y=q;return;}
    exgcd(q,p%q,x,y);
    num t=x,tt=p/q;tt=tt*y,tt.f=y.f;
    x=y;
    if(t.f&&tt.f)
    {
        if(t<tt) y=tt-t,y.f=0;
        else y=t-tt,y.f=1;
    }
    else if(t.f&&!tt.f) y=t+tt,y.f=1;
    else if(!t.f&&tt.f) y=t+tt,y.f=0;
    else if(!t.f&&!tt.f)
    {
        if(!(t<tt)) y=t-tt,y.f=0;
        else y=tt-t,y.f=1;
    }
}

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    e=num(E);
    scanf("%s%s",&a,&b);
    reverse(a,a+strlen(a));
    reverse(b,b+strlen(b));//反转！！！ 
    p=num(a),q=num(b);
    exgcd(p,q,x,y);
    x=x%b;
    if(x.f) x=b-x;
    x=x%b;
    x.print();
    return 0;
}