BZOJ1025: [SCOI2009]游戏

最新推荐文章于 2019-03-26 14:12:34 发布

DYT_B

最新推荐文章于 2019-03-26 14:12:34 发布

阅读量207

点赞数 1

CC 4.0 BY-SA版权

分类专栏：题解 BZOJ题解

本文链接：https://blog.csdn.net/DYT_B/article/details/79332288

题解同时被 2 个专栏收录

125 篇文章

订阅专栏

BZOJ题解

67 篇文章

订阅专栏

本文介绍了一种通过构建数学模型解决特定游戏循环结构的问题。针对给定数字范围内的所有可能排列，利用图论中环的概念，计算不同排列下形成特定循环结构的次数。通过将问题转化为求解一系列环的最小公倍数总数，采用素数筛法和动态规划方法给出了解决方案。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目描述：

Description
　　windy学会了一种游戏。对于1到N这N个数字，都有唯一且不同的1到N的数字与之对应。最开始windy把数字按顺序1，2，3，……，N写一排在纸上。然后再在这一排下面写上它们对应的数字。然后又在新的一排下面写上它们对应的数字。如此反复，直到序列再次变为1，2，3，……，N。如： 1 2 3 4 5 6 对应的关系为 1->2 2->3 3->1 4->5 5->4 6->6 windy的操作如下
1 2 3 4 5 6
2 3 1 5 4 6
3 1 2 4 5 6
1 2 3 5 4 6
2 3 1 4 5 6
3 1 2 5 4 6
1 2 3 4 5 6
这时，我们就有若干排1到N的排列，上例中有7排。现在windy想知道，对于所有可能的对应关系，有多少种可能的排数。
Input
包含一个整数N,1 <= N <= 1000

Output
包含一个整数，可能的排数。

Sample Input
【输入样例一】
3

【输入样例二】
10

Sample Output
【输出样例一】
3

【输出样例二】
16

题解：

这题的转化比较神= =。
首先我们发现所有的联系可以看做是一张图，把对应的点相连，最终形成一个个环，我们这时只要对每个环的大小取一个lcm就行了。但是这里有很多连的方式，所以答案就是所有lcm的数量总和。
我们把整个排列划分成多个循环节a1,a2,…ak，就可以表示为：k个数的最小公倍数的总数。这又可以表示为p1^a1*p2^a2*…*pk^ak是否为一个答案。
接下来我们就会发现一个结论：p1^a1+p2^a2+…+pk^ak<=n,则它是一个可行答案。
证明并不会。。。（详见dalao的blog）
那么我们就可以直接挖出素数来，再背包做一下就好了。
代码如下：

#include<cstdio>
#include<string>
using namespace std;
const int maxn=1005;
int n,m,p[maxn];
long long f[maxn][maxn];
bool vis[maxn];
void make_su(){
    for (int i=2;i<=n;i++){
        if (vis[i]==0) p[++m]=i;
        for (int j=1;j<=m&&p[j]*i<=n;j++) {
            vis[p[j]*i]=1;
            if (i%p[j]==0) break;
        }
    }
}
int main(){
    scanf("%d",&n);
    make_su();
    for (int i=0;i<=m;i++) f[i][0]=1ll;
    for (int j=0;j<=n;j++) f[0][j]=1ll;
    for (int i=1;i<=m;i++)
    for (int j=1;j<=n;j++) {
        f[i][j]=f[i-1][j];
        for (int k=p[i];k<=j;k*=p[i]) f[i][j]+=f[i-1][j-k];
    }
    printf("%lld\n",f[m][n]);
    return 0;
}