【题解】[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;
}
相关推荐
试试勇气15 分钟前
Linux学习笔记(八)--环境变量与进程地址空间
linux·笔记·学习
蒙奇D索大17 分钟前
【数据结构】考研数据结构核心考点:平衡二叉树(AVL树)详解——平衡因子与4大旋转操作入门指南
数据结构·笔记·学习·考研·改行学it
怎么没有名字注册了啊19 分钟前
查找成绩(数组实现)
c++·算法
沐怡旸29 分钟前
【算法】725.分割链表--通俗讲解
算法·面试
郭庆汝1 小时前
自然语言处理笔记
笔记·自然语言处理·easyui
二进制怪兽1 小时前
[笔记] 驱动开发:Virtual-Display-Driver编译过程
笔记
ouliten1 小时前
cuda编程笔记(28)-- cudaMemcpyPeer 与 P2P 访问机制
笔记·cuda
im_AMBER2 小时前
数据结构 04 栈和队列
数据结构·笔记·学习
AI+程序员在路上2 小时前
QT6中Combo Box与Combo BoxFont 功能及用法
c++·qt
L_09072 小时前
【Algorithm】Day-4
c++·算法·leetcode