「BZOJ1025」[SCOI2009] 游戏
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。
Output
包含一个整数,可能的排数。
Sample Input
「输入样例一」
3
「输入样例二」
10
3
「输入样例二」
10
Sample Output
「输出样例一」
3
「输出样例二」
16
3
「输出样例二」
16
HINT
「数据规模和约定」
100%的数据,满足 1 <= N <= 1000 。
题解
wulala:
这道题可以转换一下。
试想每一个对应关系a-b为从a->b的一条边,
那么图中一定存在n条边且每个点入度出度都为1,
易证一定存在一个或几个环。
实际上排数就是这几个环大小的最小公倍数。
即求和为n的数列的最小公倍数种数。
那么可以直接DP
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 |
#include<iostream> #include<cstdio> #include<cstring> #include<cstdlib> #include<algorithm> #include<cmath> #include<set> #include<queue> #include<map> #define pa pair<int,int> #define inf 1000000000 #define ll long long using namespace std; inline int read() { int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } int n,tot,pri[1005]; ll ans,f[1005][1005]; bool mark[1005]; void getpri() { for(int i=2;i<=1000;i++) { if(!mark[i])pri[++tot]=i; for(int j=1;j<=tot&&i*pri[j]<=1000;j++) { mark[i*pri[j]]=1; if(i%pri[j]==0)break; } } } void solve() { f[0][0]=1; for(int i=1;i<=tot;i++) { for(int j=0;j<=n;j++)f[i][j]=f[i-1][j]; for(int j=pri[i];j<=n;j*=pri[i]) for(int k=0;k<=n-j;k++) f[i][k+j]+=f[i-1][k]; } for(int i=0;i<=n;i++)ans+=f[tot][i]; } int main() { n=read(); getpri(); solve(); printf("%lld",ans); return 0; } |
Subscribe