【题解】CF1103D Professional layer
【题解】CF1103DProfessional layer
神题做前先\(orzyyb\)
一个很好的性质(之前也见过但是没有想到的)
zhengchu
\(gcd\le 10^{12}\) 所以不同的质因数\(\le 12\)
所以对这\(12\)个质因数状压。
所以答案显然小于等于\(12\)
对于每个数,只有因数和全局\(gcd\)有关的才有用,其余没用。
所以每种数(意会一下,就是和\(gcd\)关系一样的数)最多只需要\(12\)个,所以可以把\(n\)将下来。
然后就愉快地\(dp\)了(说的轻巧)
设\(dp(i,STATE)\)表示已经选用了\(i\)个数,现在\(gcd\)的质因数还剩\(STATE\)的状态
假如一个状态转移了\(\ge 13\)次,那么一定是不优的。开个数组记录顺便剪枝。转移方程
\]
\]
实际上\(13\)可以优化到\(d\)全局\(gcd\)的质因数个数\(m\)
有一个状压的技巧,如何枚举子集?$j \subset k $
- 知道\(j\)
for(register int k=(j+1)|j;k<=S;k=(k+1)|j)
- 知道\(k\)
for(register int j=(k-1)&k;j;j=(j-1)&k)
上抄zsy的代码
#include<bits/stdc++.h>
using namespace std;typedef long long ll;
#define DRP(t,a,b) for(register ll t=(a),edd=(b);t>=edd;--t)
#define RP(t,a,b) for(register ll t=(a),edd=(b);t<=edd;++t)
#define ERP(t,a) for(register ll t=head[a];t;t=e[t].nx)
#define pu(x) seg[x]=seg[(x)<<1]+seg[(x)<<1|1]
#define lef l,mid,pos<<1
#define rgt mid+1,r,pos<<1|1
#define all 1,n,1
#define re register
#define midd register int mid=(l+r)>>1
#define TMP template < class ccf >
TMP inline ccf qr(ccf b){
register char c=getchar();register int q=1;register ccf x=0;
while(c<48||c>57)q=c==45?-1:q,c=getchar();
while(c>=48&&c<=57)x=x*10+c-48,c=getchar();
return q==-1?-x:x;}
TMP inline ccf Max(ccf a,ccf b){return a<b?b:a;}
TMP inline ccf Min(ccf a,ccf b){return a<b?a:b;}
TMP inline ccf Max(ccf a,ccf b,ccf c){return Max(a,Max(b,c));}
TMP inline ccf Min(ccf a,ccf b,ccf c){return Min(a,Min(b,c));}
TMP inline void Swap(ccf& a,ccf& b){re ccf c=b;b=a;a=c;}
TMP inline void READ(ccf* _arr,int _n){RP(t,1,_n)_arr[t]=qr((ccf)1);}
//----------------------template&IO---------------------------
const int maxn=1e6+15;
//const ll mod=1e9+7;
inline ll gcd(ll x,ll y){
for(register ll t=x;y;t=y,y=x%y,x=t);
return x;
}
int n,m,S,e[maxn],vis[1<<12];
ll d,k,a[maxn],p[maxn],f[13][1<<12],ans,inf;
map < ll ,vector<int> >mp;
int main(){
#ifndef ONLINE_JUDGE
freopen("A.in","r",stdin);
freopen("A.out","w",stdout);
#endif
n=qr(1);k=qr(1ll);
RP(t,1,n) d=gcd(a[t]=qr(1ll),d);
RP(t,1,n) e[t]=qr(1ll);
for(register ll t=2;t*t<=d;++t){
if(d%t==0){
p[m++]=t;
while(d%t==0) d/=t;
}
}
if(d>1) p[m++]=d;
S=(1<<m)-1;
RP(t,1,n){
register ll temp=1;
RP(k,0,m-1){
while(a[t]%p[k]==0) a[t]/=p[k],temp*=p[k];
}
mp[temp].push_back(e[t]);
}
memset(f,63,sizeof f);
ans=inf=f[0][0];f[0][0]=0;
for(register auto pr:mp){
ll x=pr.first;sort(pr.second.begin(),pr.second.end());
if((int)pr.second.size()>m) pr.second.resize(m);
RP(t,0,S){
ll y=x,z=1;
RP(i,0,m-1)
if(t>>i&1)
while(y%p[i]==0) y/=p[i],z*=p[i];
vis[t]=(z<=k);
}
for(register auto cost:pr.second){
register bool fg=0;
for(register int i=m-1;~i;--i)
for(register int j=0;j<=S;++j)
if(f[i][j]<inf)
for(register int k=(j+1)|j;k<=S;k=(k+1)|j)
if(vis[k^j])
if(f[i+1][k]>f[i][j]+cost)
fg=1,(f[i+1][k]=f[i][j]+cost);
if(not fg) break;
}
}
RP(t,0,m)if(f[t][S]<inf) ans=Min(ans,f[t][S]*t);
if(ans==inf) puts("-1");
else cout<<ans<<endl;
return 0;
}
// orz zsy
【题解】CF1103D Professional layer的更多相关文章
- CF1103D Professional layer 状压DP
传送门 首先对于所有数求gcd并求出这个gcd含有的质因子,那么在所有数中,只有这一些质因子会对答案产生影响,而且对于所有的数,每一个质因子只会在一个数中被删去. 质因子数量不会超过\(11\),所以 ...
- CF1103D Professional layer dp
正解:dp 解题报告: 传送门! 首先不难想到求个gcd,然后把gcd质因数分解成p1w1*p2w2*p3w3*...*pmwm 显然只要满足对每个p有一个ai%pj!=0就好,也就是说对每个pj找出 ...
- CF1103D Codeforces Round #534 (Div. 1) Professional layer 状压 DP
题目传送门 https://codeforces.com/contest/1103/problem/D 题解 失去信仰的低水平选手的看题解的心路历程. 一开始看题目以为是选出一些数,每个数可以除掉一个 ...
- Professional layer CodeForces - 1103D (状压,gcd)
大意: 给定$n$元素序列$a$, 现在想要让$gcd(a_1,a_2,...,a_n)=1$. 对于每个$a_i$可以除以一个不超过$k$的因子, 代价为$e_i$, 假设一共选择了$x$个元素去除 ...
- CodeForces 1103D. Professional layer
题目简述:给定$1 \leq n \leq 10^6$个正整数$1 \leq a_i \leq 10^{12}$,修改第$i$个正整数$a_i$的花费为$1 \leq e_i \leq 10^9$,以 ...
- 题解 CF656G 【You're a Professional】
又是一道假黑题 它教会我们不要看难度标签 虽然这道题的数据范围很小,用cin能过,但我还是要讲一下快读 快读嘛,顾名思义,就是 快速读入 的意思 有同学就会问了,快速读入的原理是什么? 答:它的原理其 ...
- HDU 4374 One hundred layer DP的单调队列优化
One hundred layer Problem Description Now there is a game called the new man down 100th floor. The ...
- 【CF486E】LIS of Sequence题解
[CF486E]LIS of Sequence题解 题目链接 题意: 给你一个长度为n的序列a1,a2,...,an,你需要把这n个元素分成三类:1,2,3: 1:所有的最长上升子序列都不包含这个元素 ...
- leetcode & lintcode 题解
刷题备忘录,for bug-free 招行面试题--求无序数组最长连续序列的长度,这里连续指的是值连续--间隔为1,并不是数值的位置连续 问题: 给出一个未排序的整数数组,找出最长的连续元素序列的长度 ...
随机推荐
- springBoot 整合 mybatis+Oracle
现在比较流行的操作数据库操作层框架Mybatis,下面我们就来看看Springboot如何整合mybatis, 之前一直在用xml形式写sql,这次依然用xml的方式感觉这种还是比较灵活方便. 添加m ...
- NOI模拟题4 Problem B: 小狐狸(fox)
Solution 考虑分开统计朝向每一个方向的所有狐狸对答案的贡献. 比如说以向右为例, 我们用箭标表示每一只狐狸的方向, 用\('\)表示当前一步移动之前的每一只狐狸的位置. \[ \begin{a ...
- saltstack安装+基本命令
环境: node1:172.16.1.60 OS:centos 7.3 master hostname:centos7u3-1 node2:172.16.1.61 OS:centos 7.3 mini ...
- django导出数据到excel
import xlwt,StringIodef dumpData(request): #获取数据 activitys = Activity.objects.all().order_by('id') i ...
- BF3 里面的z cull reverse reload
Bf3 siggraph2011的 分享 http://advances.realtimerendering.com/s2011/White,%20BarreBrisebois-%20Renderin ...
- Spark原理小总结
1.spark是什么? 快速,通用,可扩展的分布式计算引擎 2.弹性分布式数据集RDD RDD(Resilient Distributed Dataset)叫做分布式数据集,是Spark中最基本的数据 ...
- 转:java的各个拓展类库的推荐方案
from: 链接:https://www.zhihu.com/question/21142149/answer/109854408 Java是世界上最强大的编程语言之一,很多开发人员和大型企业都偏爱J ...
- leetcode笔记:Ugly Number II
一. 题目描写叙述 Write a program to find the n-th ugly number. Ugly numbers are positive numbers whose prim ...
- PHP执行linux系统命令
本文是第一篇,讲述如何在PHP中执行系统命令从而实现一些特殊的目的,比如监控服务器负载,重启MySQL.更新SVN.重启Apache等.第二篇<PHP监控linux服务器负载>:http: ...
- apue学习笔记(第十四章 高级I/O)
本章涵盖了从多概念和函数:非阻塞I/O.记录锁.I/O多路转换.异步I/O.readv和writev函数以及存储映射I/O 非阻塞I/O 非阻塞I/O使我们可以发出open.read和write这样的 ...