bc#54 div2

用小号做的div2

 

A：竟然看错了排序顺序。。。白白WA了两发

注意读入一整行（包括空格）：getline(cin,st)　　　　【gets也是资瓷的

#include<iostream>
using namespace std;
string s[1000];
int d[1000];
int T,N;

void qsort(int l,int r)
{
    int i=l,j=r;
    int m=(l+r)/2;
    int mid=d[m];
    while(i<j)
    {
        while(d[i]<mid)    i++;
        while(d[j]>mid)    j--;
        if(i<=j)
        {
            int tmp=d[i];
            d[i]=d[j];
            d[j]=tmp;
            string ts=s[i];
            s[i]=s[j];
            s[j]=ts;
            i++;
            j--;
        }
    }
    if(l<j)    qsort(l,j);
    if(i<r)    qsort(i,r);
}

int main()
{
    string st;
    cin>>T;
    while(T--)
    {
        cin>>N;  
        for(int i=0;i<=N;i++)
        {
            getline(cin,st);
            //cin>>s[i]>>d[i];
            int tl=st.length();
            int la=st[tl-1]-'0';
            int lb=st[tl-2]-'0';
            int lc=st[tl-3]-'0';
            int ld=st[tl-4]-'0';
            d[i]=ld*1000+lc*100+lb*10+la;
            //cout<<i<<" "<<d[i]<<endl;
            s[i]=st;
        }

        qsort(1,N);
        for(int i=N;i>=1;i--)
        {
            //cout<<s[i]<<endl;
            int tl=s[i].length();
            string tms=s[i];
            for(int j=0;j<=tl-6;j++)
                cout<<tms[j];
            cout<<endl;
        }
    }
}

 

B：自己YY的方法

其实就是找乘积的最小非质数因子

先判断无解的情况：所有数都是1 or 只有一个非1的数且是质数，剩下的都是1

若有解：对每个a[i]都分解质因子，找出每个a[i]最小的和第二小的质因子

然后把所有的质因子们放一起排序，取最小的两个相乘既是解。

#include<iostream>
#include<algorithm>
#include<cstring>
#include<ctime>
using namespace std;
#define LL long long
int T,N;
LL a[110];
LL sm[100100];

//**********************************************
// pollard_rho算法进行质因素分解
//*********************************************
long long factor[100];//质因素分解结果（刚返回时是无序的）
//注意factor里是每个质因数。比如N=8，那么factor[0..2]={2,2,2}   (2^3=8)
int tol;//质因素的个数，编号0~tol-1


const int S = 8; //随机算法判定次数，一般8~10就够了
//计算ret  = (a*b)%c
long long mult_mod(long long a,long long b,long long  c)
{
    a %= c;
    b %= c;
    long long ret = 0;
    long long tmp = a;
    while(b)
    {
        if(b & 1)
        {
            ret += tmp;
            if(ret > c)ret -= c;//直接取模慢很多
        }
        tmp <<= 1;
        if(tmp > c)tmp -= c;
        b >>= 1;
    }
    return  ret;
}
//计算  ret = (a^n)%mod
long long pow_mod(long long a,long long n,long long  mod)
{
    long long ret = 1;
    long long temp = a%mod;
    while(n)
    {
        if(n & 1)ret = mult_mod(ret,temp,mod);
        temp = mult_mod(temp,temp,mod);
        n >>= 1;
    }
    return  ret;
}
//通过  a^(n-1)=1(mod n)来判断n是不是素数
// n-1 = x*2^t中间使用二次判断
//是合数返回true,不一定是合数返回false
bool check(long long a,long long n,long long x,long long  t)
{
    long long ret = pow_mod(a,x,n);
    long long last = ret;
    for(int i = 1; i <= t; i++)
    {
        ret = mult_mod(ret,ret,n);
        if(ret == 1 && last != 1 && last != n-1)return  true;//合数
        last = ret;
    }
    if(ret != 1)return true;
    else return false;
}
//**************************************************
bool Miller_Rabin(long long  n)
//是素数返回true,(可能是伪素数)
//不是素数返回false
{
    if(n==1)  return true;
    if( n < 2)return false;
    if( n == 2)return true;
    if( (n&1) == 0)return false;//偶数
    long long x = n - 1;
    long long t = 0;
    while( (x&1)==0 )
    {
        x >>= 1;
        t++;
    }
    rand();/* *************** */
    for(int i = 0; i < S; i++)
    {
        long long a =  rand()%(n-1) + 1;
        if( check(a,n,x,t) )
            return false;
    }
    return true;
}


long long gcd(long long a,long long  b)
{
    long long  t;
    while(b)
    {
        t = a;
        a = b;
        b = t%b;
    }
    if(a >= 0)return  a;
    else return  -a;
}
//出一个因子
long long pollard_rho(long long x,long long  c)
{
    long long i = 1, k = 2;
    srand(time(NULL));
    long long x0 = rand()%(x-1) + 1;
    long long y = x0;
    while(1)
    {
        i ++;
        x0 = (mult_mod(x0,x0,x) + c)%x;
        long long d = gcd(y - x0,x);
        if( d != 1 && d != x)return  d;
        if(y == x0)return  x;
        if(i == k)
        {
            y = x0;
            k += k;
        }
    }
}
//对 n进行素因子分解，存入factor. k设置为107左右即可
void findfac(long long n,int  k)
{
    if(n == 1)return;
    if(Miller_Rabin(n))
    {
        factor[tol++] = n;
        return;
    }
    long long p = n;
    int c = k;
    while( p >= n)
        p = pollard_rho(p,c--);//值变化，防止死循环k
    findfac(p,k);
    findfac(n/p,k);
}



int main()
{
    cin>>T;
    while(T--)
    {
        cin>>N;
        int acc=0,occ=0;
        for(int i=1;i<=N;i++)
        {
            cin>>a[i];
            if(Miller_Rabin(a[i]))     acc++;
            if(a[i]==1)                 occ++;
        }
        if(((occ==N-1)&&(acc==1))||(occ==N))
            cout<<"-1"<<endl;
        else
        {
            int tot=0;
            for(int i=1;i<=N;i++)
            {
                tol = 0;
                memset(factor,0,sizeof(factor));
                findfac(a[i],107);
                //此时factor[0..tol-1]存放的就是所有质因数（无序）
                sort(factor,factor+tol);
                //cout<<a[i]<<" "<<tol<<endl;
                //cout<<factor[0]<<factor[1]<<factor[2]<<endl;
                if(factor[0]!=0)
                {
                    tot++;
                    sm[tot]=factor[0];
                }            
                if(factor[1]!=0)
                {
                    tot++;
                    sm[tot]=factor[1];
                }
            }
            sort(sm+1,sm+tot+1);
            cout<<sm[1]*sm[2]<<endl;
        }


    }
}