hdoj 5666 Segment (大数取模)(俄罗斯乘法)

最新推荐文章于 2019-11-20 20:06:31 发布

阿.荣.

最新推荐文章于 2019-11-20 20:06:31 发布

阅读量410

点赞数

分类专栏：数学大数

本文链接：https://blog.csdn.net/bmicnj/article/details/52143221

版权

数学同时被 2 个专栏收录

67 篇文章 0 订阅

订阅专栏

大数

9 篇文章 0 订阅

订阅专栏

Segment

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 1532 Accepted Submission(s): 564

Problem Description

Silen August does not like to talk with others.She like to find some interesting problems.

Today she finds an interesting problem.She finds a segment

x+y=q .The segment intersect the axis and produce a delta.She links some line between

(0,0) and the node on the segment whose coordinate are integers.

Please calculate how many nodes are in the delta and not on the segments,output answer mod P.

Input

First line has a number,T,means testcase number.

Then,each line has two integers q,P.

q is a prime number,and

2≤q≤1018,1≤P≤1018,1≤T≤10.

Output

Output 1 number to each testcase,answer mod P.

Sample Input

Sample Output

题意：x+y==p，直线与x，y坐标围成三角形中整数坐标的点的个数，不算直线上的点。

题解：推导公式：(p-1)*(p-2)/2 %p, 数据大，超long long ,有个俄罗斯乘法可以快速相乘，（要保证除以2 的数是奇数）

俄罗斯乘法是一种计算两数相乘的算法。

举例如下：

计算 35*72

过程

35 72

17 144

8 288

4 576

2 1152

1 2304

从上到下，对每一行，若左边的数字若为奇数，则将右边的数字取出，累加。

72+144+2304=2520

累加的结果2520即为乘积。

因为该算法不需要已知九九乘法表。

所以在计算机中有应用。

代码：

#include<cstdio>
#define LL long long
LL sum;
LL multi(LL a,LL b,LL p)
{
	sum=0;
	while(b)
	{
		if(b&1)
		{
			sum+=a;
			sum%=p;
		}
		a=a*2;
		a%=p;
		b>>=1;
	}
	return sum;
}
int main()
{
	int t;
	LL q,p,a,b;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%lld%lld",&q,&p);
		a=q-1;
		b=q-2;
		if(a%2)			//保证 除以2 的数是奇数 
			printf("%lld\n",multi(b/2,a,p));
		else
			printf("%lld\n",multi(a/2,b,p));
		
	}
	return 0;
}