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
;
int cnt
;
}temp;
int n,k;
int np[4];
int isprime
;
int prime
,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;
}