【题解】[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;
}
相关推荐
小龙报2 分钟前
《算法每日一题(1)--- 连续因子》
c语言·开发语言·c++·windows·git·算法·visual studio
祁同伟.9 分钟前
【C++】异常
开发语言·c++
第七序章24 分钟前
【C + +】红黑树:全面剖析与深度学习
c语言·开发语言·数据结构·c++·人工智能
rechol26 分钟前
汇编与底层编程笔记
汇编·arm开发·笔记
夜晚中的人海28 分钟前
【C++】滑动窗口算法习题
开发语言·c++·算法
AA陈超1 小时前
虚幻引擎5 GAS开发俯视角RPG游戏 P06-09 玩家等级与战斗接口
c++·游戏·ue5·游戏引擎·虚幻
violet-lz1 小时前
数据结构四大简单排序算法详解:直接插入排序、选择排序、基数排序和冒泡排序
数据结构·算法·排序算法
·白小白1 小时前
力扣(LeetCode) ——118.杨辉三角(C++)
c++·算法·leetcode
tongsound1 小时前
libmodbus 使用示例
linux·c++
lzj_pxxw1 小时前
嵌入式开发技巧:舍弃标志位,用宏定义函数实现程序单次运行
笔记·stm32·单片机·嵌入式硬件·学习