题意： 给n个正整数，求有多少个GCD为1的子序列。答案对1e9+7取模。 1<=n<=1e5，数字ai满足1<=ai<=1e5 分析： 设f(x)表示以x为公约数的子序列个数 那么ans=f(1)- Σf(pi) + Σf(pi*pj) - Σf(pi*pj*pk) ........ 注意到对于一 ...
分类：
其他好文 时间：
2017-05-17 10:15:41
阅读次数：
174
# 线性代数# numpy.linalg模块包含线性代数的函数。使用这个模块，可以计算逆矩阵、求特征值、解线性方程组以及求解行列式等。 import numpy as np # 1. 计算逆矩阵# 创建矩阵A = np.mat("0 1 2;1 0 3;4 -3 8")print (A)#[[ 0 ...
分类：
其他好文 时间：
2017-05-10 16:32:50
阅读次数：
438
4176: Lucas的数论 Description 去年的Lucas非常喜欢数论题，但是一年以后的Lucas却不那么喜欢了。 在整理以前的试题时，发现了这样一道题目“求Sigma(f(i)),其中1<=i<=N”，其中 表示i的约数个数。他现在长大了，题目也变难了。 求如下表达式的值： 其中 表示 ...
分类：
其他好文 时间：
2017-05-09 18:55:21
阅读次数：
244
【BZOJ3781】小B的询问 题意：有一个序列，包含N个1~K之间的整数。他一共有M个询问，每个询问给定一个区间[L..R]，求Sigma(c(i)^2)的值,其中i的值从1到K，其中c(i)表示数字i在[L..R]中的重复次数 题解：初学莫队算法，差不多明白了用莫队的情况，对于这种离线的，区间长 ...
分类：
编程语言 时间：
2017-05-03 21:15:47
阅读次数：
253
题目链接 题意： n个点。依照题中给的公式能够求出随意两个点转移的概率。求从1到n的期望转移次数分析： 设dp[i]为从i到n的期望，那么能够得到公式dp[i] = sigma(dp[i + j] * p(i + j, i)) + 1。1 <= j <= m 把这个式子展开来：dp[i - m] * ...
分类：
其他好文 时间：
2017-04-22 18:50:02
阅读次数：
160
[BZOJ3781]小B的询问 试题描述 小B有一个序列，包含N个1~K之间的整数。他一共有M个询问，每个询问给定一个区间[L..R]，求Sigma(c(i)^2)的值,其中i的值从1到K，其中c(i)表示数字i在[L..R]中的重复次数。小B请你帮助他回答询问。 输入 第一行，三个整数N、M、K。 ...
分类：
其他好文 时间：
2017-04-22 10:28:17
阅读次数：
197
题目：给定n个数字a1...an。有m个询问，格式为L R X Y，意为求aL到aR之间与x的最大公因数为y的个数。 数据组数T<=20 1<=n,m<=1e5 1<=ai<=1e5 1<=L,R<=n;1<=X,Y<=1e5 分析： 考虑预处理出1~1e5所有数字的因子 然后就可以知道每个因子在1 ...
分类：
其他好文 时间：
2017-04-18 20:56:41
阅读次数：
175
2194: 快速傅立叶之二 Description 请计算C[k]=sigma(a[i]*b[i-k]) 其中 k < = i < n ，并且有 n < = 10 ^ 5。 a,b中的元素均为小于等于100的非负整数。 请计算C[k]=sigma(a[i]*b[i-k]) 其中 k < = i < ...
分类：
其他好文 时间：
2017-04-13 10:59:12
阅读次数：
195
$\Sigma_{i=1}^n\Sigma_{j=1}^mlcm(i,j)^{gcd(i,j)}$ $=\Sigma_{i=1}^n\Sigma_{j=1}^m (\frac{i*j}{gcd(i,j)})^{gcd(i,j)}$ 枚举gcd(i,j)=d $=\Sigma_{d=1}^n\Sigm ...
分类：
其他好文 时间：
2017-04-13 00:16:51
阅读次数：
238
\$w(t) \longrightarrow \bigg[\frac{\sqrt{2\sigma ^2\beta}}{s+\beta}\bigg] \longrightarrow \bigg[\frac{1}{s}\bigg] \longrightarrow y\$ $w(t) \longright ...
分类：
其他好文 时间：
2017-04-08 17:33:07
阅读次数：
141
