【题解】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,并不是数值的位置连续 问题: 给出一个未排序的整数数组,找出最长的连续元素序列的长度 ...
随机推荐
- Java 实现多线程切换等待唤醒交替打印奇偶数
引言 在日常工作生活中,可能会有用时几个人或是很多人干同一件事,在java编程中,同样也会出现类似的情况,多个线程干同样一个活儿,比如火车站买票系统不能多个人买一到的是同一张票,当某个窗口(线程)在卖 ...
- js 拦截全局 ajax 请求
你是否有过下面的需求:需要给所有ajax请求添加统一签名.需要统计某个接口被请求的次数.需要限制http请求的方法必须为get或post.需要分析别人网络协议等等,那么如何做?想想,如果能够拦截所有a ...
- 前台页面获取servlet传过来的数据
servlet中的代码: public void doGet(HttpServletRequest request, HttpServletResponse response) throws Serv ...
- sublimetext3打造pythonIDE
虽然pycharm是非常好用的pythonIDE,用来开发项目很方便,但是修改调整单个或几个小程序就显得很笨重,这时候我们可以选择使用sublime. 一般来说要开发项目我都用pycharm,开发简单 ...
- Zookeeper协调分布式节点demo
多台服务器和客户端通过第三方组件Zookeeper管理 public class DistributedServer { private static final String connectStri ...
- IDEA破解 2017 IDEA license server 激活(可用)
进入ide主页面,help-register-license server,然后输入 http://idea.iteblog.com/key.PHP(注意:php要小写)即可~
- dubbo zookeeper案例
Alibaba有好几个分布式框架,主要有:进行远程调用(类似于RMI的这种远程调用)的(dubbo.hsf),jms消息服务(napoli.notify),KV数据库(tair)等.这个框架/工具/产 ...
- 简单实现接口自动化测试(基于python+unittest)
简单实现接口自动化测试(基于python+unittest) 简介 本文通过从Postman获取基本的接口测试Code简单的接口测试入手,一步步调整优化接口调用,以及增加基本的结果判断,讲解Pytho ...
- HDU 3466 01背包变形
给出物品数量N和总钱数M 对于N个物品.每一个物品有其花费p[i], 特殊值q[i],价值v[i] q[i] 表示当手中剩余的钱数大于q[i]时,才干够买这个物品 首先对N个物品进行 q-p的排序,表 ...
- kill -signal
1. SIGHUP 启动被终止的进程,可让该PID重新读取配置文件,类似于重启服务 对应的数字为1 9.SIGTERM 以正常的结束进程来终止进程 15.SIGSTOP 暂停一个进程相当于crtl+z