Good Bye 2018
A,B都是稍微想想就可以码出来。
主要c,d吧,规律题,还是不太熟悉
C New Year and the Sphere Transmission
题意:有一个n个元素的环,编号为1~n,一个球一开始在1,每次可以顺时针传k步(1 <= k <= n),每传到一个点上就会获得这个点的编号点的价值,直到球再次到1时停止,问对于所有k,能产生哪几种不同的价值。
思路:首先上来打了个表
然后思考发现,如果k是n不是约数,那么k的价值和1相同,如果k是n的约数,那么它的循环周期 t = n / k, 从1到第t个数分别为 1 , 1 + k, …, 1 + (t - 1)k, fun = (1 + 1 + (t - 1)k) * t / 2 ,等差数列求和。然后便有了分解质因数,排序统计答案的思路。总复杂度O(nlogn)
总结:环形规律就应该往循环节,最大公约数,约数方面去想,一开始想偏好久都
没回去,然后是范围n在int 内就开int ,尽量避免过多的long long 运算,会超时。
AC code:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 10;
ll ans[maxn];
// circle = k t = n / k
// fun = (1 + 1 + (t - 1) * k ) * t / 2
inline ll getfun(int k,int n){
int t = n / k;
return t + (ll)t * (t - 1) / 2 * k;
}
int main(){
int n;
cin >> n;
int tot = 0;
ans[tot++] = getfun(1,n); //因数1,和 n 特殊处理
int up = sqrt(n);
for(int i = 2; i <= up; ++i){
if(n % i == 0){
ans[tot++] = getfun(i,n);
if(i * i != n) ans[tot++] = getfun(n / i, n);
}
}
ans[tot++] = getfun(n,n);
sort(ans,ans + tot);
for(int i = 0; i < tot; ++i){
cout << ans[i] << ' ';
}
return 0;
}
D. New Year and the Permutation Concatenation
题意:给n!个n的排列,按字典序从小到大连成一条序列,例如3的情况为:[1,2,3, 1,3,2, 2,1,3 ,2,3,1 ,3,1,2 ,3,2,1],问其中长度为n,且和为sum=n*(n+1)/2的序列有多少个?
这题是队友最后扔给我的规律,我码过了,没认真想,第二天起来补题,看到了一位博主的博客,感谢这位博主良心翻译:https://www.cnblogs.com/pkgunboat/archive/2018/12/31/10201676.html
官方题解:我们考虑一下next_perumation函数产生字典序递增的全排列的过程:
假设某一个序列长度为n,最长的递减的后缀长度k,那么它的下一个排列是这样产生的:选取序列第n-k个数,与后k个数中最小的那个交换(其实就是最后一个),然后将后k个数按从小到大排序。
例如序列1,2,5,4,3的下一个排列为1,3,2,4,5。我们观察发现:这种时候1,2,(5,4,3,1,3,)2,4,5不满足和为sum了,因为在产生下一个排列的过程中,第n-k个位置的数被替换了。
也就是说,假设一个序列存在长度为k的递减后缀,那么这个后缀不能产生一个长度为sum的序列。例如,1,2,(5,4,3,1,3,)2,4,5不行,但是1,(2,5,4,3,1,)3,2,4,5可以。
所以,我们的任务是找出每个长度为k的递减后缀有多少个?应该为C(n,n-k)*(n-k)!=A(n,n-k)=n!/k!个。因为只要选了前面n-k个数,后面长度为k的递减的序列是固定的,所以我们只需要选n-k个数全排列就行了。
我们可以得到最终的答案了:一共有n*n!-(n-1)个序列,要减去( ∑(k from 1 to n-1) n!/k! )- (n-1)个。
为什么要减去n-1个呢?我们来看最后一个排列(假设n为5)5,4,3,2,1 。5之后的序列不存在,所以要从总的序列数中减去。而这(n-1)个不存在的序列恰好会被判定为不满足题意,也应该减去。
所以总的来说,答案应该是:(所有的序列-不存在的序列)-(不满足的序列-不存在的序列)即:
ans=n*n!-∑(k from 1 to n-1) n!/k!。
AC code:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 998244353;
const int maxn = 1e6 + 5;
int f[maxn],s[maxn];
int main(){
int n;
cin >> n;
f[0] = f[1] = 1, s[n] = 1;
for(int i = 2; i <= n; ++i){
f[i] = (ll)f[i-1] * i % mod;
}
// s[ ] n ! / k !
//从后往前逆推这样好,O(n)
for(int i = n - 1; i >= 0; --i){
s[i] = (ll)s[i + 1] * (i + 1) % mod;
}
int ans = (ll)n * f[n] % mod;
for(int i = 1; i <= n - 1; ++i){
ans = (ans - s[i] + mod) % mod;
}
cout << ans <<endl;
return 0;
}
//比赛时候的AC code
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod = 998244353;
int main()
{
int n;
cin >> n;
ll fac = 1, ans = 1;
if(n==1)
{
puts("1");
return 0;
}
for(int i = 2; i <= n; ++i)
{
fac = fac * (i-1) % mod;
ans = (ans + fac - 1) % mod;
ans = ans * i % mod;
}
cout << ans << endl;
return 0;
}
Finally, happy new year!