xiefubao的专栏

题意：给   出gcd(x,y,z)=G,lcm(x,y,z)=L,求x，y，z有多少种组合（考虑顺序）。解法：首先，x，y，z都是G的倍数，L/G则会分配给他们来承担。将L/G分解质因数，例质因子a出现了n次，首先在x，y，z中一定有一个或两个没有a，而且一个或两个有n个a。将其分成两类：1、一个有n个a，一个有0个a，另一个1-n-1 ；2、一个有n个a，其余两个均为有0个a； 3、两

题意：给出n个数n (1 ≤ n ≤ 106)： x1, x2, ..., xn (2 ≤ xi ≤ 107).            然后m (1 ≤ m ≤ 50000)次查询， [Li, Ri] ；输出每次Li与Ri之间所有素数的P值之和。每个素数的p值的定义是这样的：P(x)表示在给出的n个数里面有多少个数是x的整数倍数。解法：在线性筛素数的时候，把每个素数的P值可以直

今天在北大POJ上做到一道题——青蛙的约会，连续研究了数个小时，提交了共有50多次，每次都是以wrong anwser或是Time Limit Exceeded被打了回来。后来没吃晚饭，当看到最后一次显示“Accepted”时，内心的焦躁瞬间舒展开来了。   题目原文：两只青蛙在网上相识了，它们聊得很开心，于是觉得很有必要见一面。它们很高兴地发现它们住在同一条纬度线上，于是它们约定各自朝西跳，直

题意：在n个数中，找出k个数使得他们的乘积最大。并按从大到小输出这k个数。解法：一道简单的数学题，就是把非负数和负数分开来取。如果k为奇数，则一定会取到最大的非负数。其余的按以下方法继续计算得到；将正负两个分开的数列排序，都是绝对值从大到小。然后选择两个数列中前两个乘积较大的那个数列的前两个数（此步骤迭代k/2次）。其余的，注意计算细节就好了。#include #include

关于楼教主的第一题，感觉最费时间的还是高精度的实现。实现一万位以上整数基本运算的万进制高精度模版就写了快一天（题目里用到大概也有几百位吧），只能说太容易出错了，得不断地调试修改。题目源地址是北大poj1737题，题目本身的组合数学的算法是参考某位大牛的文章，在此附上地址http://hi.baidu.com/accplaystation/item/0b3b4e1654a07548e75e06f9

题意：求组合数A(m,n)中最后一个非零位是多少；解法: 确实是一道有意思的好题。这里有份不错的题解，讲的很详细:http://blog.csdn.net/luojinping/article/details/7555293  确实不太好想。以下是我看懂后敲的代码：#include #include using namespace std;int n,m;int tool[4

题意：求for(int i=a;i!=b;i+=c,i%=(1解法：即求解C*x-(1C*x+K*y=b-a；如果g=gcd(C,K)不能被b-a整除，则说明无解。         用exgcd()求出一组C/g*x+K/g*y=1的解，然后两边乘上(b-a)/g将求出的x取最小正数输出。代码：/*************************************

题意：http://acdream.info/problem?pid=1112Problem DescriptionHere  is Alice and Bob again !Alice and Bob are playing a game. There are several numbers.First, Alice choose a number n.Then he c