【题解】[GenshinOI Round 3] P9816 少项式复合幂

题目链接

分析

首先这题给了很大的提示信息 注意 m 和 p 的范围 , 很自然的想到可以先把所有可能的 f ( x ) f(x) f(x) 算出来.

思维误区

有些人在算完 f ( x ) f(x) f(x) 之后可能就会去思考找环的问题,然后一些码力弱的大佬就会祭掉.

在经过仔细的观察之后 (大多数人其实一眼就看出来了罢 , 可以发现最终答案的计算是符合结合律的,或者说具有传递性? 所以考虑倍增.

令 f a [ i ] [ j ] fa[i][j] fa[i][j] 表示 f 1 < < j ( i ) f_{1<<j}(i) f1<<j(i) 的值,初始时把 f [ i ] [ 0 ] f[i][0] f[i][0] 算出来,后面就可以直接倍增了.

Code

cpp 复制代码
#include <bits/stdc++.h>
#define int long long
const int N = 1e5+10;

using namespace std;
int m,q,p;
int ksm(int a, int b){
	int ans = 1;
	while(b){
		if(b&1){
			ans = ans * a % p;
		}
		a = a*a%p;
		b >>= 1;
	}
	return ans;
}
int a[30],b[30];
int f[N];
int get(int x){
	int ans = 0;
	for(int i = 1; i <= m; i++){
		ans = (ans + a[i]*ksm(x,b[i])%p) % p;
	}
	return ans;
}
bool vis[N];
int belong[N];
vector<int> e[N];
int fa[N][30]; 
void init(){
	for(int i = 0; i < p; i++){
		fa[i][0] = get(i);
	}
	for(int i = 1;i <= 25; i++){
		for(int j = 0; j < p; j++){
			fa[j][i] = fa[fa[j][i-1]][i-1];
		}
	}
}
signed main(){
	cin >> m >> q >> p;
	for(int i = 1; i<= m; i++){
		cin >> a[i] >> b[i];
		a[i] %= p;
	} 
	init();
	while(q--){
		int x,y;
		cin >> x >> y;
		x %= p;
		for(int i = 25; i >= 0; i--){
			if((1 << i) <= y) x = fa[x][i],y -= (1<<i);
		}
		cout << x << endl;
	}
	return 0;
}
相关推荐
巴伦是只猫几秒前
【机器学习笔记Ⅰ】3 代价函数
人工智能·笔记·机器学习
ZZZS05162 分钟前
stack栈练习
c++·笔记·学习·算法·动态规划
位东风25 分钟前
【c++学习记录】状态模式,实现一个登陆功能
c++·学习·状态模式
hans汉斯28 分钟前
【人工智能与机器人研究】基于力传感器坐标系预标定的重力补偿算法
人工智能·算法·机器人·信号处理·深度神经网络
AI视觉网奇2 小时前
rag学习笔记
笔记·学习
vortex52 小时前
算法设计与分析:分治、动态规划与贪心算法的异同与选择
算法·贪心算法·动态规划
前端拿破轮2 小时前
🤡🤡🤡面试官:就你这还每天刷leetcode?连四数相加和四数之和都分不清!
算法·leetcode·面试
雷羿 LexChien2 小时前
C++内存泄漏排查
开发语言·c++
地平线开发者3 小时前
征程 6|工具链量化简介与代码实操
算法·自动驾驶
嘉小华3 小时前
CMake 完全指南:第一章 - 构建的烦恼 - 为什么需要CMake?
c++