HDU 4610(数论+枚举状态)

若x=a0^p0*a1^p1*……*ai^pi(a0,a1……,ai是不同的素数)

对于条件3：因子和=(a0^0+a0^1+……+a0^p0)*(a1^0+a1^1+……+a1^p1)*……

如果是素数，显然质因子个位要为1，并且(a0^(p0+1)-1)/(a0-1)为素数

对于条件4：设因子积为a0^q0*a1^q1*……*ai^qi(a0,a1……,ai是不同的素数)

则有qi=(pi*pi+1)/2*(p0+1)*(p1+1)……(p(i-1)+1),判断所有的qi是不是偶数即可

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <map>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#include <cstdlib>
#include <algorithm>
#include <fstream>
#include <ctime>
#include <set>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
#define N 2000005
struct node
{
    int a,b,p,st;//数，个数，分值，状态
}num[1005];
struct Ys//约数  
{  
    int num,len;  
    int ys[N];  
    int cnt[N];  
}temp; 
int n,k;
int np[4];
int isprime[N];
int prime[N],sum;
bool cmp(node a,node b)
{
    return a.p>b.p;
}
void initisprime()//判断素数
{
    int i,j;
    for(i=2;i<N;i++)isprime[i]=1;
    for(i=2;i<N;i++)
    {
        for(j=2;i*j<N;j++)
            isprime[i*j]=0;
    }
}
void initprime()//筛选素数
{
    int i,j,flag;  
    prime[0]=2;  
    sum=1;  
    for(i=3;i<N;i++)  
    {  
        flag=1;  
        for(j=0;j<sum&&prime[j]*prime[j]<=i;j++)  
        {  
            if(i%prime[j]==0)  
            {  
                flag=0;break;  
            }  
        }  
        if(flag)prime[sum++]=i;  
    }  
}
void fj(int num)//分解质因数  
{  
    temp.num=num;  
    temp.len=0;  
    int i;  
    for(i=0;i<sum&&prime[i]*prime[i]<=num;i++)  
    {  
        if(num%prime[i]==0)  
        {  
            temp.ys[temp.len]=prime[i];  
            temp.cnt[temp.len]=1;  
            num/=prime[i];  
            while(num%prime[i]==0)temp.cnt[temp.len]++,num/=prime[i];  
            temp.len++;  
        }  
    }  
    if(num!=1)  
    {  
        temp.ys[temp.len]=num;  
        temp.cnt[temp.len]=1;  
        temp.len++;  
    }  
}
ll pm(ll x,ll y)
{
    ll ret=1;
    while(y--)ret*=x;
    return ret;
}
int getp(int x,int mm)
{
    int ret=0,cntd=1,sumd=0,sumj,i,j,flag,rst=0;
    if(isprime[x])rst|=1,ret++;
    //第一个条件
    fj(x);
    if(temp.len==1&&isprime[temp.cnt[0]+1])rst|=2,ret++;
    //第二个条件
    if(temp.len==1&&isprime[(pm(temp.ys[0],temp.cnt[0]+1)-1)/(temp.ys[0]-1)])rst|=4,ret++;
    //第三个条件
    flag=1;
    for(i=0;i<temp.len;i++)
    {
        if(!flag)break;
        sumj=(temp.cnt[i]+1)*temp.cnt[i]/2;
        if(sumj%2==0)continue;
        for(j=0;j<temp.len;j++)
        {
            if(i!=j)
            {
                if((temp.cnt[j]+1)%2==0)break;
            }
        }
        if(j==temp.len)flag=0;
    }
    if(flag)rst|=8,ret++;
    //第四个条件
    num[mm].st=rst;
    return ret;
}
int main()
{
    initisprime();
    initprime();
    int t,i,ans,j,temp,re,ti,tt;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&k);
        for(i=0;i<n;i++)scanf("%d%d",&num[i].a,&num[i].b);
        for(i=0;i<4;i++)scanf("%d",&np[i]);
        for(i=0;i<n;i++)num[i].p=getp(num[i].a,i);
        printf("%d",num[0].p);
        for(i=1;i<n;i++)printf(" %d",num[i].p);
        printf("\n");
        sort(num,num+n,cmp);
        ans=-inf;
        for(i=0;i<16;i++)
        {
            temp=0;
            tt=0;
            re=k;
            for(j=0;j<n;j++)
            {
                if((num[j].st&i)==0)
                {
                    if(re==0)break;
                    tt|=num[j].st;
                    if(re>=num[j].b)
                    {
                        temp+=num[j].p*num[j].b;
                        re-=num[j].b;
                    }
                    else
                    {
                        temp+=num[j].p*re;
                        re=0;
                    }
                    if(re==0)break;
                }
            }
            for(j=0;j<4;j++)
            {
                if((tt&(1<<j))==0)temp+=np[j];
            }
            if(re==0)ans=max(ans,temp);
        }
        printf("%d\n",ans);
    }
    return 0;
}