青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

misschuer

常用鏈接

統計

積分與排名

百事通

最新評論

hdu 1402 A * B Problem Plus

http://acm.hdu.edu.cn/showproblem.php?pid=1402
#include<iostream>
#include
<cmath>
using namespace std;
typedef 
struct vir{
        
double re,im;
        vir()
{}
        vir(
double a,double b){re=a;im=b;}
        vir 
operator +(const vir &b)return vir(re+b.re,im+b.im);}
        vir 
operator -(const vir &b)return vir(re-b.re,im-b.im);}
        vir 
operator *(const vir &b)return vir(re*b.re-im*b.im,re*b.im+b.re*im);}
}
vir;
vir x1[
200005],x2[200005];
const double Pi = acos(-1.0);
void change(vir *x,int len,int loglen)
{
        
int i,j,k,t;
        
for(i=0;i<len;i++)
        
{
                t 
= i;
                
for(j=k=0;j<loglen;j++,t>>=1)
                        k 
= (k<<1)|(t&1);
                
if(k<i)
                
{
                        vir wt 
=  x[k];
                        x[k] 
= x[i];
                        x[i] 
= wt;
                }

        }

}

void fft(vir *x,int len,int loglen)
{
        
int i,j,t,s,e;
        change(x,len,loglen);
        t 
= 1;
        
for(i=0;i<loglen;i++,t<<=1)
        
{
                s 
= 0;
                e 
= s + t;
                
while(s<len)
                
{
                        vir a,b,wo(cos(Pi
/t),sin(Pi/t)),wn(1,0);
                        
for(j=s;j<s+t;j++)
                        
{
                                a 
= x[j];
                                b 
= x[j+t]*wn;
                                x[j] 
= a + b;
                                x[j
+t] = a - b;
                                wn 
=wn*wo;
                        }

                        s 
= e+t;
                        e 
= s+t;
                }

        }

}


void dit_fft(vir *x,int len,int loglen)
{
        
int i,j,s,e,t=1<<loglen;
        
for(i=0;i<loglen;i++)
        
{
                t
>>=1;
                s
=0;
                e
=s+t;
                
while(s<len)
                
{
                        vir a,b,wn(
1,0),wo(cos(Pi/t),-sin(Pi/t));
                        
for(j=s;j<s+t;j++)
                        
{
                                a 
= x[j]+x[j+t];
                                b 
= (x[j]-x[j+t])*wn;
                                x[j] 
= a;
                                x[j
+t] = b;
                                wn 
= wn*wo;
                        }

                        s 
= e+t;
                        e 
= s+t;
                }

        }

        change(x,len,loglen);
        
for(i=0;i<len;i++)
                x[i].re
/=len;
}



int main()
{
        
char a[100005],b[100005];
        
int i,len1,len2,t,over,len,loglen;
        
        
while(scanf("%s%s",a,b)!=EOF)
        
{
                len1 
= strlen(a)<<1;
                len2 
= strlen(b)<<1;
                len 
= 1;
                loglen 
= 0;
                
while(len<len1)
                
{
                        len
<<=1;
                        loglen
++;
                }

                
while(len<len2)
                
{
                        len
<<=1;
                        loglen
++;
                }

                
for(i=0;a[i]!='\0';i++)
                
{
                        x1[i].re 
= a[i]-'0';
                        x1[i].im 
= 0;
                }

                
for(;i<len;i++)
                        x1[i].re 
= x1[i].im = 0;
                
for(i=0;b[i]!='\0';i++)
                
{
                        x2[i].re 
= b[i]-'0';
                        x2[i].im 
= 0;
                }

                
for(;i<len;i++)
                        x2[i].re 
= x2[i].im = 0;
                fft(x1,len,loglen);
                fft(x2,len,loglen);
                
for(i=0;i<len;i++)
                        x1[i] 
= x1[i]*x2[i];
                dit_fft(x1,len,loglen);
                
for(i=(len1+len2)/2-2,over=loglen=0;i>=0;i--)
                
{
                        t 
= x1[i].re + over + 0.5;
                        a[loglen
++= t%10;
                        over 
=  t/10;
                }

                
while(over)
                
{
                        a[loglen
++= over%10;
                        over 
/= 10;
                }

                
for(loglen--;loglen>=0&&!a[loglen];loglen--);
                
if(loglen<0)
                        putchar(
'0');
                
else
                        
for(;loglen>=0;loglen--)
                                putchar(a[loglen]
+'0');
                putchar(
'\n');
        }

        
return 0;
}

posted on 2009-04-18 15:59 此最相思 閱讀(1907) 評論(6)  編輯 收藏 引用

評論

# re: hdu 1402 A * B Problem Plus 2009-10-13 20:18 zhou

這題真難呀!!!!  回復  更多評論   

# re: hdu 1402 A * B Problem Plus 2009-10-13 20:20 啊是

嗯  回復  更多評論   

# re: hdu 1402 A * B Problem Plus 2009-10-13 20:35 zhou

回復這么快!!! ORZ ORZ!!!  回復  更多評論   

# re: hdu 1402 A * B Problem Plus 2009-10-13 20:36 啊是

連接到Q了呀  回復  更多評論   

# re: hdu 1402 A * B Problem Plus 2010-01-16 00:43 abilitytao

請問有什么這方面的資料嗎?我最近也想研究這個問題 希望能和你交流^_^
我的QQ是 64076241  回復  更多評論   

# re: hdu 1402 A * B Problem Plus 2011-04-03 23:13 xiaomai

這題,太惡心了  回復  更多評論   


只有注冊用戶登錄后才能發表評論。
網站導航: 博客園   IT新聞   BlogJava   博問   Chat2DB   管理


青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            久久精品在线视频| 亚洲视频在线免费观看| 久久影视三级福利片| 一区二区三区四区国产| 久久久久久噜噜噜久久久精品| 亚洲一区二区三区在线观看视频 | 久久综合给合| 久久久91精品国产| 国产自产精品| 中文网丁香综合网| 欧美一级理论片| 韩国一区二区三区美女美女秀| 久久精品国产一区二区三区| 国产色视频一区| 亚洲视频在线观看| 久久―日本道色综合久久| 国产一区二区三区成人欧美日韩在线观看 | 欧美色大人视频| 欧美成人一区二区三区片免费| 国产欧美日韩视频一区二区| 亚洲看片免费| 亚洲天堂av综合网| 国产精品久久久久久久久久免费| 夜久久久久久| 久久精品夜色噜噜亚洲a∨| 国产欧美精品久久| 久久久久久噜噜噜久久久精品| 亚洲无线视频| 欧美日韩国产二区| 亚洲视屏在线播放| 亚洲一区二区成人| 国产精品人成在线观看免费 | 亚洲午夜影视影院在线观看| 国产欧美一区二区精品婷婷 | 亚洲在线观看| 老色鬼久久亚洲一区二区| 欧美激情按摩| 欧美一级视频| 在线观看久久av| 国产精品久久国产愉拍| 久久综合久久88| 午夜视频在线观看一区二区| 国产精品一区二区三区久久| 午夜精品久久久久久久久| 欧美aⅴ99久久黑人专区| 一本大道久久a久久综合婷婷| 国产精品区一区二区三区| 一区二区三区回区在观看免费视频| 久久精品国产99精品国产亚洲性色 | 亚洲婷婷综合色高清在线| 黄色av成人| 久久久蜜桃一区二区人| 亚洲欧美一区二区在线观看| 亚洲少妇最新在线视频| 亚洲国产精品一区| 久久综合伊人77777蜜臀| 亚洲欧美日韩国产一区二区三区 | 国内精品久久久久影院 日本资源| 欧美了一区在线观看| 亚洲肉体裸体xxxx137| 亚洲国产黄色| 最新日韩在线| 亚洲国产成人高清精品| 精品99一区二区| 精品成人在线视频| 黄色精品一二区| 日韩小视频在线观看专区| 欧美成人精品一区二区| 看片网站欧美日韩| 欧美刺激午夜性久久久久久久| 欧美一区二区三区四区在线观看| 亚洲性xxxx| 亚洲欧美激情精品一区二区| 亚洲欧美www| 一区二区三区福利| 亚洲欧美国产精品桃花| 欧美在线一二三| 久久久综合网| 91久久久久久久久久久久久| 亚洲精品视频啊美女在线直播| 一本大道久久a久久精二百| 午夜精品久久久久99热蜜桃导演| 久久国产精品99国产| 欧美激情黄色片| 亚洲欧美日韩精品久久亚洲区| 久久免费精品日本久久中文字幕| 久久免费黄色| 国产亚洲成人一区| 亚洲日韩视频| 久久美女艺术照精彩视频福利播放| 亚洲国产91精品在线观看| 亚洲欧美国产精品专区久久| 亚洲第一色中文字幕| 欧美自拍丝袜亚洲| 免费在线看一区| 欧美亚洲不卡| 国产欧美日韩综合一区在线播放| 亚洲国产视频一区| 久久久久久欧美| 亚洲男人av电影| 欧美风情在线| 99热免费精品| 亚洲毛片在线看| 欧美激情综合亚洲一二区| 91久久亚洲| 亚洲精品久久久久久久久| 欧美肥婆bbw| 欧美激情一区二区三区全黄| 亚洲国产成人在线视频| 久久久久久久久伊人| 国产区亚洲区欧美区| 欧美一区二区三区视频| 亚洲欧美日韩久久精品 | 亚洲欧美日韩精品久久| 国产一区二区三区成人欧美日韩在线观看 | 亚洲久久一区| 国产麻豆9l精品三级站| 久久中文欧美| 欧美黑人一区二区三区| 欧美一级二区| 亚洲与欧洲av电影| 国模一区二区三区| 亚洲国产婷婷香蕉久久久久久| 欧美高清视频免费观看| 亚洲图片欧美午夜| 久久www成人_看片免费不卡 | 欧美日韩一区成人| 先锋影音一区二区三区| 久久国产精品色婷婷| 亚洲视频二区| 欧美在线在线| 午夜久久久久久久久久一区二区| 久久久久久久一区| 欧美一区2区视频在线观看| 久久久噜噜噜| 久久精品人人做人人爽| 国产色产综合产在线视频| 午夜一区二区三区不卡视频| 久久疯狂做爰流白浆xx| 国产午夜精品一区二区三区视频| 亚洲在线第一页| 精品99一区二区| 亚洲尤物在线视频观看| 午夜精品三级视频福利| 国产精品美女久久久免费| 亚洲一区一卡| 美女国内精品自产拍在线播放| 激情欧美日韩| 欧美理论片在线观看| 亚洲精品乱码久久久久久日本蜜臀| 久久久久女教师免费一区| 亚洲二区视频| 久久se精品一区精品二区| 亚洲国产精品成人综合| 欧美日韩天堂| 久久综合久久久久88| 亚洲人成在线观看一区二区 | 欧美经典一区二区| 亚洲欧美日韩精品久久久久| 亚洲欧美日韩网| 亚洲国产日韩欧美综合久久| 欧美手机在线视频| 久久综合伊人77777蜜臀| 亚洲天堂偷拍| 亚洲三级电影在线观看| 麻豆久久久9性大片| 亚洲一区二区三区四区视频| 亚洲成人影音| 韩日欧美一区二区| 国产日韩欧美91| 国产精品欧美日韩久久| 欧美激情亚洲视频| 欧美国产精品久久| 久久久国产精品一区二区中文| 亚洲性视频h| 亚洲午夜av电影| 宅男在线国产精品| 亚洲性夜色噜噜噜7777| 亚洲视频在线免费观看| 亚洲深夜福利| 午夜激情亚洲| 欧美一区二区三区四区在线观看地址| 一区二区日韩精品| 中文一区在线| 男女激情久久| 亚洲欧美在线免费| 久久成人免费网| 老司机午夜精品视频在线观看| 欧美一区二区三区喷汁尤物| 欧美在线视频导航| 毛片精品免费在线观看| 亚洲黄色在线| 一本色道久久精品| 亚洲字幕在线观看| 欧美在线视频观看免费网站| 麻豆精品一区二区综合av| 欧美激情亚洲另类| 国产三区二区一区久久| 亚洲日本中文字幕免费在线不卡| 国产精品久久久久久久7电影|