【骗经验】NOI 2013 Day2 T1题解
noip吧
全部回复
仅看楼主
level 13
轩轩醉了 楼主
两场比赛蒟蒻也就这题能AC了,现场还没AC掉,哭瞎了。。。
还是说解法吧
1.暴力模拟,20分
2.矩阵+高精度+快速幂,70~80分。
发现转移方式有两种,那么分开列矩阵。
A=|a b| \n |0 1|
B=|c d| \n |0 1|
C=|1| \n |1|
其中A是第一种递推的,B是第二种递推的C是最后乘的。
那么A^(n-1) 就是一行的转移,再乘B就是一整行的转移,记为X,那么X^(m-1) * C后就能得到结果。
然后要加高精度,可得70~80分。
3.矩阵+快速幂+费马小定理
(a,p)=1则a^(p-1)≡1(mod p)
那么将n对1000000006取模,再用2中说的快速幂,可得满分。
蒟蒻当时没想到费马小定理T_T
2013年07月18日 07点07分 1
level 11
@wwwaaannngggrs UCCU带了多么不好的头。。
2013年07月18日 07点07分 3
level 9
\n。。。。
2013年07月18日 07点07分 4
没办法[委屈]
2013年07月18日 07点07分
level 10
我是直接用通项公式+快速幂+费马小定理,不需要矩阵乘法的
2013年07月18日 07点07分 5
OTL通项
2013年07月18日 07点07分
OTL
2013年07月18日 08点07分
OTL
2013年07月18日 08点07分
OTL!不过通项特别好推为什么都不推?是因为学的太多数学水平下降了?
2013年07月19日 23点07分
level 11
利用这个
构造(f[i] + A) = a * (f[i-1] + A)
A = b / (a - 1) 构造一个等比..
然后利用费马小定理乱搞..
2013年07月18日 08点07分 7
OTL
2013年07月18日 08点07分
level 11
不需要高精,不需要矩阵= =...
利用构造等比数列 (f[i] + x) = a * (f[i-1] + x) -> x = b / (a - 1) 这个方法..
2013年07月18日 08点07分 8
回复 yangff2 :这个显然好算通项啊..
2013年07月18日 09点07分
回复 _喔嘞嘞 :看到递推果断去算矩阵了
2013年07月18日 09点07分
我就是这么做的。话说出题人方法是Jordan块的Fermat小定理
——来自 诺基亚 Lumia 928
2013年07月20日 12点07分
level 8
求通项
等比数列求和 + 费马小定理+欧几里得扩展求逆元
2013年07月18日 08点07分 9
还有快速幂
2013年07月18日 08点07分
level 8
直接矩阵快速幂就好分成十份而不是二分快速幂处理十进制数超爽[吐舌]
2013年07月18日 09点07分 10
level 11
可以推通项公式+费马小定理高精度取模没开long long送了50分结果就只有Cu了
2013年07月18日 09点07分 11
Cu都没有哭瞎了
2013年07月18日 09点07分
level 9
强势围观各路大神
2013年07月18日 10点07分 13
level 13
[鬼脸]蒟蒻表示矩阵最后优化成两个变量{a,b} 转移也是{x.a*y.a,x.b*y.a+y.b} 你还在用矩阵么
2013年07月19日 17点07分 14
ans = a + b
2013年07月19日 17点07分
优化下来确实是这么回事啊。不管咋说还是费马小定理靠谱
2013年07月20日 01点07分
回复 轩轩醉了 :我去构造一个题好了
2013年07月20日 01点07分
回复 帅芞の面包 :去吧面包大神。对了,在火车上思考下比赛题目啊。我准备放到28号。
2013年07月20日 01点07分
level 13
问!
为神马费马小对矩阵乘法成立啊[疑问]
考场上就是没推出这个,我就去推通项了……
2013年07月20日 01点07分 15
为什么不成立。。。
2013年07月20日 01点07分
回复 轩轩醉了 :费马小定理不是基于完全剩余系的吗?[疑问]
2013年07月20日 02点07分
回复 pyx1997 :互素就可以了啊。这里模的是素数啊
2013年07月20日 02点07分
回复 轩轩醉了 :您说的是不是对于矩阵A,A^m = A^(m%(P-1)) (mod P), (P in Prime),同余指矩阵的对应元素同余……请问矩阵和数的互质是怎么做到的……或者我语文能力又拙计了- -||
2013年07月20日 03点07分
level 13
横排O(1),竖排O(1)
2013年07月20日 04点07分 16
level 9
推通项不是很欢快?。
2013年07月20日 12点07分 17
level 10
直接快速幂+各种常数优化AC怎么破。。
2013年07月20日 12点07分 18
太神
2013年07月21日 01点07分
1