hdoj 4686 矩阵相乘

感觉最近刷题不是练算法的，而是找bug的，今天做的矩阵相乘，明明是long long 的，但是我却用int型的，再哪里弄了半天,一直tle，感觉越来越对自己写的代码没有信心了。昨天和钢牛做武大上的比赛，不仅没有作对一题，而且贡献了无数次wa，一个简单的背包问题，初始化的问题让我很纠结，就这样简单题一直做不对让我对自己很没有信心，一到比赛都不敢写题了，感觉自己的都是错的，这样的恶性循环很不好的，每次我总说我这样做的好的话 就可以提高很多，但是都没有提高，集训也要结束了，没有多少时间了，该是对自己下点决心了，昨天多校，我就做了那个最没技术含量又最苦逼的题，还wa了好多次，今天又在没有必要的部分弄了那么久，无语啦！下面贴上本本的代码
：

/////////////////////////////////////////////////////////////////////////
// File Name: 4686.cpp
// Author: wang
// mail:
// Created Time: 2013-8-21 9:31:25
/////////////////////////////////////////////////////////////////////////
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <cstdio>
#include <cstdlib>
#include <climits>
#include <cstring>
#include <cmath>

#include <algorithm>
#include<iostream>
#include<queue>
#include <map>
using namespace std;
typedef long long ll;
#define INF (INT_MAX/10)
#define SQR(x) ((x)*(x))
#define rep(i, n) for (int i=0; i<(n); ++i)
#define repf(i, a, b) for (int i=(a); i<=(b); ++i)
#define repd(i, a, b) for (int i=(a); i>=(b); --i)
#define clr(ar,val) memset(ar, val, sizeof(ar))
#define mod 1000000007
#define M 5
#define N
ll n;
struct node
{
ll mp[M][M];
node()
{
rep(i,M) rep(j,M) mp[i][j]=0;
}
node(ll mp[M][M])
{
rep(i,M) rep(j,M) this->mp[i][j]=mp[i][j];
}
node operator * (node &a)
{
node r;
rep(i,M)
rep(j,M)
{
r.mp[i][j]=0;
rep(k,M)
{
r.mp[i][j]+=mp[i][k]*a.mp[k][j]%mod;
r.mp[i][j]%=mod;
}
}
return r;
}
};

node mul(node r,ll n)
{
if(n==1) return r;
node p;
rep(i,M) p.mp[i][i]=1;
while(n!=1)
{
if(n%2!=0) p=p*r;
r=r*r;
n/=2;
}
return r*p;
}
int main()
{
ll ax,ay,bx,by,a0,b0;
while(cin>>n)
{
cin>>a0>>ax>>ay;
cin>>b0>>bx>>by;
if(n==0)
{
cout<<0<<endl;
continue;
}
if(n==1)
{
cout<<a0*b0%mod<<endl;
continue;
}
ll a[M][M];
rep(i,M) rep(j,M) a[i][j]=0;
a[0][0]=ax; a[0][3]=ay;
a[1][1]=bx; a[1][3]=by;
a[2][0]=ax*by%mod; a[2][1]=bx*ay%mod;
a[2][4]=(ax*bx)%mod; a[2][3]=ay*by%mod;
a[2][2]=1;
a[3][3]=1;
a[4][0]=ax*by%mod; a[4][1]=bx*ay%mod;
a[4][4]=ax*bx%mod; a[4][3]=ay*by%mod;
//a0 b0 sum0 1为另一个坐标的值的
node r(a);
r=mul(r,n-1);
ll sum=0;
sum=(r.mp[2][0]*a0%mod+r.mp[2][1]*b0%mod+r.mp[2][2]*a0%mod*b0%mod+r.mp[2][3]%mod+r.mp[2][4]*a0%mod*b0%mod)%mod;
cout<<sum<<endl;
}
return 0;
}