QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#692874#9248. An Easy Math ProblemWaO_o#WA 638ms5984kbC++202.4kb2024-10-31 15:09:182024-10-31 15:09:26

Judging History

This is the latest submission verdict.

  • [2024-10-31 22:36:43]
  • hack成功,自动添加数据
  • (/hack/1098)
  • [2024-10-31 22:13:58]
  • hack成功,自动添加数据
  • (/hack/1096)
  • [2024-10-31 22:00:43]
  • hack成功,自动添加数据
  • (/hack/1095)
  • [2024-10-31 15:09:26]
  • Judged
  • Verdict: WA
  • Time: 638ms
  • Memory: 5984kb
  • [2024-10-31 15:09:18]
  • Submitted

answer

#include <bits/stdc++.h>

using namespace std;

#define int long long
#define deg( x ) cout<<""#x"="<<x<<endl
#define endl '\n'
#define pll pair<int,int>
#define se second
#define fr first
#define double long double

const int N=1e5+10;

int gcd( int a , int b ){
    if( b==0 ) return a;
    return gcd( b , a%b );
}

int idx=0 , idn=0;

int fac[ N ] , num[ N ] , fan[ N ];

void dec( int x ){
    idx=0;

    //fac[ ++idx ]=1 , num[ idx ]=1;
    for( int i=2; i*i<=x; i++ ){
        if( x%i==0 ) fac[ ++idx ]=i , num[ idx ]=0;
        while( x%i==0 ){
            num[ idx ]+=1;
            x/=i;
        }
    }

    if( x!=1 ) fac[ ++idx ]=x , num[ idx ]=1;

    idn=1; fan[ 1 ]=1;

    for( int i=1; i<=idx; i++ ){
        int en=idn , te=1;
        int now=fac[ i ];

        for( int k=1; k<=num[ i ]; k++ ){
            te*=now;
            for( int j=1; j<=en; j++ ) fan[ ++idn ]=fan[ j ]*te;
        }
    }


    sort( fan+1 , fan+1+idn );

}

const int mod=998244353;

int exgcd( int a , int b , int &x, int &y ) //扩展欧几里得求同余方程解; x,y及为 x*a+y*b=1 的一组可行解;
{
    if( b==0 ){
        x=1,y=0;
        return a;
    }

    int com=exgcd( b , a%b , x , y );

    int _x=x , _y=y ;
    y=_x-(a/b)*_y; x=_y;

    return com;
}

int getinv( int x , int m ) // exgcd获取x在模m情况下的逆元
{
    int ans,y;
    exgcd( x , m , ans , y );
    ans=( ( ans%m )+m )%m;
    return ans;
}

int inv[ N ];

void init(  ){
    for( int i=1; i<=idn; i++ ){
        inv[ i ]=getinv( fan[ i ] , mod );
    }

}

void solve(){
    int n;
    cin>>n;
    //n=1e10;

    dec( n );

    init();
    int ans=0;

    bool ok=false;

    map< int , int > mp;
    for( int i=1; i<=idn; i++ ){
        int p=fan[ i ];
        //cout<<p<<" ";
        for( int j=i; j<=idn; j++ ){
            int q=fan[ j ];
            if( p*q > n ){
                ok=true;
                break;
            }
            else if( n%( p*q )==0 ){
                int te=p*inv[ j ]; te%=mod;
                if( mp[ te ]==0 ){
                    mp[ te ]=1;
                    ans++;
                }
            }
        }
        if( ok ) break;
    }

    cout<<ans<<endl;

}

signed main() {
    ios::sync_with_stdio( 0 );
    cout.tie( 0 ); cin.tie( 0 );

    int T=1;
    cin>>T;
    while( T-- ){
        solve();
    }

    return 0;
}

详细

Test #1:

score: 100
Accepted
time: 1ms
memory: 5676kb

input:

10
1
2
3
4
5
6
7
8
9
10

output:

1
2
2
3
2
5
2
4
3
5

result:

ok 10 lines

Test #2:

score: -100
Wrong Answer
time: 638ms
memory: 5984kb

input:

2000
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
6469693230
646969323...

output:

1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
1535
...

result:

wrong answer 1st lines differ - expected: '29525', found: '1535'