qingyezhu

摘要： /* * timus1748.c * * Created on: 2011-10-6 * Author: bjfuwangzhu *//* 这里面我们有一个prime【16】的数组，为什么只要这几个素数呢，因为这几个素数的乘积大于10^16, 而且就反素数的性质来说 比如2^t1*3^t2*5^t3*...p1^x***p2^y,假设p1是大于prime[]中所有的素数的， 因为这几个素数的乘积大于10^16,如果 我们添加p1在这个连乘积式子里面，那么必然有至少一个prime[i]不在这个连乘积式子里面， 但是对于因子的总数目而言 我们在乎的是幂的大小，而非素因子的大小，也就是说... 阅读全文

摘要： 以下转自于：http://blog.sina.com.cn/s/blog_59e67e2c0100a84o.html题目意思不难已知给定k,x,y求 1<=a<=x 1<=b<=y中满足gcd(a,b)=k 的(a,b)对数。（注意数对是无序的）。 1<=x,y<=10w, 0<=k<=10w题目有比较恶心的一点，数据有k==0的，这时显然答案是0，没有2个数的gcd为0。首先，gcd是没啥用的。因为约掉gcd后两个数互质。于是我们可以让x/=k y/=k并且假设 x<=y然后题目变成了 2个数分别在区间[1..x]和[1..y]中的互质数 阅读全文

摘要： 以下转载于：http://blog.csdn.net/xieshimao/article/details/6840731数论题！求与N不互质的数的K次方（K=4），反过来想若知道与N互质的K次方和，那所求就容易多了哦。观察到与n互质的数的性质比如12=2*2*3那么与12不互质的数就有2,3,4,6,8,9,10,12其实就是2的所有倍数，以及3的所有倍数所以可以先求一个1到12的所有数的四次方和。这个有公式：n*(n+1)*(2*n+1)*(3*n*n+3*n-1)/30注意对与除以30可以看成是乘以30的逆元（对于1000000007的逆元）。求的所有的四次方和之后当然要减去那些不互质的数 阅读全文

摘要： /* * sgu231.c * * Created on: 2011-10-1 * Author: bjfuwangzhu */#include<stdio.h>#include<string.h>#include<stdlib.h>#include<math.h>#define nmax 8169int num[nmax] = {5, 7, 13, 19, 31, 43, 61, 73, 103, 109, 139, 151, 181, 193, 199, 229, 241, 271, 283, 313, 349, 421, 433, 463, 阅读全文

摘要： import java.io.BufferedInputStream;import java.math.BigInteger;import java.util.Scanner;public class Main { public static void main(String[] args) { Scanner cin = new Scanner(new BufferedInputStream(System.in)); BigInteger num, two = BigInteger.valueOf(2), one = BigInteger.ONE; ... 阅读全文

摘要： /* * sgu154.c * * Created on: 2011-10-1 * Author: bjfuwangzhu */#include<stdio.h>#define nnum 5#define nmax 0x7fffffffint getNum(int n) { int res; res = 0; while (n) { res += n / nnum; n /= nnum; } return res;}void solve(int n) { int left, right, mid, temp; l... 阅读全文

摘要： pku 1804归并排序/* * poj1804.c * * Created on: 2011-9-29 * Author: bjfuwangzhu */#include<stdio.h>#define nmax 1010int num[nmax], temp[nmax];int res;void merge(int left, int mid, int right) { int i, j, k; i = left, j = mid + 1, k = 0; while (i <= mid && j <= right) { if (num[j] < 阅读全文

摘要： zjut1005 做幻方此题是求解奇幻方记：关键是确定下一个数填的位置，规律是：1以后的每一个数只能填在它所在位置的下一行的下一列（记为A），如果A已经有数，就把下一个数填在当前这个数的上面（即与这个数同一列的上一个位置）#include<stdio.h>#include<string.h>#define nmax 101int num[nmax][nmax];int main() {#ifndef ONLINE_JUDGE freopen("data.in", "r", stdin);#endif int n, nn, i, j 阅读全文

摘要： 扇形面积=R*R*Hudu(扇形弧度)/2;弧长L=R*Hudu(扇形弧度）/* * zoj1597.c * * Created on: 2011-9-21 * Author: bjfuwangzhu */#include<stdio.h>#include<math.h>#define pi acos(-1.0)#define eps 1.0e-8double distance(double x1, double y1, double x2, double y2) { return sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 阅读全文

摘要： 使用矩阵的方法解决zoj2105结果出现了一下两种情况：仅仅是改了一下求第n项，变为先求第n-1,n-2项，再求第n项！结果竟是不一样~~路过的大牛指导一下，谢谢！！Wrong的代码：/* * zoj2105.c * * Created on: 2011-9-20 * Author: bjfuwangzhu */#include<stdio.h>#define nmax 2#define nnum 7typedef struct matrix { int num[nmax][nmax];} matrix;matrix mul_matrix(matrix a, matrix b... 阅读全文