洛谷 P7800 [COCI2015-2016#6] PAROVI

PS：如果读过题了可以跳过题目描述直接到题解部分

提交链接：洛谷 P7800 [COCI2015-2016#6] PAROVI

题目

题目描述

Mirko 和 Slavko 在玩一个游戏，先由 Mirko 在 1…N 中选出几组互质的数。例如当 N=5 时，Slavko 可以选择 {{1,2},{3,4},{2,5},{3,5},…} 中的几组。

然后轮到 Slavko。他需要找到一个 x∈[2,n] 使得对于每组 {a,b} 都满足以下两个条件之一：

• a，b<x
• a，b≥x

例如，如果 Mirko 选了 {{1,2},{3,4}}，那么 x 可以等于 3。

如果 Slavko 找不到满足条件的 x 值，则表示 Mirko 获得胜利。现在请你求出 Mirko 获胜的不同情况的总数，在对 10^9 取模后告诉他。

输入格式

第一行包含一个整数 N。

输出格式

第一行输出一个整数，为 Mirko 获胜的不同情况的总数对 10^9 取模后的值。

样例1

样例输入

2

样例输出

1

样例2

样例输入

3

样例输出

5

样例3

样例输入

4

样例输出

21

提示

样例 1 解释

Slavko 只有一种取法 {{1,2}}。

样例 2 解释

Slavko 的其中一种取法为 {{1,2},{1,3}}。

数据范围

对于 100% 的数据，1≤N≤20。

题解

最大公因数

求gcd（最大公因数）是这道题的前置知识点，我们一般用辗转相除法来求解。

证明如下：

若 x 与 y 的最大公因数为 k ，则 |x-y| 中也可被 k 整除。因此，我们可以让两数相减，再取两较小数继续相减，直到其中一个数为 0 为止，此时另一个数就是最大公因数。

我们可以进行类推，直接用 x 与 y 取模的时候得到的余数与 x 和 y 中的较小数继续取模，也可得到最大公因数，并且效率更高。

动态规划

我先来简单翻译一下题面。

就是当我们把一组里的数看做坐标轴上的一个线段，那么我们需要找出一些线段的选法，使得这些线段能够覆盖整个 1 到 n 的区间。

而且数据范围也不太大，所以显然，我们可以用dp来做这道题。（而且几乎不需要加优化）

我们定义 f[i][j] 表示在前 i 个线段中，能覆盖 1 到 j 区间的选法。

那么转移的之后就需要分三种情况：

1. 不选当前的线段
2. 选当前的线段并且可以和之前的连接起来覆盖更大的区间
3. 选当前的线段但不能和之前的连接起来

使用dp的时候需要注意后效性的问题，所以我们可以先O（n^2）地按右端点的顺序预处理出右端点在 n 及其左边的所有线段，后面就可以直接按顺序用了。

然后我们再来依次看上面的三种情况。

1. 不选当前的线段：直接从上一条线段覆盖相同长度转移。f[i][j]=(f[i][j]+f[i-1][j])%mo
2. 选当前的线段并且可以和之前的连接起来覆盖更大的区间：更新覆盖更大的区间的选法。f[i][a[i].r]=(f[i][a[i].r]+f[i-1][j])%mo
3. 选当前的线段但不能和之前的连接起来：还是更新当前覆盖的长度。f[i][j]=(f[i][j]+f[i-1][j])%mo

代码实现

//洛谷 P7800 [COCI2015-2016#6] PAROVI
#include<iostream>
#include<cstdio>
using namespace std;
const long long mo=1000000000;
int n;
long long f[400][25];//注意数组的大小
long long x;
int cnt;

struct edge{
	int l,r;
}a[400];

int gcd(int a,int b){//求最大公因数
	if(b==0){
		return a;
	}
	return gcd(b,a%b);
}

int main(){
	scanf("%d",&n);
	if(n==2){
		printf("1\n");
	}
	else{
		for(int i=1;i<=n;++i){
			for(int j=1;j<i;++j){
				if(gcd(i,j)==1){
					a[++cnt].l=j;
					a[cnt].r=i;
				}
			}
		}
		f[1][1]=1;
		f[1][2]=1;
		for(int i=2;i<=cnt;++i){
			for(int j=1;j<=n;++j){
				f[i][j]=(f[i][j]+f[i-1][j])%mo;//不选
				if(a[i].l<=j){
					f[i][a[i].r]=(f[i][a[i].r]+f[i-1][j])%mo;//选且可以连接
				}
				else{
					f[i][j]=(f[i][j]+f[i-1][j])%mo;//选但不可连接
				}
			}
		}
		printf("%lld\n",f[cnt][n]);
	}
	return 0;
}