Vijos 1921 严厉的班长 【状态压缩动态规划】

严厉的班长

描述

木姑娘在班级里面是班长。虽然是副班长,却有着比正班长更高的威信,并深受小朋友们的爱戴。

每天眼保健操时间,木姑娘都要监督所有小朋友认真做眼保健操。整个过程被描述为n个时间段,第i个时间段内对于揉捏穴位的速度被要求为a[i]。

如果某一位小朋友在n个时间段内,速度分别为b[1],b[2],...,b[n],则他的最终得分为 score = |a[1]-b[1]| + |a[2]-b[2]| + ... + |a[n]-b[n]|。

分数越大,则说明表现越差,如果b[1]到b[n]中,又满足任意两个数字均互素,则会受到木姑娘的严厉批评。

对于爱慕木姑娘的孩子来说,能和木姑娘有更多时间待在一起,才是最关键的。

格式

输入格式

第一行一个整数n,1<=n<=100。

第二行有n个整数,分别为a[1],a[2],...,a[n]且对于每一个a[i]满足1<=a[i]<=30。

输出格式

希望可以找出来一组b[1],...,b[n]满足score最小化,且b数组中的n个元素两两互素,并输出此刻的分数score。

样例1

样例输入1

5
1 1 1 1 1

Copy

样例输出1

0

Copy

样例2

样例输入2

5
1 6 4 2 8

Copy

样例输出2

3

Copy

限制

25%的数据,n<=10。
60%的数据,n<=30。
100%的数据,n<=100。

来源

感谢 小岛

题目链接:

  https://www.vijos.org/p/1921

题目大意:

  题目给一个N个数的数列Ai,数大小不超过30,要求给出一组解Xi,使得解Xi Xj两两互质的情况下和原先数列的绝对值的差的和最小(min Σ|Ai-Xi|),求这个和。

题目思路:

  【状压DP】

  首先考虑,因为要求两两互质,所以1肯定满足,那么取1<=Xi<2Ai-1(取2Ai-1和取1效果相同且更不易满足互质条件)

  然后看2x30-1=59以内的质数,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53 共16个。

  也就是说当n>16时,多出来的必须用1,那么考虑数的大小可知尽量用1取对应小的数,于是排序的时候按从大到小排序,只做前16个,后面的全是1(当然前面也可以用1而不用质数)

  所以问题转化为用16个质数的倍数去满足前16个数即可。将1~59质因数分解,记下每个数含有的质因数的二进制码,例如1为0,2为1,3为10,4为1,5为100,6为11 以此类推

  DP的时候f[i][j]表示前i个数,质因数已经被选走的状态为j的最小和(j为16位二进制数)

  考虑当前状态如果可行,就枚举当前的Xi,判断是否冲突,不冲突就更新新的f

  

 1 /****************************************************
 2
 3     Author : Coolxxx
 4     Copyright 2017 by Coolxxx. All rights reserved.
 5     BLOG : http://blog.csdn.net/u010568270
 6
 7 ****************************************************/
 8 #include<bits/stdc++.h>
 9 #pragma comment(linker,"/STACK:1024000000,1024000000")
10 #define abs(a) ((a)>0?(a):(-(a)))
11 #define lowbit(a) (a&(-a))
12 #define sqr(a) ((a)*(a))
13 #define mem(a,b) memset(a,b,sizeof(a))
14 const double EPS=0.00001;
15 const int J=10;
16 const int MOD=100000007;
17 const int MAX=0x7f7f7f7f;
18 const double PI=3.14159265358979323;
19 const int N=104;
20 const int M=65544;
21 using namespace std;
22 typedef long long LL;
23 double anss;
24 LL aans;
25 int cas,cass;
26 int n,m,lll,ans;
27 int p[16]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};
28 int a[N],fj[N];
29 int f[N][M];
30 void init()
31 {
32     int i,j,x;
33     for(i=1;i<59;i++)
34     {
35         x=i;
36         fj[i]=0;
37         for(j=0;j<16;j++)
38         {
39             while(x%p[j]==0)
40             {
41                 x/=p[j];
42                 fj[i]|=(1<<j);
43             }
44             if(x==1)break;
45         }
46     }
47 }
48 int main()
49 {
50     #ifndef ONLINE_JUDGE
51 //    freopen("1.txt","r",stdin);
52 //    freopen("2.txt","w",stdout);
53     #endif
54     int i,j,k;
55     int x,y,z;
56 //    for(scanf("%d",&cass);cass;cass--)
57 //    for(scanf("%d",&cas),cass=1;cass<=cas;cass++)
58     init();
59     while(~scanf("%d",&n))
60     {
61         mem(f,0x7f);
62         for(i=1;i<=n;i++)
63             scanf("%d",&a[i]);
64         sort(a+1,a+1+n,greater<int>());
65         f[0][0]=0;
66         for(i=1;i<=min(n,16);i++)
67         {
68             for(j=0;j<(1<<16);j++)
69             {
70                 if(f[i-1][j]==MAX)continue;
71                 for(k=1;k<=a[i]+a[i]-1;k++)
72                 {
73                     if(fj[k] & j)continue;
74                     f[i][j|fj[k]]=min(f[i][j|fj[k]],f[i-1][j]+abs(a[i]-k));
75                 }
76             }
77         }
78         ans=MAX;
79         for(i=0;i<(1<<16);i++)
80             ans=min(ans,f[min(n,16)][i]);
81         if(n>16)
82         {
83             for(i=17;i<=n;i++)
84                 ans+=a[i]-1;
85         }
86         printf("%d\n",ans);
87     }
88     return 0;
89 }
90 /*
91 //
92
93 //
94 */

时间: 2024-08-29 08:21:36

Vijos 1921 严厉的班长 【状态压缩动态规划】的相关文章

奇妙的算法—状态压缩动态规划

华电北风吹 天津大学认知计算与应用重点实验室 日期:2015/8/27 由于代码未调试完全正确论文草稿呈现 poj上一道需要用到状态压缩动态规划,链接http://poj.org/problem?id=3254 网上看到有很多人写出了代码,参考了一个带备忘的自顶向下的动态规划解法,自己写了一个有底向上的动态规划解法 有底向上: #include<iostream> #include<math.h> #include<ostream> #include<fstrea

状态压缩动态规划总结

状态压缩是一个很广的概念,在OI中也有很多的应用,当我们把他应用到动态规划中,可以用来精简状态,节约空间,也方便转移.最常见的就是用二进制来表是状态,利用各种位移运算,就可以实现\(O(1)\)的转移.状压DP适用于“窄棋盘”上的DP,否则状态太多无法存下. POJ1185 炮兵阵地 题意:给一个\(N \times M\)的地形盘,有平原和山坡,要求在平原上部署尽量多的炮(攻击范围2),使得不互相攻击. 数据范围:N <= 100:M <= 10,符合条件.如何表示状态?按行DP,一个二进制

状态压缩动态规划 -- 炮兵阵地

司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队,一个N*M的地图由N行M列组成,地图的每一格可能是山地(用"H" 表示),也可能是平原(用"P"表示),如下图.在每一格平原地形上最多可以布置一支炮兵部队(山地上不能够部署炮兵部队):一支炮兵部队在地图上的攻击范围如图中黑色区域所示: 如果在地图中的灰色所标识的平原上部署一支炮兵部队,则图中的黑色的网格表示它能够攻击到的区域:沿横向左右各两格,沿纵向上下各两格.图上其它白色网格均攻击不到.从图上可见炮兵的攻击

状态压缩动态规划 -- 棋盘问题 POJ 1321

一个 N * N 的棋盘上面,有些格子不能放,放置 M 的棋子, 每两个棋子不能在同一行或者同一列,问有多少种放法 DFS太慢,用SCR好点点 Python 只有 22 行,其实可以更短,但是得排成很长很长的一行 while True: table = [ [ 0 for j in range( 300 ) ] for i in range( 12 ) ] table[0][0] = 1 boardsize, chessnum = map( int, raw_input().split() )

状态压缩动态规划 -- 旅行商问题

旅行商问题: N个点(N<16)的带权有向图D,求一条路径,使得这条路经过每个点恰好一次, 并且路径上边的权值和最小(或者最大),或者求一条具有这样性质的回路. 状态压缩: 将二进制表示十进制数N的点集,比如: 10 = 0000000000001010 代表第1和3个点已经路过 18 = 0000000000010010 代表第1和4个点已经路过 一个整数就是一个点集, dp_arr[binary][to_]代表经过点集 binary 中,当前终点为to_, 且路径最短的值,若该状态不存在就是

状态压缩动态规划 -- 多米诺骨牌

用1*2 的骨牌通过组合拼成 m * n 的大矩形,问有几种拼法. 题目链接:http://poj.org/problem?id=2411 状态转移: 1.由于上一行的该列竖直放置骨牌为 0,影响到当前行的该列,当前行的该列为 1 2.当前行骨牌横放,上一行骨牌横放, 都为11 3.上一行该列置为 1,当前行当前列立着放为 0 #include <iostream> #include <cstring> using namespace std; #define MAXSIZE 12

状态压缩-动态规划

状态压缩-动态规划 ---By蒟蒻鱼 用二进制表示城市的到达(规划的)状态 每一个二进制数都代表一个唯一的十进制数 预备知识 位运算 优先级 ~ << 和 >> & ^ | 按位与 & 全一则一,否则为零 按位或 | 有一则一,否则为零 按位取反 ~ 是零则一,是一为零 按位异或 ^ 不同则一,相同为零 移位 >> << 基本运算: 集合取并 A|B 集合取交 A&B 集合相减 A&~B 集合相减 ALLBITS~A 置位 A

vijos - P1286座位安排 (DP状态压缩 + 组合数 + python)

P1286座位安排 Accepted 标签:[显示标签] 背景 快要期中考试了!老师需要hzy帮他排考试的座位... 描述 考场里的座位恰好有n行m列,并且恰好有n*m位考生在这个考场里面考试,也就是说,所有的座位上都有考生.hzy根据学校记载,有k位考生可能作弊,因此hzy不能让他们之中的任何两个人做在相邻的座位上!所谓相邻的座位,即在同一行相邻列或者在同一列的相邻行的座位.hzy准备这样安排座位,首先随机选择一种方案,如果这种方案是合法的,就用这种方案,否则重新选择.你的任务是计算,他得到一

vijos - P1456最小总代价 (状态压缩DP + 记忆化搜索)

P1456最小总代价 Accepted 标签:[显示标签] 描述 n个人在做传递物品的游戏,编号为1-n. 游戏规则是这样的:开始时物品可以在任意一人手上,他可把物品传递给其他人中的任意一位:下一个人可以传递给未接过物品的任意一人. 即物品只能经过同一个人一次,而且每次传递过程都有一个代价:不同的人传给不同的人的代价值之间没有联系: 求当物品经过所有n个人后,整个过程的总代价是多少. 格式 输入格式 第一行为n,表示共有n个人(16>=n>=2): 以下为n*n的矩阵,第i+1行.第j列表示物