Problem of Precision hdu2256

最新推荐文章于 2010-08-10 20:55:00 发布

huixisheng

最新推荐文章于 2010-08-10 20:55:00 发布

阅读量433

点赞数

分类专栏： ACM 矩阵二分文章标签： struct 算法

本文链接：https://blog.csdn.net/huixisheng/article/details/5405257

版权

ACM 矩阵二分专栏收录该内容

17 篇文章 0 订阅

订阅专栏

http://acm.hdu.edu.cn/showproblem.php?pid=2256 //运行15ms的时候出错了，哪位大牛帮我看下！

/* 这个矩阵的构造有点稍微不一样的 */ #include <iostream> #include <cstdio> using namespace std; const int mod = 1024; struct Mat { int a; int b; }; Mat init, unit; Mat mul(Mat init, Mat unit) { Mat r; r.a = (init.a * unit.a + 6 * init.b * unit.b) % mod; r.b = (init.a * unit.b + init.b * unit.a) % mod; return r; } Mat mal(int n) { while(n) { if(n & 1) { init = mul(init , unit); n--; } else { unit = mul(unit, unit); n >>= 1; } } return init; } int main() { int t; cin >> t; while(t--) { int n; cin >> n; init.a = 1; init.b = 0; unit.a = 5; unit.b = 2; Mat ans; ans = mal(n); printf("%d/n", (2 * ans.a - 1) % mod); } }

#include<iostream>//为何这样就过了。感觉两个算法没什么不同啊？？想不明白！！
using namespace std;
const int mod = 1024;

struct Mat
{
   int n;
   int m;
}kk,rr;

Mat Mul(Mat p,Mat q)
{
  Mat r;
  r.n = p.n*q.n + 6 * p.m*q.m;
  r.m = p.n * q.m + p.m*q.n;
  r.n %= mod;
  r.m %= mod;
  return r;
}

Mat mal(int n)
{
  Mat p1;
  p1 = kk;
  if(n == 1) return p1;
  else
  {
   if(n %2 ==0)
    {
        p1 = mal(n/2);
        return Mul(p1,p1);
    }
   else
    return Mul(mal(n-1),kk);
  }
}//这里我真的想不明白是为何？递归的可以，为何下面的为何不可以啊？？
  //while(n != 1)
// {
  // if(n %2 == 0)
   // {
     // p1 = Mul(p1,p1);
     // n /= 2;
    //}
  // else
   // {
   // rr = Mul(rr,kk);
   // n--;;
   // }
// }
  //return Mul(rr,p1);

//}

int main()
{
int n, t;
cin>>t;
while(t--)
{
  cin>>n;
  kk.n = 5;
  kk.m = 2;
  //rr.n = 1;
  //rr.m = 0;
  Mat answer;
  answer = mal(n);
  cout<<(2*answer.n - 1)%mod<<endl;

}
}

#include<iostream>
using namespace std;
const int mod = 1024;

struct Mat
{
   int n;
   int m;
}kk,rr;

Mat Mul(Mat p,Mat q)
{
  Mat r;
  r.n = p.n*q.n + 6 * p.m*q.m;
  r.m = p.n * q.m + p.m*q.n;
  r.n %= mod;
  r.m %= mod;
  return r;
}

Mat mal(int n)
{
  Mat p1;
  p1 = kk;
  while(n != 1)
  {
   if(n %2 == 0)
    {
       p1 = Mul(p1,p1);
       n /= 2;
    }
   else
    {
     rr = Mul(rr,kk);
     n--;;
    }
  }
  return Mul(p1,rr);

}

int main()
{
int n, t;
cin>>t;
while(t--)
{
  cin>>n;
  kk.n = 5;
  kk.m = 2;
  rr.n = 1;
  rr.m = 0;
  Mat answer;
  answer = mal(n);
  cout<<(2*answer.n - 1)%mod<<endl;

}
}