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

misschuer

常用鏈接

統(tǒng)計(jì)

積分與排名

百事通

最新評論

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

這題真難呀!!!!  回復(fù)  更多評論   

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

嗯  回復(fù)  更多評論   

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

回復(fù)這么快!!! ORZ ORZ!!!  回復(fù)  更多評論   

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

連接到Q了呀  回復(fù)  更多評論   

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

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

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

這題,太惡心了  回復(fù)  更多評論   


只有注冊用戶登錄后才能發(fā)表評論。
網(wǎng)站導(dǎo)航: 博客園   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>
            欧美电影电视剧在线观看| 亚洲日本久久| 久久久久久久久岛国免费| 亚洲欧美精品在线观看| 亚洲欧美在线视频观看| 性做久久久久久免费观看欧美| 一本不卡影院| 欧美一区三区二区在线观看| 久久久久久国产精品一区| 久久婷婷国产综合国色天香| 女人色偷偷aa久久天堂| 欧美日韩成人一区二区三区| 国产精品免费看片| 海角社区69精品视频| 在线播放精品| 亚洲小说春色综合另类电影| 久久精品论坛| 亚洲伦理精品| 久久国产欧美日韩精品| 欧美本精品男人aⅴ天堂| 国产精品xxxxx| 狠狠综合久久av一区二区老牛| 亚洲精品日韩精品| 欧美在线你懂的| 亚洲激情第一页| 一区二区三区日韩在线观看| 久久精品久久综合| 国产精品久久久久77777| 136国产福利精品导航| 亚洲欧美www| 亚洲国产天堂久久国产91| 亚洲一二三区在线观看| 乱码第一页成人| 国产伦精品一区二区三区视频孕妇| 亚洲国产精品99久久久久久久久| 亚洲欧美日韩人成在线播放| 亚洲丰满在线| 久久婷婷一区| 国产亚洲精品资源在线26u| 中文欧美在线视频| 欧美不卡视频| 久久久久久一区| 久久久久国产一区二区三区| 欧美国产日韩一区| 亚洲在线一区二区三区| 每日更新成人在线视频| 国产欧美日韩一区| 亚洲香蕉网站| 亚洲日本中文字幕区| 久久全球大尺度高清视频| 国产人久久人人人人爽| 亚洲在线一区二区三区| 亚洲欧洲在线观看| 欧美肥婆在线| 亚洲精品免费网站| 亚洲黄色在线视频| 欧美激情一区二区三区高清视频| 亚洲观看高清完整版在线观看| 久久午夜精品一区二区| 久久av资源网站| 黄色成人片子| 免费视频亚洲| 麻豆国产精品777777在线| 国内精品国语自产拍在线观看| 久久高清一区| 久久精品国产视频| 精久久久久久| 老司机免费视频一区二区三区| 香蕉久久夜色精品| 国产自产高清不卡| 久久深夜福利| 久久综合九色欧美综合狠狠| 1769国内精品视频在线播放| 女同性一区二区三区人了人一| 久久久精彩视频| 91久久精品网| 亚洲人成亚洲人成在线观看图片| 欧美精品尤物在线| 亚洲欧美日韩在线不卡| 欧美一区二区高清| 亚洲黄色免费电影| 亚洲免费激情| 国产丝袜一区二区三区| 欧美88av| 欧美先锋影音| 久久免费观看视频| 欧美国产第二页| 亚洲欧洲av一区二区| 久久国产精品久久国产精品| 亚洲国产你懂的| 亚洲无玛一区| 91久久久久久| 亚洲欧美999| 91久久精品网| 香蕉久久夜色精品国产使用方法| 1769国产精品| 亚洲一本大道在线| 亚洲激情成人网| 亚洲尤物视频网| 亚洲精品乱码久久久久久按摩观 | 性亚洲最疯狂xxxx高清| 欧美在线高清| 亚洲国产成人av| 亚洲视频电影图片偷拍一区| 国产日韩欧美麻豆| 亚洲国产一成人久久精品| 国产精品女人网站| 欧美激情免费在线| 国产日本欧美视频| 亚洲激情在线观看视频免费| 国产欧美大片| 日韩午夜在线播放| 一区二区在线观看视频| 中文亚洲字幕| 日韩一级大片| 免费成人av资源网| 久久久av毛片精品| 国产精品另类一区| 亚洲精品综合久久中文字幕| 在线观看av不卡| 性欧美18~19sex高清播放| 国产精品99久久久久久白浆小说| 久久久久久久网| 久久久av水蜜桃| 国产亚洲欧美一区| 亚洲欧美另类在线观看| 亚洲一区二区三区视频播放| 欧美成人精品一区| 欧美成年人视频| 伊人久久综合| 久久综合狠狠综合久久综青草| 久久米奇亚洲| 伊人影院久久| 麻豆精品视频在线观看视频| 狼人社综合社区| 黄色精品一区| 久久尤物视频| 亚洲欧洲日韩综合二区| 日韩视频在线观看免费| 欧美激情1区2区| 日韩视频中午一区| 午夜精品久久久久久| 国产精品日韩精品欧美在线 | 久久只有精品| 精品二区视频| 老司机午夜免费精品视频| 免费一级欧美片在线播放| 亚洲电影在线免费观看| 奶水喷射视频一区| 最新国产乱人伦偷精品免费网站| 日韩午夜视频在线观看| 欧美日韩免费观看一区| 妖精视频成人观看www| 午夜精品影院在线观看| 国产婷婷一区二区| 麻豆久久久9性大片| 亚洲第一精品夜夜躁人人爽| 99精品国产热久久91蜜凸| 欧美日韩大陆在线| 亚洲欧美成人一区二区在线电影| 久久福利一区| 亚洲精选一区| 国产老肥熟一区二区三区| 欧美一区二区三区啪啪| 欧美激情一级片一区二区| 一区二区激情小说| 国产精品一区二区三区四区五区 | 久久精品国产综合| 欧美国产日韩亚洲一区| 9l国产精品久久久久麻豆| 国产精品国产三级欧美二区 | 久久综合色8888| 亚洲精品综合久久中文字幕| 欧美在线播放| 亚洲免费激情| 国产日韩在线视频| 欧美精品久久一区| 午夜精品久久久久久久99热浪潮| 久久婷婷国产综合尤物精品| 日韩视频永久免费| 国产一区999| 欧美日韩在线免费视频| 久久精品一级爱片| 一本一本久久a久久精品综合麻豆 一本一本久久a久久精品牛牛影视 | 欧美激情成人在线| 香蕉乱码成人久久天堂爱免费| 亚洲第一精品在线| 国产女主播一区二区三区| 欧美国产综合| 久久嫩草精品久久久久| 亚洲欧美另类久久久精品2019| 最新69国产成人精品视频免费| 久久精品一区二区三区中文字幕| 亚洲免费播放| 亚洲激情网站免费观看| 国内免费精品永久在线视频| 欧美日韩免费在线视频| 欧美1区免费| 久久久久综合| 久久久国产亚洲精品|