ljpadamaa-CSDN博客

原创 c++里面计算时间差

#include clock_t start, finish; start=clock(); finish=clock(); double duration = (double)(finish - start) / CLOCKS_PER_SEC; printf( "%f seconds\n", duration );

2015-04-19 16:11:23 1049

原创 qtcreator 里面使用openmp,c++11

QMAKE_CXXFLAGS += -std=c++11QMAKE_CXXFLAGS += -fopenmp -O2LIBS+= -lgomp -lpthreadadd these lines to the .pro file

2015-04-19 16:04:55 1136

原创 qtcreator里使用opencv 的各种坑

INCLUDEPATH += D:\software\opencv\mingw\install\include\opencv\D:\software\opencv\mingw\install\include\opencv2\D:\software\opencv\mingw\install\includeLIBS+=D:\software\opencv\mingw\lib\lib

2015-04-19 16:01:36 701

原创在windows下编译opencv

cmake的版本要3.0以上，否则出问题，环境变量Path要加入qt里的desianger.exe的目录。

2015-04-19 15:54:41 198

原创 poj 3735 Training little cats 构造矩阵

g i 让第i只猫得到一个豆e i 让第i只猫吃掉所有的豆s i j 让第i，j只猫交换豆给出n只猫，k条指令，所有指令重复m次因为是反复重复，所以可以用矩阵乘法，构造矩阵，利用矩阵的快速幂因为有加法（得到一个豆），所以得构造（n+1）*（n+1）的矩阵，最后一阶为1注意：矩阵乘法时，矩阵内有很多0，所有可以判断一下，在决定是否计算，否者会超时#include#

2013-10-30 21:19:01 249

原创 poj 1026 Cipher 置换（未完）

#include#include#includeusing namespace std;int num[300],n,pos[300],len[300];void circle(){ for(int i=1;i<=n;i++) { int j=i; len[i]=1; while(i!=pos[j]) {

2013-10-23 22:25:37 231

原创 poj 3270 Cow Sorting 置换问题

有n头牛，没头牛都有自己的暴躁指数，要把这些牛按暴躁指数从小到大排列，每次交换两头牛，消耗的体力为两头牛的暴躁指数和。应用的知识为置换。如把 1 2 3 4 5 6 变成 6 3 4 2 1 5 发现可以分成两组来交换，第一组是（6 1 5）第二组为（2，3，4）只需交换组内的元素就可以完成变换。所以第一步就先分组，在这道题中采用的方法是把元素排序，然后和原来的序列对比。那么每组内

2013-10-23 19:22:20 247

原创 poj 2154 Color polya 定理的优化

#includeusing namespace std;int mod;int eulr(int n){ int ans=n; for(int i=2;i*i<=n;i++) { if(n%i==0) { ans=ans/i*(i-1); while(n%i==0) n/=i;

2013-09-26 21:34:25 246

原创 hdu 1812 Count the Tetris polya 定理 java 大数

已知c种颜色和n＊n格的正方形，问能形成多少种不同的正方形旋转只有 0，90，180，270度三种旋法。旋0度，则置换的轮换数为n*n旋90度，n为偶数时，则置换的轮换数为n*n/4，n为奇数，则置换的轮换数为(n*n-1)/4+1旋180度，n为偶数时，则置换的轮换数为n*n/2，n为奇数，则置换的轮换数为(n*n-1)/2+1旋270度，n为偶数时，则置换的轮换数为n*n

2013-09-26 19:59:06 358

原创 poj 2409 Let it Bead polya 定理

#includeusing namespace std;long long poww(long long a,long long b){ long long ans=1; while(b--) ans*=a; return ans;}int gcd(int a,int b){ int r=0; while(b!=0) {

2013-09-26 18:59:32 319

原创 poj 1286 Necklace of Beads polya定

裸polya定理，已知有3种颜色的珠子，问用m个珠子能形成多少种项链（可旋转，可对称）polya定理：ans=(m^c1+m^c2+m^c3.....)/|G|m为颜色数，c为每种置换的循环节数，|G|为置换的个数具体到项链：项链顺时针旋转i个，循环节数为gcd(n,i) 翻转时，n为奇数，有n个循环节为（n+1）／2的循环群

2013-09-26 18:43:49 343

原创 poj 2661 Factstone Benchmark 算指数

题意：1960年可以用4位存储，以后每10年翻一倍，已知年份，问那年可以存储的最大n！是多少？算出最大的ln（n！）c语言算log（）用#include 中的log（n） c语言中 log(n) 相当与ln（n）所以这道题应转换成log（n）／log（2）计算#include#includeusing namespace std;int main(){

2013-09-08 20:34:54 368

原创 hdu 1521 排列组合指数型母函数

遇到的第一道指数型母函数，因为指数型母函数中系数是小数，所以用doulbe存，但是不知道是不是所有的题用都不了都能过。。。因为到最后还得转换成整数题意：共有n个数，a1重复n1次，a2重复n2次，a3重复n3次。。。问从这n个数中取r个数的排列数是多少#include#include#includeusing namespace std;double a[20],b[20],q

2013-09-08 19:41:59 397

原创 poj 1014 Dividing 母函数

题意：每行输入6个数，第i个数代表第i种珠宝的个数，第i种珠宝的价值为i问：能否把这些珠宝按价值平均分，也就是问这些珠宝是否能组合成总价值的1／2我是按母函数做的，也可以按多重背包做如果不对数据处理也可以直接过，不过很费时（900多毫秒。。。。。）还是应该对数据做一下处理，假设某种珠宝非常多，其实可以把它的数量减少到一定的量，但是依然不影响问题的答案，我们要做的就是找到减少到多少合

2013-09-08 17:17:23 416

原创 poj 1430 Binary Stirling Numbers 判断第二类stirling数的奇偶性

s（n,k）%2=c( n-m, n-m+(m-1)/2 )%2c（a,b）=a!/b!*(a-b)!分别判断a!,b!,(a-b)!中2的个数，然后算出c（a，b）中2的个数，如果大于0，就是偶数，等于0是奇数a！中2的个数=a/2+a/4+a/8.....(a／2计算出1，2...a中有多少个偶数，a／4计算出有多少个4的倍数，2的倍数有1个二，4的倍数有2个二)#includ

2013-09-05 21:26:57 468

原创 hdu 1711 Big Event in HDU 母函数

应该能用背包解，母函数有点慢用n种物品，已知每种物品的价值和数量问怎样把这些物品分成两堆，才能使两堆的价值差最小很好的母函数题，要不是归类到母函数里面，肯定想不到用母函数解法：先用母函数得出这些物品能组成的不同价值数（这里只要x^i的系数不为0，就能组合出价值i，不用考虑用系数＊i，因为这种方案会在后面的x^j里面包含了），然后找到和总价值的一半最接近的那个价值即可#inclu

2013-09-05 19:52:32 371

原创 hdu 1085 Holding Bin-Laden Captive! 母函数

拉登。。。。已知面值1，2，5这三种硬币的个数。问不能组成的最小面值是多少#include#includeusing namespace std;int a[9001],b[9001];int main(){ int n[3],w[3]; w[0]=1; w[1]=2; w[2]=5; while(cin>>n[0]>>n[1]

2013-09-05 19:26:15 335

原创 hdu 1398 Square Coins 母函数

有17种硬币每种硬币的价值为 i＊i，问n可以拆分成多少种硬币的组合#include#includeusing namespace std;int a[301],b[301];int main(){ int n; while(cin>>n) { if(n==0) break; memset(b,0,sizeof(b));

2013-09-05 19:05:29 326

原创 hdu 1028 Ignatius and the Princess III 整数的拆分母函数

问n可以拆分成多少种整数的和的形式 (5=2+3 5=3+2 算一种）有数论的做法更快，这里介绍的是母函数的方法#include#includeusing namespace std;int a[121],b[121];//x^i 的系数int main(){ int n; while(cin>>n) { memset(b,0,sizeo

2013-09-05 16:58:12 518

原创 hdu 2082 找单词母函数

基础的母函数已知每个字母的价值和个数，问这些字母组成的单词的价值小于50的有多少个#include#includeusing namespace std;int a[51],b[51];//记录x^0 x^1 x^2 ....的系数int main(){ int n; cin>>n; while(n--) { memset(a,0,siz

2013-09-05 16:38:49 358

原创 zoj 3113 John 反sg函数

和zoj2507相同题意：取石子游戏，不能取的人为胜对于任意一个Anti-SG游戏，如果我们规定当局面中所有的单一游戏的SG值为0时，游戏结束，则先手必胜当且仅当：（1）游戏的SG函数不为0且游戏中某个单一游戏的SG函数大于1；（2）游戏的SG函数为0且游戏中没有单一游戏的SG函数大于1。题中 sg（x）=x #includeusing namespace s

2013-08-23 14:05:27 314

原创 zoj 2507 Let's play a game Anti-SG游戏

题意：取石子游戏，不能取的人为胜对于任意一个Anti-SG游戏，如果我们规定当局面中所有的单一游戏的SG值为0时，游戏结束，则先手必胜当且仅当：（1）游戏的SG函数不为0且游戏中某个单一游戏的SG函数大于1；（2）游戏的SG函数为0且游戏中没有单一游戏的SG函数大于1。题中 sg（x）=x #includeusing namespace std;int main(){

2013-08-22 11:33:28 319

原创 zoj 3084 S-Nim 博弈

题意：还是取石子游戏，不过取的个数有限制用sg函数的定义求出每堆的sg函数然后^#include#include#includeusing namespace std;int sg[10010], w[10010];int s[200];int k;void get(){ memset(sg,0,sizeof(sg)); for(int i=1;i<10

2013-08-21 16:02:47 257

原创 zoj 3057 Beans Game 博弈直接推理

题意：有三堆石子每次可以从一堆里取，也可以从两堆里取同样多的，不能取的输不能用sg函数，只能用输赢的定义推注意数组不能开成int 否则超内存#include#includeusing namespace std;bool dp[301][301][301];void ok(int a,int b,int c){ if(a<=300&&b<=300&&c<=300)

2013-08-21 15:09:01 270

原创 hdu 4686 构造矩阵乘法待更新

这题给我坑坏了构造矩阵乘法，时刻要注意%mod 还有n==0时输出0#include#include#includeusing namespace std;const long long maxn=10;long long n,m,k;const long long mod=1000000007;typedef struct{ long long m[maxn][

2013-08-20 20:15:22 288

原创 zoj 2083 Win the Game 博弈 sg函数

题意：n个木条，每次在可以选一根在上面涂上2米长的颜色，不能涂的人输输入：木条的根数和长度思路：可以看成取石子，不过一堆石子在取完后变成两堆，这样求sg(x)时，就应用sg函数的定义sg（x）=mex{S} S为x的次态的sg的函数的集合 mex表示不在集合中的第一个数这道题中，每一个堆减少两个元素后变成两个堆，求出这两个堆的所有不同情况，然后用这两个堆的sg函数相互^ ,然

2013-08-17 20:30:24 344

原创 zoj 1827 The Game of 31 博弈 dfs

题意：有点数为1～6的扑克，每种点数有4张，每次每个人抽出一张，放到已知的牌堆中，牌堆中点数的总和不能超过31，无法出牌的人输输入：牌堆中已经有的牌没用sg函数用的还是博弈的基本结论：1.如果当前状态的次态有一个必败态，则当前状态为必胜态 2.如果当前状态的全部次态都为必胜态，则当前状态为

2013-08-17 13:56:19 309

原创 zoj 2290 博弈斐波那契数列

题意：取石子游戏，第二个人取的石子用dfs（）打表发现必败的状态满足斐波那契数列求第一步最少取几个石子：不是 n-最近的斐波那契数，因为这样不是最小的，发现第一阶段应取完 n-最近的斐波那契数这么多的石子，这就成了一个递归的问题，出口为 n-最近的斐波那契数==一个斐波那契数（

2013-08-17 11:30:48 371

原创 hdu 4678 Mine 博弈 sg函数

#include#include#include#include#includeusing namespace std;int tu[1010][1010];int aa[8]={1,1,1,0,0,-1,-1,-1},bb[8]={1,-1,0,1,-1,1,-1,0};int n,m,k;int ans,ff,ansk,ansf;queue q;vectorv;int

2013-08-16 14:40:12 369

原创 poj 3233 Matrix Power Series 矩阵快速幂二分

计算 s(k)= a+a^2+a^3+.....a^ka^k 用矩阵的快速幂因为k比较大，不能直接循环相加，因吃应用二分的方法s(k)=s(k/2)+s(k/2)*a^k/2 (k为奇，偶略有不同)#include#includeusing namespace std;const int maxn=40;int n,m,k;typedef struct{ i

2013-08-14 15:12:54 278

原创 hdu 4666 Hyperspace

可以添加和删除点，问这些点中两点间的最大曼哈顿距离曼哈顿距离：(x0,y0,z0) (x1,y1,z1) d=|x0-x1|+|y0-y1|+|z0-z1|去掉绝对值符号，整理后发现 d=(+-x0+-y0+-z0)-(+-x1+-y1+-z1) ＊＊＊前后括号中对应的符号相同＊＊＊因此只需记录相应符号下的最大，最小值即可应用set ，set中的元素默认是按从小

2013-08-14 13:31:39 285

原创 hdu 4671 Backup Plan

题意：有n个服务器，m个数据库，每个服务器可以给多个数据库服务，记录每个服务器服务的数据库的个数为ai ，|ai-aj|我是用贪心做的，只需考虑每个数据库的前两个服务器具体过程看代码 #include#include#includeusing namespace std;int a[120][120],b[120][120];//a[i][j]表示第i个数据库的第j个服务器，b[

2013-08-13 19:17:53 418

原创 poj 2649 素因子分解 ------ n！的素因子

n!的素因子分解中的素数p的个数为n/p+n/p^2+n/p^3+n/p^4+......原理：n/p 计算出了1～n中含有素因子p的个数 n/p^2计算出了1～n中含有因子p^2的个数 ......... 求和就是n！中素数p的个数注意：n或m为0时，都不能整除

2013-08-09 14:43:30 343

原创 poj 2689计算素数间距离

利用两次素数筛选，第一次筛选 1--500000，第二次利用第一次的结果筛选出l--r间的素数，注意如果l是1,，则筛选不出，应单独判断。最后穷举最大最小距离#include#includeusing namespace std;const int N=50001;bool prime[N+1];long long p[N+1],pri[1000005];int w;void g

2013-08-08 10:15:04 263