[2018.05.05] T1 互质

题目描述

Shy 有 n 个数，问这 n 个数里最多有几个数两两之间互质。

输入输出格式

输入格式

第一行一个整数 n。

第二行 n 个数表示数组。

输出格式

输出一个数表示答案。

输入样例

5
1 2 3 4 5

输出样例

4

解题分析

感觉这道题好妙啊…根据数据范围得到算法…

考虑到最大的数据只有1000， 而每个数据只可能有一个33以上的因数。我们考虑分情况转移：

1.该数字不含33以上的质数： 直接暴力转移即可。

2.该数字含有33以上的质数： 对于这个质数建立一个分组背包再转移……

代码如下：

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cctype>
#include <cstring>
#include <cstdlib>
#define R register
#define IN inline
#define gc getchar()
#define W while
#define MX 10005
int pri[] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31};
int dp1[MX + 10], dp2[MX + 10], ext[MX + 10][50], cnt[MX + 10];
int dot;
template <class T>
IN void in(T &x)
{
    x = 0; R char c = gc;
    W (!isdigit(c)) c = gc;
    W (isdigit(c))
    x = (x << 1) + (x << 3) + c - 48, c = gc;
}
int main(void)
{
    int a, b, state;
    R int i, j, k;
    in(dot);
    for (i = 1; i <= dot; ++i)
    {
        in(a); state = 0;
        for (j = 0; j <= 10; ++j)
        {
            if(!(a % pri[j]))
            {
                state |= 1 << j;
                W (!(a % pri[j])) a /= pri[j];
            }
        }
        if(a == 1)//没有33以上的质因数， 暴力转移
        {
            for (j = 1; j <= 4500; ++j)
            if(!(j & state)) dp1[j | state] = std::max(dp1[j | state], dp1[j] + 1);
        }
        else ext[a][++cnt[a]] = state;//存入分组背包
    }
    for (i = 1; i <= 1000; ++i)//两个dp数组滚动一下
    {
        for (j = 1; j <= MX; ++j) dp2[j] = dp1[j];
        for (j = 1; j <= cnt[i]; ++j)
        for (k = 1; k <= 4500; ++k)
        if(!(k & ext[i][j])) dp2[k | ext[i][j]] = std::max(dp2[k | ext[i][j]], dp1[k] + 1);
        std::swap(dp1, dp2);
    }
    int mx = -1;
    for (i = 1; i <= MX; ++i) mx = std::max(mx, dp1[i]);
    printf("%d", mx);
}