刚开学没多久就打了一个网络赛，通过这次网络赛我是发现我是真的菜...
放假前校赛的排名让我有些自满，寒假丝毫没有接触ACM，一直沉迷于Steam，这个真的值得好好反省。
虽然现在大一课有点多，在学校也有些事务，但是这些都不是我松懈的理由，
在此写下这篇博客就是为了提醒自己：Why（为什么别人进科协，我要打 ACM ），How，What
这次比赛的反思：
数论的学习实在是太过于薄弱，要加强，对数字的规律不够敏感，要锻炼，
数据结构最常用的树，不会，要学，
这次题目总体不算特别难，题目的灵活度不大，关键自己会的知识太少。
题目链接
http://acm.hdu.edu.cn/search.php?field=problem&key=%A1%B0%D7%D6%BD%DA%CC%F8%B6%AF-%CE%C4%D4%B6%D6%AA%D0%D0%B1%AD%A1%B1%B9%E3%B6%AB%B9%A4%D2%B5%B4%F3%D1%A7%B5%DA%CA%AE%CB%C4%BD%EC%B3%CC%D0%F2%C9%E8%BC%C6%BE%BA%C8%FC&source=1&searchmode=source
HUD 6461（概率的问题）（答案4行1.00，//题目有点没说明白）
HDU 6462（斐波那契数列平方和）
题意有点看不太清楚，但是所给的数据还是很容易发现规律，
斐波那契数列每一项的平方和，把所给的4进制a,b,换成10进制就行，打个表就OK了
#include <iostream>
#include <cstdio>
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)>(y)?(y):(x))
using namespace std;
const int mod=192600817;
const int maxn=4e4+5;
long long fio[maxn]={0,1};
long long sum[maxn]={0,1};
void make_table()
{
for(int i=2; i<maxn; i++)
{
fio[i]=(fio[i-1]+fio[i-2])%mod;
sum[i]=(fio[i]*fio[i])% mod+sum[i-1];
//后面是要求 sum的差值，这里的sum不要求余，可大致看一下也能发现long long 不会溢出
}
}
int main()
{
make_table();
/*
for(int i=0; i<10; i++)
{
printf("%d %I64d %I64d\n", i , fio[i], sum[i]);
}
*/
int Q;
while(cin>>Q)
{
while(Q--)
{
int a,b,c,d;
cin>>a>>b>>c>>d;
a=a*4+b+1;
c=c*4+d+1;
b=max(a,c);
d=min(a,c);
cout<<(sum[b]-sum[d-1])%mod<<endl;
//取余别忘记了
}
}
return 0;
}
HDU 6463（查询第k个鸽子数）
一看就是肯定有规律，
方法1：
随便取几个数字出来，打个表（我傻傻的那笔在那里找规律，把电脑当计算器??）看一下，都会出现那些重复的数字，
发现要么会出现4要么会出现1，然后打个表就OK了；
1 #include <iostream>
2 using namespace std;
3
4 const int maxn=150005;
5 long long table[maxn];
6
7 bool myJudge(int n)
8 {
9 while(n!=1&&n!=4)
10 {
11 int sum=0;
12 while(n)
13 {
14 sum+=(n%10)*(n%10);
15 n/=10;
16 }
17 n=sum;
18 }
19 return n==1;
20 }
21
22 void make_table()
23 {
24 int cnt=0;
25 for(int i=1; cnt<maxn; i++)
26 {
27 if(myJudge(i))
28 table[cnt++]=i;
29 }
30 }
31
32 int main()
33 {
34 make_table();
35 int Q; cin>>Q;
36 while(Q--)
37 {
38 int k; cin>>k;
39 cout<<table[k-1]<<endl;
40 }
41
42 return 0;
43 }
方法2：
当时感觉有循环节，确定一个大概的循环次数，超出这个次数还未有1即false，不过当时很快就找出来规律，这个就没写了
方法3：
记忆化搜索（配合set）
HDU 6464
看了不少题解都说要用 带权值的线段树，或者主席树，这个还不会，先放着，等学到这个在A了
HDU 6465（三角形的变换）
找初始三角形边的向量和变换后的边的向量，在确定一下矩阵的类型，求矩阵。
求出这个矩阵之后，对于每个给你的点，先求出它关于三角形第一个点变换之前的向量，乘上这个矩阵得到变换之后的向量，再加上三角形第一个点变换之后的坐标即可。
1 #include <iostream>
2 #include <cstdio>
3 using namespace std;
4
5 int main ()
6 {
7 int T; cin>>T;
8 while(T--)
9 {
10 double x[10],y[10];
11 for(int i=1;i<=6;i++)
12 cin>>x[i]>>y[i];
13
14 for(int i=2;i<=3;i++)
15 x[i]-=x[1],y[i]-=y[1];
16 for(int i=5;i<=6;i++)
17 x[i]-=x[4],y[i]-=y[4];
18
19 double a,b,c,d;
20 c=(x[5]*x[3]-x[6]*x[2])/(y[2]*x[3]-x[2]*y[3]);
21 a=(x[5]-y[2]*c)/x[2];
22 d=(y[5]*x[3]-y[6]*x[2])/(y[2]*x[3]-x[2]*y[3]);
23 b=(y[5]-y[2]*d)/x[2];
24
25 int Q; cin>>Q;
26 while(Q--)
27 {
28 cin>>x[2]>>y[2];
29 x[2]-=x[1],y[2]-=y[1];
30 printf("%.2f %.2f\n",x[2]*a+y[2]*c+x[4],x[2]*b+y[2]*d+y[4]);
31 }
32 }
33 return 0;
34 }
HDU 6466
超出的我的能力范围...
HDU 6467（数学式的简化+快速幂(这种式子化简的方法有待总结)）
F(n)=∑i=1n(i×∑j=inCij)
难点在于这种式子式如何化简的；
化简得F(n)=（n-1）*2^n+1;
#include <iostream>
#include <cstdio>
using namespace std;
long long mod=1000000007;
long long poww(long long a, long long b)
{
long long ans=1,base=a;
while(b)
{
if(b&1)
ans=(ans*base)%mod;
base=(base*base)%mod;
b>>=1;
}
return ans;
}
int main ()
{
long long n;
while(cin>>n)
{
cout<<(poww(2,n)*((n-1)%mod)+1)%mod<<endl;
}
return 0;
}
HDU 6468 (应该是线段树划分区间或者字典树 暂时还不会，有待解决)
HDU 6469 不会
HDU 6470 （给F（n）的递推式求F（n））
给了fn的递推式，确定一下如何用矩阵表示，矩阵快速幂即可。（下图引用自csdn https://blog.csdn.net/winter2121/article/details/88612870）
发现自己有必要学习如何定义矩阵（class）以及矩阵重载运算符这些，这样的话会方便很多。
#include <iostream>
#include <cstdio>
using namespace std;
const int mod=123456789;
struct Matrix{
long long m[6][6];
};
Matrix Mul(Matrix A, Matrix B)
{
Matrix ans;
for(int i=0; i<6; i++)
for(int j=0; j<6; j++)
ans.m[i][j]=0;
for(int i=0; i<6; i++)
for(int j=0; j<6; j++)
for(int k=0; k<6; k++)
ans.m[i][j]=((A.m[i][k]*B.m[k][j])%mod+ans.m[i][j])%mod;
return ans;
}
Matrix Quick_Pow(Matrix A, long long n)
{
Matrix ans;
for(int i=0; i<6; i++)
for(int j=0; j<6; j++)
{
if(i==j)
ans.m[i][j]=1;
else
ans.m[i][j]=0;
}
while(n)
{
if(n&1)
ans=Mul(ans,A);
A=Mul(A, A);
n>>=1;
}
return ans;
}
int main()
{
int T; cin>>T;
while(T--)
{
long long n; cin>>n;
Matrix A;
for(int i=0; i<6; i++)
for(int j=0; j<6; j++)
A.m[i][j]=0;
A.m[0][0]=1; A.m[0][1]=2; A.m[0][2]=1; A.m[0][3]=3; A.m[0][4]=3; A.m[0][5]=1;
A.m[1][0]=1;
A.m[2][2]=1;
A.m[3][2]=1; A.m[3][3]=1;
A.m[4][2]=1; A.m[4][3]=2; A.m[4][4]=1;
A.m[5][2]=1; A.m[5][3]=3, A.m[5][4]=3; A.m[5][5]=1;
Matrix ans=Quick_Pow(A, n-2);
long long Ans=0;
long long a[6]={2,1,1,2,4,8};
/*
for(int i=0; i<6; i++)
for(int j=0; j<6; j++)
{
printf("%I64d ",A.m[i][j]);
if(j==5)
cout<<endl;
}
cout<<endl;
*/
for(int i=0; i<6; i++)
{
Ans=(Ans+(a[i]*ans.m[0][i])%mod)%mod;
}
cout<<Ans<<endl;
/*
for(int i=0; i<6; i++)
for(int j=0; j<6; j++)
{
printf("%I64d ",ans.m[i][j]);
if(j==5)
cout<<endl;
}
*/
}
return 0;
}
原文地址：https://www.cnblogs.com/Yokel062/p/10577872.html