chenjunjie1994

2014年9月4日

摘要：求本原毕达三原组的个数以及不是毕达三元组的数。这道题需要把规模降下来，由x=m^2-n^2y=2mnz=m^2+n^2由有x+y#include #include #include #include using namespace std;const int Maxn=1000010;bool fl... 阅读全文

posted @ 2014-09-04 19:56 chenjunjie1994 阅读(112) 评论(0) 推荐(0) 编辑

POJ 1091

摘要：这题确实是好。其实是求x1*a1+x2*a2+....M*xn+1=1有解的条件。很明显，就是(a1,a2,...M)=1了。然后，可以想象，直接求有多少种，很难，所以，求出选择哪些数一起会不与M互质。。。好吧，思路就到这里了。。。T_T经过人提示，若（a1,a2,,,,an）与M不互质，则最大公约... 阅读全文

posted @ 2014-09-04 15:57 chenjunjie1994 阅读(267) 评论(0) 推荐(0) 编辑

POJ 2142

摘要：这道题的难点在于求|x|+|y|的为最小的值吧。想了好久才想出来，发现自己的数学能力确实跟不上。可知。x=x0+b/d*t;y=y0-a/d*t;则为|x0+b/d*t|+|y0-a/d*t|，仔细想想，可以看成的是两条直线方程y绝对值之和。那么，必然最小值只能出现在两条直线方程的两个零点之间，则枚... 阅读全文

posted @ 2014-09-04 11:33 chenjunjie1994 阅读(152) 评论(0) 推荐(0) 编辑

2014年9月3日

HDU 4686

摘要：再不能直视这道题，换INT64就过了。。。。。。。同样可以使用矩阵的方法。构造1*5的D[N],a[n],b[n],a【n】*b[n],1接着你应该就会了。#include #include #include #define LL __int64using namespace std;const L... 阅读全文

posted @ 2014-09-03 20:05 chenjunjie1994 阅读(160) 评论(0) 推荐(0) 编辑

HDU 4767

摘要：昨晚苦恼了一晚，因为即将大三了，必须要准备实习什么的事了。一般都会去公司实习吧，但是看看自己的简历，实在拿不出手，因为大一大二都在搞ACM（虽然真正搞的只有大二一年）,但却没有什么成绩，又不愿意做项目，因为只想一心搞ACM。可到了现在，如果真的没点什么的话，可能就去不了做实习生了。但是，想想自己，这... 阅读全文

posted @ 2014-09-03 16:30 chenjunjie1994 阅读(222) 评论(0) 推荐(0) 编辑

2014年9月2日

HDU 1757

摘要：简单的矩阵构造题，参看我前几篇的谈到的矩阵的构造法。#include #include #include #include using namespace std;int Mod;struct Matrax { int m[15][15];};Matrax a,per;int ats[15],an[... 阅读全文

posted @ 2014-09-02 16:41 chenjunjie1994 阅读(285) 评论(0) 推荐(0) 编辑

POJ 3613

摘要：可以利用DP的思想来做，不过是在DP时加上了矩阵乘法的思想而已，但乘法不是真的乘法，而是mp[a][i]+mp[i][b]I=0的话，则会一直不动，所以要初始化到I->I=INF;#include #include #include using namespace std;const int inf... 阅读全文

posted @ 2014-09-02 16:04 chenjunjie1994 阅读(206) 评论(0) 推荐(0) 编辑

HDU 2157

摘要：呃，离散数学课上直接就有这样的题#include #include #include using namespace std;const int Mod=1000;struct Matrax { int m[25][25];};Matrax a,per;int n,m;void initial(){... 阅读全文

posted @ 2014-09-02 13:36 chenjunjie1994 阅读(124) 评论(0) 推荐(0) 编辑

HDU 2371

摘要：知道了怎么置换之后，就可以用矩阵来置换了，但这道题一直关于置换的地方读不明白。#include #include #include using namespace std;const int Maxn=100;int pn[Maxn],xn[Maxn],bn[Maxn];int ansp[Maxn]... 阅读全文

posted @ 2014-09-02 10:18 chenjunjie1994 阅读(173) 评论(0) 推荐(0) 编辑

NYOJ 298

摘要：利用矩阵来做变换，参考Max大神的思想的，虽然不是同一道题。-----------给定n个点，m个操作，构造O(m+n)的算法输出m个操作后各点的位置。操作有平移、缩放、翻转和旋转这里的操作是对所有点同时进行的。其中翻转是以坐标轴为对称轴进行翻转（两种情况），旋转则以原点为中心。如果对每个点分别进行... 阅读全文

posted @ 2014-09-02 08:51 chenjunjie1994 阅读(207) 评论(0) 推荐(0) 编辑

2014年9月1日

HDU 3306

摘要：先转一些 http://www.cnblogs.com/frog112111/archive/2013/05/19/3087648.htmlFibonacci数列：F(0)=1 , F(1)=1 , F(n)=F(n-1)+F(n-2)我们以前快速求Fibonacci数列第n项的方法是构造常系数矩... 阅读全文

posted @ 2014-09-01 21:10 chenjunjie1994 阅读(194) 评论(0) 推荐(0) 编辑

HDU 1788

摘要：必须MARK下：任何时候都要保持清醒头脑，不要被题目绕了。。其实就是求最小公倍数。#include #include #include using namespace std;__int64 v[20];__int64 gcd(__int64 a,__int64 b){ if(b==0) retur... 阅读全文

posted @ 2014-09-01 11:02 chenjunjie1994 阅读(199) 评论(0) 推荐(0) 编辑

POJ 1006

摘要：直接使用中国剩余定理#include #include #include using namespace std;int r=3;int m[4]={0,23,28,33};int a[4]={0};void exgcd(int a,int b,int &d,int &x,int &y){ if(b... 阅读全文

posted @ 2014-09-01 10:25 chenjunjie1994 阅读(134) 评论(0) 推荐(0) 编辑

2014年8月31日

HDU 4365

摘要：把涂色的格子按对称旋转至左上角。当未涂色时，若要符合要求，则必须要求每一圈矩形都是上下左右对称的。注意是一圈的小矩形。对于N*N的阵，若最外层一圈的小矩形要符合要求，则（假设N%2==0）可以涂色的种数为K^（N/2）种。全个矩阵可涂色数为K^(N/2)*(N/2+1)/2。接第一段，(N/2)*(... 阅读全文

posted @ 2014-08-31 21:55 chenjunjie1994 阅读(179) 评论(0) 推荐(0) 编辑

HDU 3221

摘要：不知道哪里错了，测试了几十组数据均正确。。。可以找出规律，指数的增长是兔子数列。这个数列，是可以用矩阵快速幂得到的，见POJ 3070然后，竟然有一条公式：A^B%C = A^(B%phi(C)+phi(C))%C然后就可以求解了#include #include #include #include... 阅读全文

posted @ 2014-08-31 21:09 chenjunjie1994 阅读(194) 评论(0) 推荐(0) 编辑

勤能补拙是良训，一分辛苦一分才。