D Interesting series


#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef double db;
typedef complex<db> C;
const int N = 1e5 + 5;
const ll mod = 100003;
const db pi = acos(-1);

int n, Q;
ll A, a[N];
vector<ll>w[N << 2];
ll fac[N], rev[N], ans[N];

ll qpow(ll a, ll b) {
	ll r = 1;
	while (b) {
		if (b & 1)r = r * a % mod;
		b >>= 1; a = a * a % mod;
	}
	return r;
}

void preC() {
	fac[0] = 1;
	for (int i = 1; i <= n; ++i)fac[i] = fac[i - 1] * i % mod;
	rev[n] = qpow(fac[n], mod - 2);
	for (int i = n - 1; ~i; --i)rev[i] = rev[i + 1] * (i + 1) % mod;
}

struct FFT {
	int n, m, rev[N << 2];
	C a[N << 2], b[N << 2];
	void init(int len) {
		for (n = 1, m = 0; n <= len; n <<= 1, m++);
		for (int i = 0; i < n; ++i) {
			rev[i] = rev[i >> 1] >> 1 | (i & 1) << (m - 1);
			a[i] = C(0, 0);
			b[i] = C(0, 0);
		}
	}
	void fft(C *a, int f) {
		for (int i = 0; i < n; ++i)if (i < rev[i])swap(a[i], a[rev[i]]);
		for (int i = 1; i < n; i <<= 1) {
			db alpha = pi / i;
			if (f == -1)alpha = -pi / i;
			for (int k = 0; k < i; ++k) {
				C w = exp(C(0, alpha * k));
				for (int j = k; j < n; j += (i << 1)) {
					C x = w * a[j + i];
					a[j + i] = a[j] - x;
					a[j] += x;
				}
			}
		}
		return;
	}
	void Calculate() {
		fft(a, 1); fft(b, 1);
		for (int i = 0; i < n; ++i)a[i] *= b[i];
		fft(a, -1);
	}
} F;

void Solve(int x, int l, int r) {
	if (l == r) {
		w[x].push_back(1);
		w[x].push_back(a[l]);
		return;
	}
	int mid = (l + r ) >> 1;
	Solve(x << 1, l, mid);
	Solve(x << 1 | 1, mid + 1, r);
	F.init(r - l + 1);
	for (int i = 0; i <= mid - l + 1; ++i)F.a[i] = C(w[x << 1][i], 0);
	for (int i = 0; i <= r - mid; ++i)F.b[i] = C(w[x << 1 | 1][i], 0);
	F.Calculate();
	for (int i = 0; i <= r - l + 1; ++i) {
		ll val = floor(F.a[i].real() / F.n + 0.5);
		w[x].push_back(val % mod);
	}
}

ll Com(int n, int m) {
	return n < m ? 0 : fac[n] * rev[m] % mod * rev[n - m] % mod;
}

int main() {
	scanf("%d%lld%d", &n, &A, &Q);
	for (int i = 1; i <= n; ++i) {
		scanf("%lld", &a[i]);
		a[i] = qpow(A, a[i] % (mod - 1));
	}
	Solve(1, 1, n);
	preC();
	for (int k = 1; k <= n; ++k) {
		ans[k] = (w[1][k] - Com(n, k) + mod) % mod * qpow(A - 1, mod - 2) % mod;
	}
	for (int i = 1, k; i <= Q; ++i) {
		scanf("%d", &k);
		printf("%lld\n", ans[k]);
	}
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
辽B代驾管理系统对代驾订单管理、用户咨询管理、代驾订单评价管理、代驾订单投诉管理、字典管理、论坛管理、公告管理、新闻信息管理、司机管理、用户管理、管理员管理等进行集中化处理。经过前面自己查阅的网络知识,加上自己在学校课堂上学习的知识,决定开发系统选择小程序模式这种高效率的模式完成系统功能开发。这种模式让操作员基于浏览器的方式进行网站访问,采用的主流的Java语言这种面向对象的语言进行辽B代驾管理系统程序的开发,在数据库的选择上面,选择功能强大的Mysql数据库进行数据的存放操作。辽B代驾管理系统的开发让用户查看代驾订单信息变得容易,让管理员高效管理代驾订单信息。 辽B代驾管理系统具有管理员角色,用户角色,这几个操作权限。 辽B代驾管理系统针对管理员设置的功能有:添加并管理各种类型信息,管理用户账户信息,管理代驾订单信息,管理公告信息等内容。 辽B代驾管理系统针对用户设置的功能有:查看并修改个人信息,查看代驾订单信息,查看公告信息等内容。 辽B代驾管理系统针对管理员设置的功能有:添加并管理各种类型信息,管理用户账户信息,管理代驾订单信息,管理公告信息等内容。 辽B代驾管理系统针对用户设置的功能有:查看并修改个人信息,查看代驾订单信息,查看公告信息等内容。 系统登录功能是程序必不可少的功能,在登录页面必填的数据有两项,一项就是账号,另一项数据就是密码,当管理员正确填写并提交这二者数据之后,管理员就可以进入系统后台功能操作区。项目管理页面提供的功能操作有:查看代驾订单,删除代驾订单操作,新增代驾订单操作,修改代驾订单操作。公告信息管理页面提供的功能操作有:新增公告,修改公告,删除公告操作。公告类型管理页面显示所有公告类型,在此页面既可以让管理员添加新的公告信息类型,也能对已有的公告类型信息执行编辑更新,失效的公告类型信息也能让管理员快速删除。新闻管理页面,此页面提供给管理员的功能有:新增新闻,修改新闻,删除新闻。新闻类型管理页面,此页面提供给管理员的功能有:新增新闻类型,修改新闻类型,删除新闻类型。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值