A PERFECT NUMBER PROBLEM 简单数论
求因子和(除自身)等于自身的前五个数。
第五个数是3e8级别的,比赛的时候我用了奇技淫巧,但是这不是长久之计。
现在认真的做一下:
质因数分解 n=a1k1a2k2……aiki
f(n)=(1+a1+a12+……+a1k1)(1+a2+a22+……+a2k2)……(1+ai+ai2+……+aiki) - n
欧拉线性筛(讲解求约数和传送门
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
const int N = 4e6 + 100;
int prime[N], f[N], fu[N], tot;
bool is_prime[N];
void init() {
memset(is_prime, true , sizeof(is_prime));
for (int i = 2; i < N; i++) {
if (is_prime[i]) {
prime[tot++] = i;
f[i] = 1 + i;
fu[i] = i;
}
for (int j = 0; j < tot && i * prime[j] < N; j++) {
is_prime[i * prime[j]] = false;
if (i % prime[j] == 0) {
fu[i * prime[j]] = fu[i] * prime[j];
if (fu[i] == i) {
f[i * prime[j]] = f[i] + i * prime[j];
}else {
f[i * prime[j]] = f[i / fu[i]] * f[fu[i] * prime[j]];
}
break;
}else {
fu[i * prime[j]] = prime[j];
f[i * prime[j]] = f[i] * f[prime[j]];
}
}
if (f[i] == i * 2)
printf("%d\n", i);
}
}
int main() {
puts("6");
puts("28");
puts("496");
puts("8128");
puts("33550336");
return 0;
}
法二:
发现
n = 1+21+ 22 + 23 + 2z + (n/21)+(n/22)+…+n/(2z).
然后化简这个式子,就是n=2z*2z-1。然后枚举z就可以了,看看有没有满足的。
//标程
#include <stdio.h>
int main(){
int x[5] = {1, 2, 4, 6, 12};
int i;
for(i = 0; i < 5; i++) {
printf("%d\n", (1 << x[i]) * ((1 << x[i] + 1) - 1));
}
return 0; }
H Coloring Game 思维
求从左上到右下的相连方块的方案总数。
一开始把这题想成递推了,推出来才发现 发现对于1e9的数据根本出不来。
失误,以后做题要先看数据再下手
两端左上和右下为黑,其余可黑可白,对答案贡献是* 2 * 2.
中间不能为全白,故每列贡献为*3
最终答案= 4 * 3n-2
快速幂即可
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <cstdlib>
#include <string.h>
#include<iomanip>
#define M 1000000007
#define INF 10000000
#define LL long long
using namespace std;
int n;
LL quick_mod(LL a,LL x)
{
LL res=1;
while(x)
{
if(x&1) res=res*a%M;
a=a*a%M;
x>>=1;
}
return res;
int main(){
scanf("%d",&n);
if (n==1) printf("1\n");
else
{
LL sum = quick_mod(3,n-2);
sum = (sum * 4)%M;
printf("%lld\n", sum);
}
return 0;
}
I Max answer 单调栈+线段树
找区间最小值*区间和的最大值。
首先题目原型是一道单调栈的提目,给定正数,只需要以每个数为最小值寻找左右边界即可。
这题还加了负数进去,那就需要讨论了
M Subsequence 思维
给一个主串S,询问n次,每次给出一个串t,问t是否为S的子序列。
赛后发现一个一个字符暴力往后遍历有两组数据过不了:正确的算法是对字符串t,计算当前字母j在主串的i位置后的第一次出现的位置,最后判断i是否跳出主串,没跳出则t为其子串。
//暴力遍历
#include<bits/stdc++.h>
#define M 1000000007
#define N 100010
#define ll long long
using namespace std;
int main() {
char s[N],t[N];
int n,len;
scanf("%s",s);
scanf("%d",&n);
len = strlen(s);
while(n--)
{
scanf("%s",t);
int q=0,p=0,l = strlen(t);
while(p < len && q < l)
{
if(s[p] == t[q])
{
p++;
q++;
}
else p++;
}
if(q==l) cout<<"YES\n";
else cout<<"NO\n";
}
return 0;
}
用nxt[i][c]记录位置i后c第一次出现的位置。
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
const int N = 1e5 + 100;
const int M = 26;
int n, m, L;
int nex[N][M];
char str[N];
void init() {
memset(nex, 0x3f3f3f3f, sizeof(nex));
L = strlen(str);
for (int i = L - 1; i >= 0; i--) {
for (int j = 0; j < 26; j++) {
if (str[i] - 'a' == j) nex[i][j] = i;
else nex[i][j] = nex[i + 1][j];
}
}
}
int main() {
scanf("%s", str);
init();
scanf("%d", &n);
while (n--) {
scanf("%s", str);
int len = strlen(str);
int p = -1;
for (int i = 0; i < len && p < L; i++) {
p = nex[p + 1][str[i] - 'a'];
}
if (p < L) puts("YES");
else puts("NO");
}
return 0;
}